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1. Introduction 


In classical continuum mechanics, the motion of a material medium relative 
to a rectangular Cartesian system of spatial gbordsriat es Bi U1 Os Oy Mise hee 


presented by a set of three single-valued functions Zee ('T), K=1, 2, 3, of the 
K 


coordinates z* and the time T. The functions Z representing a real continuous 
motion are subject to the condition 


K 
Dadeta,Z+0, a=. (1.1) 


eee K 
Every set of three values Z* = Z(z', T) for the functions Z serves as names for 
the material points P of the continuous medium*, Every set of values for the 
four space-time coordinates (z', T) serves as names for an event E. 


Kae 
* We shall use the convention that a set of functions like Z (z’, T) will be written 
with the labeling index set over the kernel letter. When we wish to regard the values 
of such a set of functions as coordinates we write ZX. Thus, for example, the partial 


K Kon : 
derivatives of the functions Z are again functions 0; 2 (z*, T) of the coordinates 2’, T. 
However, from the condition (1.1), it follows that we can solve for the z* as functions 
of the ZX and T. We write the solution in the form #=2(Z¥*, T). The gradients 
6; Z* can be regarded as Senor of the Z¥ and T, in which case we denote these func- 
tions by 0,2*(Z*, T)=9; Pe (z AZK, T), T). This convention avoids oe RENe use of 


new symbols and gives a precise meaning to the two different quantities 0; Vs and 0; 2%. 
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The functions Z representing the motion specify the material point Z* which 
experiences the event (?', T). We restrict attention ae what follows to local 


properties of motions. We assume that the functions ZS are analytic at an event 
E with coordinates es T). Thus we assume the existence of a power series 
0 o 60 


Zé, 1) =Z (7, T) tZ, (2 


C : : ‘ ; K 
A) +Z,(T-1)+42Z,,(¢—#) a) +--4 240, (1.2) 


convergent at every event E in some neighborhood N (E) of the event E. In 
all that follows, by ‘‘analytic’’ we shall mean analytic at the event E and all 


functions of the coordinates that appear will be assumed analytic at BE. From 
(1.1) follows the existence of a unique inverse 


A= 2(ZK,T). (1.3) 


Ke 
The functions z : will be analytic at the point with coordinate values 2 = Jb ( z, T), 
Di Tr. 


us Ai(T) be any set of nine analytic functions of T such that 
OAL AL = 6". (1.4) 
A‘(T) is an orthogonal matrix. Let d'(T) be any three analytic functions of 


J 
Ue 
the time T. Then in terms of a given motion Z(z', 7) we define the class of all 


Kes 
motions Z ‘(z', T) which differ from Z(z',T) by a rigid motion as follows: 
Kee Ke ee eee 
Z'(#,T) =Z(Aj2'+d', T + constant). (1.5) 


Conversely, any two motions (Z’, 4 gy ») between which there exists a relation of 
the form (1.5) are said to differ from each other by a rigid motion. 


An important problem in classical continuum mechanics is the construction 
of constitutive equations for the stress, internal energy, and heat flux. In general, 
these constitutive equations depend on a motion and are required to transform 
in a definite way when a motion is replaced by any motion differing from it by 
a rigid motion. For example, the theory of finite deformations of elastic media 
is based on constitutive equations for the stress tensor ¢7 having the general 
form [12, 8, 11] 


aS 


. ij M 
vi(2 T) = T (Z (zt, yn Skeet (1.6) 
where, for definiteness, we may ea we the functions ‘ are peta and 


single valued in all 12 arguments Z and ¢; he The functional form of T depends 
on the ae properties of the material medium. However, for all elastic 


materials the fe are required to satisfy the invariance condition 


ij K M ki K 


T(2!, OZ!) = Ay ACh e(Zae wat (4.7) 
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aes 
where (Z Pes ) is any pair of motions differing from each other by a rigid motion. 
Using known methods of invariant theory one can prove that any set of func- 


aj 
tions T satisfying (1.7) is reducible to the form 


a 


=; 


ie 


pe! Mal Riu PO ep es ON 
Pease (ZsGaD), atliGes 04 O:7,.0,Z. 4 (1.8) 


In all that precedes, we have regarded a motion as being relative to a fixed 
rectangular Cartesian system of spatial coordinates z’. However, the invariant- 
theoretic problems of continuum mechanics involving the class of all motions 
differing from each other by a rigid motion can be attacked from another point 
of view which is more convenient for our present purposes. Consider the group 
of coordinate transformations 


f= At (Te ba" (7), 


AG 
aie COnsranite ~ 


i : : Lae, 
where A; (7) is an orthogonal matrix. Let Z(z’, T) represent a motion relative 


; Le as 
to the frame (2,7). Let Z’(2’,T’) be the scalar transform of Z(z', T). Thus 
% Ie 
L423 FAT) (1.10) 


according to the definition of a scalar in tensor calculus. The functions Zz ee da) 
define a motion relative to the frame (z’, 7’). The transformation law (1.7) 
implies that 

ie Cees a Wa oe We Os bao CI I 4 Cale (1.11) 


If we define a set of 16 quantites t””, w,v=1, 2,3, 4 by setting 


LID 
TY = | ; ee (1.12) 


we see that the law of transformation (1.11) for ¢‘/ is precisely the law of trans- 
formation implied by transforming the quantities t”” as a 4-dimensional tensor 
under the group of coordinate transformations (1.9) in four variables, where the 


time T is regarded as the 4" coordinate. Thus, transforming the functions 2 
representing a motion as a set of three scalars under the group (1.9), we see 
that the fundamental assumption (1.11) of continuum mechanics assumes a 
simple and familiar form in terms of a 4-dimensional tensor law of transformation 
under the group (1.9). 

In presenting these familiar ideas we have attempted to indicate the im- 
portant role in mechanics of the theory of invariants of a motion under a group 
of coordinate transformations in a 4-dimensional space. In this work we shall 
introduce three such groups of transformations: the Euclidean group, the Galilean 
group, and the Lorentzian group. Each of these groups of coordinate trans- 
formations on 4 real variables defines a type of geometric space. In each of these 
spaces, we shall define a motion of a continuous medium. We then define 
Euclidean kinematics, Galilean kinematics, and Lorentzian kinematics as the 


Be 
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theory of the invariants of a motion in these 4-dimensional spaces. We have 
attempted to design a formalism which treats these three kinematical theories 
on a parallel basis. This has been done so that a comparison of classical kine- 
matical concepts, definitions, and theorems with their relativistic counterparts 
is made easier. Another advantage of this symmetrical treatment is that we can 
borrow ideas from the more familiar classical invariant theory of a motion and 
transfer them by analogy to the relativistic case. The Euclidean and Galilean 
groups are firmly interlocked with classical mechanics; the Lorentzian group, 
with electromagnetic theory and relativistic mechanics. However, kinematics 
being the science of motion in itself, independent of the natural laws presumed 
to govern the motion, we have considered here only briefly the application of 
kinematical results to these more restricted theories. 


2. Euclidean, Galilean, and Lorentzian space-time 


So as to fix the meaning we attach to the words space, geometry, field, and 
invariant, allow us to describe briefly how they shall be used. 

An n-dimensional geometric space is a set of points p such that to every 
point p there corresponds a subset of points X¢(p) containing p which can be 
placed into one to one correspondence with all the ordered sets of 7 real numbers 
it 07, ..., &°) Lyme im somevinterval 4" = pr xP RE, h“>0 and such 


0 
that p corresponds to aM, together with a group & of allowable coordinate trans- 


formations x =x (x"), x=x(x). The class of all coordinate systems related 
by elements of the group & is called the class of allowable coordinate systems 
for the points 9t(p) in the neighborhood of p. The characterization of a space 
may involve also a set of functions ®,, ®,..., ®y of the coordinates x“ having 
an assigned transformation law under the group &. By choosing different groups 
® and different sets of functions ® and different transformation laws for the 
set ®, we obtain various examples of geometric spaces. Thus ordinary 3-dimen- 
sional Euclidean metric space corresponds to letting & be the group of orthogonal 
transformations, where with this choice of &, the set of functions ® is empty. 
However, we can also let & be the group of general analytic transformations ®, 
provided we append a Euclidean metric tensor field g, ;(«*). The two spaces so 
defined are regarded as equivalent. Curved Riemannian spaces, affinely con- 
nected spaces, conformal spaces, efc., correspond to various other choices for the 
group & and the functions ® together with their transformation law [8, 10]. 
The foregoing example of Euclidean metric space shows that different choices 
of the pair of objects (6, ®) may serve to define spaces which are regarded as 
equivalent. 

By a field ® in a space Y with group & we shall mean a set of functions of 
the coordinates Dp (x"), Q=1, 2,..., N having an assigned law of transformation 
under the group &. By “an assigned law of transformation” we mean a rule 
such that when the representation of the field ® by functions @.(x") in any 
one allowable coordinate system x" is given, the representation Oy (x) of 
the field in any other allowable coordinate system is uniquely determined 


by the functions ®,(x“) and the coordinate transformation relating the x“ 
and 52". 
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A tensor field in a space Y with group © is a set of functions of the coordinates 
having a transformation law under & of the general form 


Oy (0) (2 x)? (sen ey EAS a4 DH (x), 
ee eee 6x Ox” ex? a eee (2.1) 
(x' x) = det 
Ox” 


If y =0, ® is called a tensor field of weight w. If y=1, ® is called an axial tensor 
field of weight w. If w=y=0, ® is called an absolute tensor field. The number 
of superscripts and subscripts on a tensor field determine its contravariant and 
covariant rank, respectively. The rank of a tensor field is the sum of its contra- 
variant and covariant ranks. Tensors of rank zero are called scalars, and tensors 
of rank one are called vectors. 


An affine connection is a field T having a law of transformation of the form [10] 


; xt ax! ax” Ox ax” 
CR ra EE apltiesietes = re: D2 
cae Ox” ax” ax ' ax ax¥ ay” *” ne 


We shall have occasion to consider only symmetric affine connections: Le =e 


If g,, is any symmetric non-singular absolute tensor field, the Christoffel symbols 
based on g,,, are defined by 


Neh == = 27 (8,844 = On8ya — OB uy) , (2.3) 


where g®’ is the inverse of g,,,, g°* g,,=6%. Christoffel symbols have the trans- 
formation law (2.2) of an affine connection. The Riemann curvature tensor based 
on an affine connection TI is defined by* 


Rigs 1) = 20, Nh, + 2th (2.4) 


If Rj;;"(T) =0, the affine connection I is said to be flat or integrable. We denote 
the Riemann curvature tensor based on the Christoffel symbols {7,}, by Ri5;“(g). 
If Rj;;"(g) =0, and if g,, is positive definite, the field g,, is called a Euclidean 
metric tensor. 
The covariant derivative of a tensor field based on an affine connection I is 
the tensor field defined by 
rel, Oe nee Once a Oo 
A 


As. 

Le Ou D, : (2.5) 
If the components of the affine connection are Christoffel symbols based on a 
tensor g,,, theri we write V,, for the operator of covariant differentiation. 


An invariant property of a field ® ina spaces Y with group & is a property 
possessed in common by each of its representations @g(x"). The invariants of 


* Square brackets enclosing a set of indices denote the alternating sum over all 
permutations of the enclosed indices divided by #!, where k is the number of such 
indices. Round brackets denote the sum over all permutations divided by k!. Thus, 
for example, 4f;,;=2(4;;— i), and ag; =2 (4; + ji). 
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a field consist in all of its invariant properties and in other fields which can 
be defined in terms of it. The joint invariants of a set of fields ®,, ®,, ... consist 
in all of the invariant properties of the fields held singly and jointly. Dzfferentval 
invariants or joint differential invariants of a set of fields are joint invariant 
algebraic relations between the components of the fields and their partial deri- 
vatives of all orders with respect to the coordinates. 

The geometry of a set of fields in a space with group @ is the theory of the 
joint invariants of the fields under the group ©. 

Consider now the three 4-dimensional geometric spaces Ye, %, and SF, 
defined by the following groups of allowable coordinate transformations on four 
real variables 2”, uw =1,2,3,4*. 


I. Euclidean space-time SF, and the group Sy of Euclidean transformations : 


2g = AV (24) ef + dé (24), 
OLE SAG Pe 


4’ — 74 1 constant, 


Z 
where A” and d” are analytic functions of z4 and A¥ is an orthogonal matrix. 


II. Galilean space-time SF, and the group We of Galilean transformations: 


te At zi + u* 24+ constant, 
(2.7) 


2* = 24+ constant, 


where A* is a constant orthogonal matrix and the wu‘ are constants. 


III. Lorentzian space-time SF, and the group &y of Lorentz transformations: 
gh LH 2 (2.8) 


where L” is a constant matrix satisfying the equations 


ie Le be os) ries bh (2.9) 
{02 FOr 20 
iy ry’ O 1 O 0) 
py pry ee 
Ninemeil 00 1 0 
FO) WSs 


A motion of a material medium M in Euclidean or Galilean space-time is defined 
K 
by any three absolute scalar fields Z (z), K =1, 2,3 such that the matrix Ca (2) 
defined by pel pl 
is positive definite. 


* From this point on, Greek lower case indices will always range over the four 
values 1, 2, 3, 4. Greek upper case indices will be reserved for a variable range depend- 
ing on the context. Latin lower and upper case indices will always range over the 
three values 1, 2, 3. The summation convention applies to all types of indices. 
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A motion of a material medium M in Lorentzian space-time is defined by any 


K KL 
three absolute scalar fields Z(z”) such that the matrix X1(z“) defined by 
KL Seat, 
ar 710, Z (2.11) 


KL 
is positive definite, t.e., X+41VxV,>0 for all V~+0*. 
The geometry of a motion in “%, %, and HS will be called Euclidean 
kinematics, Galilean kinematics, and Lorentzian kinematics, respectively. 


3. Klein’s principle and general coordinates in Euclidean, 
Galilean, and Lorentzian space-time 


The use of curvilinear coordinates for 3-dimensional Euclidean space is 
familiar in mechanics. Many authors in continuum mechanics, especially in 
finite elasticity theory [11], use curvilinear and deforming spatial coordinates x‘ 
in Galilean space-time. This type of coordinate system is related to an inertial 
rectangular Cartesian coordinate system 2 by a general analytic transformation 
of the form Peay ier . 

x* = 244+ constant. 6-1) 
Unless the transformation (3.1) is a Galilean transformation, the spatial coordinate 
system x‘ is said to be non-inertial or curvilinear, or both non-inertial and curvi- 
linear. Non-inertial spatial coordinate systems in classical mechanics are further 
classified as rigid, deforming, accelerated, rotating, etc. Though the use of non- 
inertial curvilinear spatial coordinates is accepted practice in classical mechanics, 
the fourth coordinate (time) is rarely transformed more generally than in a 
Galilean transformation. Thus there has arisen a large body of literature [4, 7] 
concerned with the invariants of a motion under the more general group of 
transformations (3.1). Now the Euclidean and Galilean groups are subgroups 
of the more general transformations (3.1); however, the Lorentz group is not a 
subgroup of (3.1) since the fourth coordinate in a Lorentz transformation is 
transformed more generally than in (3.1),. The utility of introducing a more 
general class of coordinates than the z“ in space-time once accepted, there seems 
little motivation for giving undue special attention to the group (3.1) in this 
work, which attempts a uniform treatment of Euclidean, Galilean, and Lorentzian 
kinematics. What we shall do is to develop a formalism for kinematics in terms 
of invariants under the group Wx of unrestricted analytic transformations on all 
four coordinates of events. General coordinates in space-time will be denoted by 
x” and a typical element of the group Gg, is written in the form 
wt HE (x), aE (x). (3.2) 

* In §7 we introduce a group of transformations ZRLH (Z*) of the material 

coordinates Z*. At the appropriate point in the discussion, it will be shown how a 


motion in any of these spaces determines an inverse relation Mak (Z*, t) between 
the space-time coordinate 2“, the material coordinates Z* and a suitable fourth 


KL KL 
variable t. When the scalar fields C-! and X+ are considered as functions of the 
ZX and t, we shall write (C)*4, (X1)K£ consistent with the fact that these quantities, 
so regarded, transform as tensors under transformations of the material coordinates. 
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The groups ®,, Gg, and Gy, are all subgroups of &,. Once such a formalism 
for the three kinematical systems has been developed, the classical problem of 
introducing moving and deforming coordinates is of course solved since (3.1) is 
a subgroup of (3.2). That is, any invariant of a motion under (3.2) is automati- 
cally an invariant under the subgroup (3.1) of these more general transformations. 
The concepts needed to construct such a formalism for kinematics are embodied 
in KxeIn’s principle [10, p. 65]: 

If in any space with group ©, the subgroup &, is introduced, consisting in all 
transformations which leave a figure (field) ®, invariant, then the geometry of a 
figure B, with respect to WG, is identical with the geometry of the set of figures 
(B,, B,) with respect to &,. 

Let us illustrate the application of KLEIN’s principle that we intend to make 
by the following familiar example. Suppose we have given a tensor field fe 
in ordinary 3-dimensional Euclidean metric space where the group @, is the 
orthogonal group, 7.é., Bs is a Cartesian tensor. Let ©, be the group of general 
analytic coordinate transformations in 3-dimensions. (, is a subgroup of @,. 
Let g;;(x*) be an absolute symmetric positive definite tensor field under ©, such 
that its Riemann curvature tensor vanishes. Then in the space with group ©, 
we know [6, §10] that there exist preferred coordinate systems z‘ such that 
&; ;(2*) =0;;- Furthermore, any such pair of coordinate systems are related to 
each other by an orthogonal transformation. Thus the group ®, can be defined 
as the subgroup of &, which leaves the canonical form 6;; of the Euclidean metric 
field g,; invariant. Let gi-:(x*) be any field in the space with group ©, having 
any law of transformation under @, such that 


G2 (2) = FG") (3.3) 


in every preferred coordinate system in which g;;=6;;. According to KLEIN’s 
principle, the theory of the invariants of the field fi}: under the group @, is 
identical with the theory of the joint invariants of the fields (gi, g,,,) under 
the group ©, of general analytic transformations or under any group containing 
@®, as a subgroup. 

Consider the group ®, of general analytic transformations (3.2) of the four 
coordinates of events. Our objective is to define three sets of fields {©}, A =E, 
G, L, having an assigned law of transformation under ®, such that (1) there 
exists a subclass of preferred coordinates 2“ in which the fields {@}, assume 
certain canonical forms and (2) the subgroups ®;, ®¢, and G,, consist in all the 
transformations of ®, which leave invariant the canonical forms of the sets of 
fields {®},, {®},, and {@},, respectively. Once we determine such a set of 
fields we invoke KiE1n’s principle and give new but equivalent definitions of 
Euclidean, Galilean, and Lorentzian kinematics. That is, these three theories 
can then be defined as the theories of the joint invariants of the combined sets 


K 
of fields Z(x"), {B},(x") under the group G,. 


Case I. Euclidean space-time. Let t, (x) +=0 be an absolute covariant vector 
field under &, such that 


Aut = 0. (3.4) 
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Let g"”(x") be a symmetric contravariant singular absolute tensor field under 
®, such that 


REO PRE, 8, 0% (3.5) 
for all v,,==0 and not parallel to ¢,. The condition (3.4) is necessary and sufficient 
for the existence of a scalar field ¢(x") such that 

ty = yt. (3.6) 


Moreover, the field ¢(*”) is uniquely determined by (3.6) to within an additive 


constant. Let ¥ (xe), 1 =1, 2, 3 be any three analytic functions of the coordinates 
such that * Ape) 
O = &"°7 0, 90,9 0,9 Ot + 0. (3.7) 


From (3.7) it follows that we can solve for the coordinates x“ as functions of 
the variables y' and T =¢(x*) bagel 
ii (Ae Te) (3.8) 
Consider the functions g‘/(y*, T) defined by 
gi (y*, T) = g'*(%) a,9° 8,97. (3.9) 


We assume that the Riemann curvature tensor based on the positive definite 
symmetric g*/(y*, T) vanishes for all values of y’ and T corresponding to the 


events in N(E), 
(F) R;;i'(g) =0. (3.10) 
These conditions are necessary and sufficient that we be able to eee functions 


vy" r= (yi, T) such that the g” defined with respect to the y” have values 
g'7 —6'T. Thus there exists a coordinate transformation 


(3741) 


such that in the coordinate system z“ the fields g"”(z”) and t,,(2”) have the 
canonical form 


MPs 
g = 


o = So © 


O 
: t, = (0, 0, 0, 1 2 
0 ) w= (0, ? A ie G: ) 
0 


ooo) = 
ye VS 


Applying the assumed tensor law of transformation to these canonical forms we 
then see that the Euclidean group of transformations (2.6) can be defined as 
the subgroup of &, which leaves these canonical forms invariant. 


Thus Euclidean kinematics is the theory of the joint invariants under the group 


(4 of the set of frelds K 
Z (x), ghr(x™), #,(%”), (3.13) 
* vet and «,,,, are the completely antisymmetric axial tensor fields of weights 
+4 and —1, respectively, whose components are +1, —1, or 0 in every coordinate 


system and ¢1?34— ¢534= +1. 
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where the condition (2.10) invariant under ®y is replaced by the condition 
KL KL eee 

iets Ve V,> Oe Vig = 0, C1= gt? aZ ine ae (3.14) 

A coordinate system z“ in Euclidean space-time such that (3.12) holds will be 

called a Euclidean frame. We shall call gi” the space metric, and we shall call t, 

the covariant space normal. 


Case II. Galilean space-time. Let I}7,(x") be an affine connection under ®,. 
We assume that TI is a flat or integrable connection. 


BijeB iP), ==0.- (3.15) 


Let g” and ¢, be tensors under G, having the same properties assigned to these 
fields as in Case I above, but which in addition satisfy the conditions 


We — ge” +1 ge” + Ij’, ep =a 


(3.16) 
2 
Paty = Oty — Tu =0, 


jointly with the connection ite That is, the covariant derivatives of g“” and 
t, based on the connection I)7, vanish identically. 

From (3.15) follows the existence of preferred coordinate systems in which 
all of the components of the connection vanish [7, § 29]. Any two such systems 
are related by a linear transformation. From (3.16) it follows that, in any of 
the coordinate systems in which the connection vanishes, the components of g“” 
and ¢, are constants. Set z4=t(x"). This will be a linear transformation leaving 
the connection zero, and ¢, will assume its canonical form 


t= (0,0;0,A). (3.17) 
From (3.5) it then follows that g“” is reduced to the form 
mate 0 
= A As 
g ae (3-48) 


where g’ is a constant symmetric positive definite matrix. Thus by a further 
linear transformation of the first three coordinates not involving z*, preserving 
the condition (3.17), and the vanishing of the connection, we can reduce g“” to 


its canonical form 
OPO 0 


NO. -@ 
OO. 4 -@ 
OOO 


oO Oo = 


= (3.19) 


jo) 


It is then an easy matter to verify that the Galilean group (2.7) is the subgroup 
of &, which leaves invariant all three canonical forms (3.17), (3.18), and TAO 
Thus Galilean kinematics is the theory of the joint invariants under ®, of the 

set of fields K 
Z (x0) eee" (0 Vim eee incre aca (3.20) 


K 
where the Z satisfy the invariant condition (3.14). The preferred coordinate 
systems z“ in Galilean space-time in which we have (3.47)1(3a40), and bE? ag 
will be called Galilean frames, and I), will be called the Galilean connection. 
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Case III. Lorentzian space-time. Let y“’(x") be a non-singular symmetric 
tensor field with signature 2 whose Riemann curvature tensor vanishes. These 
are necessary and sufficient conditions for the existence of preferred coordinates 
2" such that [6, § 27] _ 


yi = (3.21) 


oO © 


0 
0 
1 
0 


eae i 
SO 2) 


== 


The Lorentz transformations (2.8) consists in the group of all transformations 
of G, which leave the canonical form (3.24) for y”” invariant. 


Thus Lorentzian kinematics is the theory of the joint invariants under &, of 
the set of fields K 
Z(x"),  yl(x"), (3.22) 


where the condition (2.11) invariant under Gy, is replaced by the condition 


KL _ KL K L 
X11 Veg Vi > O, Ve O, PGW A OnLy (3.23) 


invariant under &,. The preferred coordinates z in Lorentzian space-time in 
which we have (3.21) will be called Lorentz frames. The tensor field y"” and its 
inverse y,, Will be called the Lorentz metric. 


Thus we have succeeded in formulating all three kinematical systems as 
theories of the joint invariants of a motion and a suitable set of fields under a 
common group of coordinate transformations G4. Quantities transforming as 
a tensor under ®, will be called world tensors. An affine connection under &, 
will be called a world connection. The Galilean connection is a world connection, 


K 
g”, t,, and y"” are world tensors, and the Z are world scalars. 


2 "I? 


4. Euclidean and Galilean kinematics 


The results we present in this section and the following one are not intended 
to represent an exhaustive, systematic study of the invariants of a motion in 
Sz, Fo, and “,. Rather, the remarks and equations in these sections are 
intended merely to illustrate the world invariant formalism and to show the ease 
with which familiar results, often obtained otherwise by cumbersome methods, 
follow easily and elegantly. 

Since the Galilean group is a subgroup of the Euclidean group and a motion 
is defined in just the same way in both spaces, it is clear that any Euclidean 
invariant of a motion is also a Galilean invariant. Comparing the lists of fields 
(3.13) and (3.20) we see that, in the world invariant formalism, this means that 
any joint invariant of the fields (3.13) is also a joint invariant of the fields (3.20). 
Thus it is appropriate that we discuss these two kinematical theories simul- 
taneously. Any invariant of a motion in %, or %, which does not depend on 
the Galilean connection is a Euclidean invariant of the motion. However, if an 
invariant depends explicitly on the Galilean connection, then this invariant will 
be a Galilean invariant of the motion but not a Euclidean invariant of the motion. 
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Euclidean kinematics as we have defined it here is equivalent to an invariant 
theory proposed by DEFRISE [20]. We have announced the problem somewhat 
differently, but the two theories are in fact the same. The paper by DEFRISE 
contains a number of geometrical results pertaining to a motion in Euclidean 
space-time. Some of these will be included in the discussion here. The references 
contain other sources of related material. 


Consider the axial world scalar of weight 1 defined by 
1 i ape HS 
= 71 eh"0? OZ 0,L 0,2 Ob exiy (4.1) 
and the axial world vector of weight 1 defined by 
1 Bey ey. Os 
Way cee 0,2 0,2 0,2 Exrm- (4.2) 


K KL 
In a Euclidean frame, D=det 0,2 —+ Vaet CL. Since the law of trans- 
formation for D is D’=(x’ x) 4D and (x’ x) is never zero, D=-0 in any coordinate 
system. Thus we can define the absolute world contravariant vector field 


peat 


a ey (4.3) 


called the world velocity vector of the motion. The form which any world tensor 
or other type of world invariant takes in every Euclidean, Galilean, or Lorentz 
frame, depending on the context, will be called its canonical form. The canonical 
form of the world velocity vector v“ is 


ye == (v" A). (4.4) 


Since D is never zero, we can always solve for any system of general coordinates 


K 
x" in terms of the material coordinates Z*=Z and the time T =t(x"). Thus 
we always have relations of the form 


x" =x (Z¥, T), (4.5) 


where the functions % are single-valued and analytic. In terms of the x the 
world velocity vector v" is given by 


; K Lu - 
vl =37 =2(Z,%), MCAS Z) aps (4.6) 


The result (4.5) serves to promote the geometric interpretation of a motion in 
terms of a congruence of lines in space-time which are nowhere tangent to the 
surfaces ¢(x“) = T=constant. Such a surface is called an instantaneous space. 
The material coordinates Z* serve as names for the lines of the congruence and 
TY is an admissible parameter whose value is never stationary as one moves 
along a line of the congruence. A line of the congruence (4.5) is called the world 
line of the corresponding material point Z*. Each surface t(x") = T is an ordinary 
3-dimensional Euclidean metric space imbedded in 4-dimensional Euclidean or 


Galilean space-time. One can introduce a general system of parameters or 
instantaneous space coordinates rae 
ys YK) (4.7) 
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on each of the one parameter family of surfaces ¢(x”) =T. We can also arrange 
matters so that the x“ are given in terms of the y’ and T by functions 


xt = % (y*, T) (4.8) 


analytic in all four variables y’ and 7. The material coordinates Z* of a material 
medium constitute one such set of instantaneous space coordinates. The induced 
surface metric g*/(y*, T) defined by 


gil (y*, T) = gt" (#) 0,97 8,9 (4.9) 
is always Euclidean. Let the functions (C)X4(Z™, T) be defined by (cf. the 
remarks in the footnote, page 186) 

KL 
(CE ZT 1G (2))S (4.10) 


KL 

where C~1 is the set of scalar functions defined in (3.14) and the ¥ are the functions 
occurring in (4.5). If we choose material coordinates Z* for the instantaneous 
space coordinates y’, the components of the instantaneous space metric g**(Z™, T) 


are obviously given by ekL(ZM, T) = (C-)KE (2M, T). (4.14) 
Thus for any motion in Euclidean space-time or Galilean space-time we have 
Reims’ (Cpo) =0. (4.12) 


The quantities (C1)** or the inverse Cx;, Cgz(C)*™=6¥#, are called material 
measures of deformation [8, p. 140]. 

A motion in Euclidean or Galilean space-time is called vigid if and only if 
the functions (C-1)*” are independent of the time T, i.e., 
a (Co eer He — 
ae Base 0,C* =0. (4.13) 
Let gx,(Z™, T) denote the inverse of g*”. The distance between two neighboring 
material points Z* and Z*+dZ* at time T is given by 


(Gr 5e2 


dS? = gx, dZ* aZ* =CxaZ* aZ". (4.14) 


Thus a motion is rigid if and only if dS =0 for every pair of neighboring material 
points. 

The Lie derivative [10, p.106] of a tensor field o-: with respect to an 
absolute contravariant vector field v“ is a tensor field of the same type as gf: 
defined by 


£ ph =v" 0,082 — Gh et — tpt de +--+ wav phn. (4.15) 
Consider then the contravariant absolute symmetric world tensor A“” defined by 
AY’ = (ea 4£ gh” 


4.16 
= — 1 (v7 0,4” — 27” 0,0" — gt* a,v"), ( ) 
where v" is the world velocity vector of the motion. The tensor A”” has the 
canonical form Zi 0 


Aer = 
oO O 


| » @=4(0,0i+ 0,v'). (4.17) 
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The quantities d'’ are the familiar Cartesian components of the rate of deformation 
tensor [8, p. 150]. Since A”” has the canonical form (4.17), we shall call it the 
world rate of deformation tensor. Since A” is defined independently of the Galilean 
connection, it is a Euclidean invariant as well as a Galilean invariant of a motion*. 


The world scalar invariant equation 
L 


KD K 

v 0,0 = — 2A 0 Z.G 7. (4.18) 

can be easily verified by referring all quantities to a Euclidean frame. Thus a 

sufficient condition for a motion to be rigid is that 4””=0. This condition is 
also necessary, for on referring all quantities to a Euclidean frame we get 


KL ES K L K 
vd,C1=0=—2d0;2 0,2, detdO,Z+0, (4.19) 


from which it follows that d'‘7=0. Thus every component of A“” vanishes in a 
Euclidean frame if the motion is rigid. Since A”” transforms as a tensor under 
general transformations of the coordinates, its components will vanish in every 
system of coordinates if the motion is rigid. Thus the vanishing of the world rate 
of deformation tensor A” is a necessary and sufficient condition that a motion in 
SF; or So be rigid. 

The field g defined by 


B= Suede torr ke eB oe (4.20) 


is a world scalar of weight — 2 having the constant value g =1 in every Euclidean 
frame. The familiar absolute scalar invariants of the rate of deformation tensor 
d'! are given in world invariant form by the formulae 


al 


0 = SF tpoas frmtg Beg AY 0" vF, (4.21) 

0; = E - Eneas Siate gee arr oe (4.22) 
-1 

On =a Sage Bopp eA SAP a oe. (4.23) 


The canonical form of these absolute scalars is 
O, = traced’’, @, = sum of the principal minors of d‘/, @, = det d'7. (4.24) 


As an illustration of how the theory may be applied to problems in mechanics, 
consider the world tensor t#” defined by 


THY — 10, gH” + 2u AeY — peg, (4.25) 
The canonical form of this tensor is 
io 


of b= Ado + au aii — p oi, (4.26) 


THY = | 


* In pure geometry (see, e.g., [10, p. 346]), a motion in a space characterized by 
a set of fields yf: is a vector field w” satisfying one or more equations of the type 
£9. = 0. Thus we are using the word motion in quite a different sense here. 


However, much of the theory of motions in the sense of pure geometry can be applied 
to the study of motions of continuous media. 
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The quantities ¢’’ transform according to required law (1.11) under the group Gp 
relating the Euclidean frames. Equations (4.26), are the familiar constitutive equa- 
tions for the stress tensor of a classical linear viscous fluid, where A and yu are the 
viscosities and p is the pressure [8, p. 126]. The pressure is assumed to be some 


function of the world ‘scalars C-! and the temperature. From the point of view of 
world invariant kinematics, the stress t“” is a world tensor differential invariant of 
a motion in Euclidean space-time satisfying the invariant equations t#”= 1”, t#4,= 0. 

The Lie derivative of the world rate of deformation tensor A”” is again a 
tensor given by 


* 
At? == £ AMY =o 0, AM? — A? 0" + AP? 3,0" (4.27) 


The canonical form of A is 


ee. 
ste OA x. ad’! Bes ; i j i 
WA Oras A a re isla 0,41 v' — d?? 6,u' — dl? d,v". (4.28) 
*.. 
The quantities d’’ transform as a 3-dimensional tensor under &,. We see that 
the above process may be continued indefinitely to obtain an infinite sequence of 
* 


(”) 
world tensors A”’, A“”,..., A4“”,.... All of these fields have a canonical form 


similar to (4.28) involving a sequence of 3-dimensional tensor fields d’’, reel 
(n). 
d’!,.... Such a sequence of differential invariants of a motion under the group 


(, has been considered by ERICKSEN & RIVLIN [/2] in the formulation of consti- 
tutive equations for the stress in a visco-elastic material. If t”” is amy tensor 
invariant of a motion in Euclidean space-time such that it has the canonical 
form (4.26),, then ¢’? is an admissible stress tensor defined in terms of the motion 
and satisfying the transformation law (1.11) under the group Gy. Thus we see 
how the present formalism can be put to use in the problem of formulating 
admissible consititutive relations for the stress tensor in classical continuum 
mechanics. The Lie derivative of t“” with respect to the world velocity v“ has 
the canonical form 


uy pe aie ote if ast a i jt b 
m= | a Deeg Ort Ue OE Ou, (4.29) 


The quantities /*i transform as a tensor under (, and have been used in formu- 
lating the constitutive equations of a class of materials called hypo-elastic [16,27]. 


The tensor f“/ is called the convected time flux of the stress tensor. The proof of 
its invariance under ®, has been discussed from numerous points of view 


[7, 8, 13, 14, 21]. 
In terms of the world tensors g“” and v" we can define the non-singular sym- 
metric contravariant world tensor p”” given by 
pi? = gh? — yh y?, (4.30) 
Let #,, denote the inverse of ””. The canonical form of #,, is 
0;; 


—y 4+ v%0* 


buy = ed det p,,=—1. (4.31) 
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Let Sy, denote the Christoffel symbols based on the tensor Pav: Except for 
special motions, these Christoffel symbols define a world connection in Galilean 
space-time which is distinct from the Galilean connection. The Riemannian 
curvature tensor based on the {ote does not vanish except for special motions. 
The Galilean connection is independent of a motion; the Christoffel symbols 
based on #,, are a type of differential invariant of a motion in Y~; or %,. Before 
proceeding, we list some invariant algebraic relations satisfied by #,,: 


per t,.=— ve, Puvt = —t,, bio Prikl =—Pao t ye (4.32) 
The canonical form of the Christoffel symbols is 
out 


Te ae A aHE \es tg @ AS ISG. 7a eee es $s Jst Sy.5% 4 gS as¥ JST 
{jb =v d! ; ie}, ——@ 1— yr ys as? {aah er a= ose a, 


{ait a ae vu d's, {ashs = vu vs Hes. {i= hig 


where w'? = }(6,v’— 6,v'). 


(4.33) 


Since adding any world tensor S?,—S°, to the Christoffel symbols yields 
another set of quantities transforming as the components of a symmetric world 
connection [1, p. 6], we see that one can construct an infinity of world connec- 
tions all of which are part of the Euclidean and Galilean geometry of a motion. 
All that is needed is to be able to define a tensor of rank 3, S®%,, in terms of the 
motion, g“” and ¢,. The tensor #,,,v° is one such admissible tensor, there being 
an infinite number of others. Without further motivation derived from physical 
applications or intuition, there seems little to recommend an intensive study of 
this variety of differential invariants of a motion. DEFRISE [20] has based the 
determination of a world connection in Euclidean space-time upon the intuitive 
notion that the world lines of the particles of the medium shall be a system 
of “parallel straight lines” in the 4-dimensional sense, plus other ideas based 
on the parallel transport of tensor fields. Having determined such a connection, 
DEFRISE proceeds to study in detail various other invariants of a motion that 
can be defined with his connection, such as its Riemann tensor. The equations 
used by DEFRISE to determine a world connection °, (hereafter referred to as 
the Defrise connection) are equivalent to the following: 


Y vf = 0, YM [ees 2 AK” th iG ty pk (4.34) 


where the operator v, denotes covariant differentiation based on the Defrise 


connection $2,,. These equations have a unique solution for all 64 components 
of 2%, in terms of the fields g"”, t,, v“ and their derivatives. The canonical form 
of the Defrise connection is 


tiles 4 F . Ovi é 3 
OQ, ee 0, oy =0, aj ae ee 0,0 ’ 44 aime sa oe v! d,v". (4.35) 


The difference between the Defrise connection and the Christoffel symbols based 
on p,, 18 a world tensor given by 


Se eee ri Lo, a Pua pyr A** ve a Pua At, fe Poa A*et,.. (4.36) 
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The canonical form of S8, is 


EN IES aS 2 RE er RL 
Sig =v! vs vo" ds" — 205 d*®, Sie UCU ae. 
As a variant of DEFRISE’s procedure, we can solve the equations 
aR o) i = [4 9g] 4 —_ 
Ng Os snort Y, v=, Yat 0, YU, ult g 0 (4.38) 


for the components of yet another symmetric world connection Y,. According 
to (4.38),, the world lines of the material points undergoing the motion will be 
paths of the connection W% [1, § 22]. The canonical form of this connection is 


Bi=0, Yi=0, Bix, Biy=(o—ayv— (439) 


Oger’ 


and the difference between this connection and the Defrise connection is the 
world tensor Uf, given by 


UE a= PE Oe = 7 At, + PA, (4.40) 


The formulae (4.36) and (4.40) are convenient to have when one considers a 
question of the type: What is the invariant significance of the absolute derivative 
of the velocity field with respect to the Christoffel symbols Seater The answer 
follows simply from the result (4.36). Let V,, denote covariant differentiation 
based on the {%,},. We then have e 


V vt 0,08 + {Eh ve = d,0" + OF ve + Sh. (4.41) 


The first two terms on the right hand side of (4.41) cancel each other by (4.34),. 
Substituting from (4.36) for the Sf, yields 


Vet =$,,m" ee (4.42) 
If we multiply this equation by #”? and sum, we get 
pe prey, — Ae (4.43) 
Since 4“” satisfies the invariant equation 4“”4,=0, we have from (4.42) 
68 wi =v" p,, AM* = — t, AX*==0. (4.44) 


Thus the world lines of the motion are paths of the connection teehee 
Let V,, denote covariant differentiation based on the Galilean connection N Bde 
zB 


The world acceleration vector of a motion in Galilean space-time is defined by 
— vl, ve", (4.45) 
The canonical form of the world acceleration vector is 
eee rhe 
C—18t, O) ns) == a ae ce (4.46) 
The transformation law relating the components a‘ and a’ in two Galilean 


frames is eA é ai, (4.47) 
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where At is a constant orthogonal matrix. Note that the world acceleration 
vector of a motion in Euclidean space-time is undefined. 
Consider the world tensor W“” defined by 


we= git Ag v’. (4.48) 
The canonical form of this tensor 1s 
wt Me 
Ve him LC OS ow. (4.49) 
ORO 
The symmetric and antisymmetric parts of W”” are world tensors 
Av = we), Quy — wer) wer = Av’ ee (4.50) 
where the canonical form of Q“” is 
a7 
Que? ON okt 4 Oh Oy (4.51) 
@) © 


Since w‘/ is the classical measure of vorticity, we call Q"” the world vorticity 
tensor. Now the world rate of deformation tensor A“”’ was defined in (4.16) 
independently of the Galilean connection; whereas, from (4.50) it is not im- 
mediately apparent that the symmetric part of W“’ can be expressed solely in 
terms of the fields v“, g“” and the derivatives of these fields. That this is the 
case follows from (3.16). 


we) — glu a,v” + gi lu hy ve, (4.52) 
and from (3.16); we get 

tT ene (4.53) 
Substituting this last result into (4.52) we verify the identity W“”) = — 2 £ en. 


The components of the world vorticity tensor Q”” cannot in this same way be 
expressed solely in terms of the fields v”, g””, ¢, and their derivatives. That is, 
vorticity is not a Euclidean invariant of a motion. But this is clear from an 
intuitive point of view since vorticity measures a rate of rotation and rotation 
does not have an invariant significance under the Euclidean group, but does 
have an invariant significance under the smaller Galilean group. 


5. Lorentzian kinematics 

The study of invariants of a motion in Lorentzian space-time has important 
applications in relativistic mechanics and electromagnetic theory. The sym- 
metrical treatment of all four coordinates for events has found greater usage 
and favor in relativity theory than in classical mechanics. However, as we have 
seen, the formulation of classical invariant-theoretic problems in terms of a 
4-dimensional geometry is easily accomplished and has some formal manipulative 
advantages. There appears still to exist some misunderstanding in applied work 
concerning the introduction of general coordinates in space-time and the physical 
implications of such a process. The equations of classical mechanics, as well as 
all of the purely kinematical considerations we have given here, can be phrased 
and presented in terms of world invariants of a suitable set of fields in space-time. 
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This does not alter the physical hypotheses or interpretations constituting the 
theory any more or less than the familiar use of curvilinear coordinates for posi- 
tions in space. McVi1rTik [5] has considered certain of the Galilean invariants 
of a motion as limiting cases of Lorentz invariants with infinitesimal velocities. 
Here we have preferred an independent development of each of these theories 
of motion treating each as an exact science. One of the best sources for ideas 
and results in Lorentzian kinematics is M@LLER’s book on relativity theory [9]. 
Following is a brief discussion of a few Lorentz invariants of a motion whose 
importance, relative to others we might consider, is suggested by applications 
in relativistic theories of elastic bodies and fluids. 


Consider the world axial vector field of weight 1 defined by 
D 1 vau e y es 
py! ee b 0,2 Onn 0,2 EKLM: (5.1) 
The world scalar of weight 2 given by 
C= 0. (5.2) 


is always negative. This follows on substituting the definition of v” into (5.2) 
and deriving the identity RL 
6@ =dety,,, det X4 <0. (5.3) 
The inequality holds since the determinant of the Lorentz metric is always 
KL 

negative and the determinant of X! is by assumption always positive. The 
absolute world vector defined by 

pb 


wi, ty hw = —1, (5.4) 


is called the relativistic world velocity vector of the motion. 
Let the motion be referred to an arbitrary Lorentz frame. We then have 


KL Oke ED Ve ee 
X41 = 69,2 0,2 — 0,2 a,Z, 


wwe Ke LE I agrmas (5.5) 
07 0;2 02=X4- 0,2 OZ. 
Taking the determinant of both sides of this last equation yields 
K KL K £ 
D* = (det 2,2)? = det X* (1+ XZ %Z) >0, (5.6) 


KL LEQe 
where X is the inverse of the matrix X~. Since X~ is positive definite, so is its 
KL 


inverse, and from (5.6) we see that D* is never zero for any motion of a material 
medium in Lorentzian space-time. This means that we can always solve for the 
first three Lorentz coordinates z’ of a Lorentz frame in terms of the Z* and 24 


=i (ZX, aA), (5.7) 
The canonical form of the axial vector ” is 


as 
Oz4’ 


b= D/ 1) =D(o',1), (5.8) 


14% 
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where v' is the “classical velocity” of a material point relative to the Lorentz 
frame 2", where we think of 24 as classical time. This is not a Lorentz invariant 
notion since z4 and T have different transformation laws. The canonical form 


Saliais @ = D2 (v' v'— 1) <0, (5.9) 
where the inequality follows from (5.3). Thus we see that the hypothesis that 


x1 is positive definite for the motion of a material medium in Lorentzian space- 
time leads to the familiar result that the “classical velocity’ of the motion 
relative to any Lorentz frame is always less than 1, where this upper limit for 
the speed of any motion is identified with the speed of light in the chosen system 
of units. 

The components of the relativistic world velocity vector w“ in a Lorentz 
frame have the values 

ie ut 1 2 ae ae 
w aoe eal vs vv". (5.10) 

In relativistic mechanics, the fourth component of w“ in a Lorentz frame less 1 
is called the kinetic energy per unit of mass [4]. 

An affine connection in “, is determined by the Christoffel symbols based 
on the Lorentz metric. All of these symbols vanish in a Lorentz frame. Let 
Yi denote covariant differentiation based on the Christoffel symbols of y,,. 


The world velocity gradients are defined by Vw”. The relativistic counterparts 
of the world rate of deformation tensor and the world vorticity tensor of Euclidean 
and Galilean kinematics can be defined as follows: 


At? = A (Ee v) ir Alu ny)? | 3 
May Vw, Oia Y, w (5.11) 


where we shall have the identity 
AN = Sees (5.12) 
analogous to the classical case (4.16), (4.50). 


In the case of Euclidean or Galilean kinematics, a rigid motion can be defined 


. . . KL 
by either of the conditions v“d,C-'=0, or A“”=0. In the case of a motion in 
Lorentzian space-time, however, the two analogous conditions 
KL 


wid,X+=0, AM =0 (5.13) 


are not equivalent. There exist motions for which we have (5.13), and do not 
have (5.13)2. In fact, if w satisfies (5.13),, then it is a translation [10, p. 349]. 


Its components in a Lorentz frame will have the form (5.10), where the v‘ will 
be constants. The motion 


Tey : : ; 
Z (2, 24) = AK (cA) zi, AK AL = 6X4, (5.14) 


representing a rotation about the origin of the Lorentz frame z‘ is an admissible 
motion in “, for all values of z and z4 such that 
KL 


At ORE AR AE 7 zh (5.15) 
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is positive definite. This will be true for all sufficiently small values of z’. If 
the time derivatives of the A* are constants, then the first of the conditions 
(5.13) is satisfied. The motion (5.14) is a rigid motion in Euclidean or Galilean 
space-time. Thus it is reasonable to define a rigid motion in Lorentzian space- 
time by the first of the conditions (5.13). We are aware that there has been con- 
siderable debate as to what a useful and appropriate definition of a rigid motion 
in relativistic kinematics should be. We see that the conditions (5.13), would 
be too restrictive since they rule out all but the uniform translations. Pursuing 
the analogy with Euclidean kinematics, consider the functions (X71)*"(Z*, 24) 
defined by rae 
(CS Conran (e(2 es"), 24), (5.16) 


where the functions 2 are those occurring in (5.7). Thus, treating z4as a parameter, 
we can construct the Riemann curvature tensor Rz yn" (Xx1), Xx, (X24 =O", 
MOLLER [9] calls a motion in .%, satisfying the conditions 


Re wn* (Xpe) =0 (5.17) 


a Euclidean motion. The rigid motion (5.14) is not a Euclidean motion. Equation 
(5.17) is to be compared with its classical counterpart (4.12), which holds for 
any motion in “; or “%,. Moreover, (4.12) is a Euclidean and Galilean invariant 
property of any motion in “, or %; whereas, if (5.17) is satisfied by a motion 
relative to one Lorentz frame, it need not be satisfied by the same motion 
referred to another Lorentz frame. That is, since z4 is not an invariant para- 
meter in (5.47) under Lorentz transformations, this condition is not Lorentz 
invariant. 

In special relativity theory, the equations representing conservation of energy, 
conservation of momentum, conservation of angular momentum, and the equi- 
valence of momentum and energy flux take the form [4] 


V, Pe = (0; Umea hey (5.18) 
y 
where P“” is the stress-energy-momentum tensor. It is customary to write P”” 
in the form 
Pie 7h” — o.u0t wo", (5.19) 


where t“” embodies the relativistic counterparts of the classical stress, internal 
energy, and heat flux. In classical mechanics, it is customary to provide con- 
stitutive equations for all of these quantities in terms of the motion, the tem- 
perature, the electromagnetic field, e¢c., and to require their invariance under 
various groups of transformations such as the rigid motions. But in relativity 
theory, owing to the open question of a proper definition of a rigid motion which 
is not a uniform translation, a clear and concise statement of the class of ad- 
missible constitutive equations for t“” depending on a motion cannot, to our 
knowledge, be found in the literature. The relativistic counterpart of the funda- 
mental assumption (1.11) of classical continuum mechanics has not been enuncrated. 
It would seem that a sound relativistic generalization of, say, classical elasticity 
theory rests on questions of this nature. 
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6. Convected coordinates in Euclidean and Galilean space-time 

Given a motion in Galilean or Euclidean space-time, there are, in addition 
to the Galilean and Euclidean frames, other classes of preferred coordinates 
defined in a natural way in terms of the motion [12, 14]. One such class of 
coordinate systems are the convected coordinates belonging to a motion (JJ, 7]. 
A convected coordinate system of a motion in -%; or “%, can be defined as any 
coordinate system x® (we use upper case Greek indices to indicate convected 
coordinates) in which the world velocity vector v“ and the covariant space nor- 
mal ¢,, have the convected form 


1-210, 0/0) “to =(0,0, 0A (6.1) 


The existence of such coordinate systems is easy to perceive. Let x (Z*) be any 
three functionally independent functions of the Z*. If we transform the com- 
ponents v“ and ¢, in an arbitrary system of coordinates x” to the coordinate 
system x® determined by ze 
ies XZ, (0), 


Meat (4h) aed, Ce 


we conclude immediately that the coordinate system (x*, T) is a convected 
coordinate system. Any two convected coordinate systems x” and x® are 
related by a transformation having the general form 


x — & (x*) , 
ee (63) 
x” =x*-+ constant. 
From (3.5), it follows immediately that in a convected frame 
eee (6.4) 
and from the definition of the scalars Ca we have 
2 K B 
C-1— pun _0Z 02 6.5) 


ax! axN - 


. . . . K K 
If we consider the special case in which X=Z, we have the simpler relation 
ene 
Ca ae (6.6) 
Thus the non-vanishing components of g"” in a convected frame x* =Z* of a motion 


are equal respectively to the six material measures of deformation Co, Some workers 
in elasticity theory [17] who employ convected coordinates almost exclusively 
refer to the quantities Cx, as the metric tensor. The result (6.6) provides the 
principal motivation for the use of such terminology. However, it must be realized 
that the equality (6.6) holds only in a restricted class of coordinate systems. 

Consider next the rate of deformation tensor A“’, whose components in a 


general frame are given by (4.16). Substituting the convected form of the velocity 
vector into this formula, we conclude that 


PHL elie) — core Ao gie 
A Ala == 05" ioe = as hae (6.7) 
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n) 
More generally, we have for the series of world tensors Ae, wat Ae, eee (4:27) 
(ny (n) (n) n oK L 
24 74a _ Kiseree AC e 
7A fied base Ue = oT: (6.8) 


Since the covariant derivatives of g“” and ¢, with respect to the Galilean 
connection vanish in an arbitrary coordinate system, we have 


Jota — Tee ty =O, (6.9) 
yg? ? + Teg e4? + Ty, e°4 = 0. (6.10) 


From (6.9) and (6.1), it follows that, in a convected frame, 


VD sails (6.11) 
From (6.10) we obtain the equations 
1 @ekKLl 
Ty = (8), Tie —— 4 28E 8 ee, (6.12) 


where the I74, are the Christoffel symbols based on gx7, gx, g°"=6™. 
The components of the world vorticity tensor 2°” in a convected frame are 


given by 204 = —Ot9=09, QKL = Mik PU, (6.13) 
The world acceleration vector has components given by 
@=0, ak =k. (6.14) 


The formulae of this section are useful for the interpretation of any system 
of deforming and accelerated coordinates in Euclidean or Galilean space-time 
not necessarily associated with the motion of a material medium. That is, consider 
the class of all coordinate systems in “,; or Y, in which we have ¢,= (0, 0, 0, 1). 
Equation (6.7) then gives an interpretation of the time dependence of the non- 


vanishing components of the space metric g“” in one of these general types of 


KL 
coordinate systems. —+ a is a measure of the rate of deformation of the 


coordinate system. Similarly, the formulae (6.12), (6.13), and (6.14) provide an 
interpretation of the non-vanishing components of the Galilean connection in 
such a system of coordinates. 


7. The geometry of the space of material points and material symmetry 

In continuum mechanics we assign a geometry to the 3-dimensional space 
Sy, with coordinates Z* by introducing a group Gy, of allowable material coordinate 
transformations. A typical element of this group has the form 


se cea Bi (7.1) 


where, in all of the applications with which we are familiar, it is sufficiently 
general to assume that ®y is some subgroup of the 3-dimensional orthogonal 
group. Thus the matrix of coefficients in (7.1) is an orthogonal matrix. By 
demanding invariance of constitutive equations under © we obtain further 
restrictions on the form of these equations. The group ©, determines the material 
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symmetry of the medium. Let us see how this relation between the group of 
material coordinate transformations ®, and the idea of material symmetry 


nee 
arises. It is customary to require that, for some instant 7), the functions Z (z’, T) 
representing a motion relative to a Euclidean or Galilean frame (2, I) reduce 


to the form Pee 
=O Z12, Ip) = Ogee (7-2) 


That is, the material coordinates Z* at the instant J) coincide with the Cartesian 
coordinates z’. Now if the material points are all identical and arranged in space 
in a uniform homogeneous array, the intuitive notion is that we cannot in this 
way ascribe unique names to each material point, but that the names (material 
coordinates Z*) are determined only to within an arbitrary orthogonal trans- 
formation. This is the intuitive picture of an isotropic homogeneous material 
medium. Thus the appropriate group ®,, for a material with this symmetry is 
the complete orthogonal group. If some of the material points have different 
properties than others, such as in a crystal, and if they are arranged in space 
in some non-uniform or inhomogeneous array at the instant 7), the class of 
equivalent (allowable) material coordinate systems Z* will be smaller than the 
corresponding class for isotropic homogeneous materials. Thus the group Gy 
in the general case will be some proper subgroup of the orthogonal group. The 
continuum theory of elastic homogeneous crystals is based on constitutive 
equations for the stress which are invariant under a group of material coordinate 
transformations ®y,, where the set of matrices A*’ in (7.1) constitute the elements 
of one of the 32 crystallographic subgroups of the orthogonal group characteristic 
of the point symmetry of the crystal. In the theory of finite elastic deformations 
of crystals and in the classical linear theory of elastic crystals the constitutive 
equations are assumed invariant to arbitrary translations D* in (7.1). This last 
assumption represents an approximation to a more detailed description of the 
symmetry of a crystal in which one would require invariance only to a group 
of discrete translations DX. However, the present formalism does not rule out 
the possibility of an accurate and detailed continuum description of such ‘‘micro- 
scopic”’ structure or symmetry. Other choices for the group &,, describe materials 
having transverse isotropic symmetry, orthotropic symmetry, eéc., [19]. 

Let us illustrate the way in which the invariance of constitutive equations 
under the group Gy restricts their general form by considering the case of the 
stress tensor in finite elasticity theory. In §4 we remarked that if the stress 
tensor ¢"’ in an elastic body of any symmetry whatever were invariant under the 
group of rigid motions, then it must reduce to the form 


PZ Oy 7.28) 


If we now demand that ¢’/ transform as an absolute scalar under the group Gy, 
of all orthogonal transformations, whereby we assume that the material is 
. . . . to 

isotropic and homogeneous, then it is known that the functions P are expressible 
in the special form [72] 


IQUE, 
Pa FRU + FCH4 FC, 73) 
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. eZ . . . KL 
where the ¥’s are functions of the scalar invariants of the matrix C7. 
A7o= 7p (Leal; U1); O==0,41,2; 


KL KL KL 
I=trace C~*, Il =sum of the principal minors of C-1, III = det C-!'=92. 


(7.4) 


The constitutive equations for the stress in a material with less symmetry than 


. . . . . . KL K KL 
an isotropic homogeneous medium will involve functions P of Z, C71, and D 
more complicated than (7.3). 


In presenting these few remarks and examples concerning the relation between 
material symmetry and the invariant-theoretic problems encountered in the 
formulation of constitutive relations in continuum mechanics, we have attempted 
to make clear that the problem is one of invariance under at least two distinct 
groups of transformations: 1) invariance to rigid motions, 2) invariance under 
a group of material coordinate transformations. A third group not considered 
here is the group of wnit transformations. It is important that these three 
demands for invariance not be confused and that each of them be satisfied. A 
confusion of this sort is the apparent source of difficulty in some recent attacks 
on the foundations of classical elasticity theory [18, et al.}|. As we have seen, 
these invariances are not equivalent. 

In the 3-dimensional space Y, with group &\y we can introduce general 
coordinates X* and appeal once again to KLEIN’s principle so as to obtain an 
equivalent statement of the invariant theoretic problem. Let &,, denote the 
group of unrestricted analytic transformations on the material coordinates X* 
with typical element * 2 Satie? Rial ets: 

ee Y. ( Xe eee A Ie (75) 
Let H,(X*), Q=1,2,...,N bea set of fields having an assigned law of trans- 
formation under &,,4 such that the fields Hg possess certain canonical forms 
Hy (Z*) in a subclass Z* of the general coordinate systems X* and such that the 
group &, consists in all the transformations of ©, which leave these canonical 
forms invariant. By KLEIN’s principle, the geometry of a field or set of fields 
such as Cx, (Z™, T) with respect to the group ®y is equivalent to the geometry 
of the fields Cx, (X™, T) together with the fields Hy (X*) under the group Gya. 


Since we assume that Gy is some subgroup of the orthogonal group, we can 
always choose for one of the fields Hg a symmetric positive definite Euchdean 
metric tensor Gx (X™) and identify the frames Z* with some subclass of the coordinate 
systems in which we have Gx,(Z“) =dx 7. If Gy is the complete orthogonal 
group, Gxz is the only field in the set Hg since the orthogonal group consists 
in all the transformations which leave its canonical form 6x; invariant. The 
procedure of introducing general coordinates in the space Ay of an isotropic 
homogeneous material can be illustrated by writing down the form of (7.3) which 
is invariant under ®y,. It is as follows [8]: 


PKL (XM, T) = F,GRE(X™) + FCA) (KT) FACIE, (7.6) 


K 
* Here we must use a different symbol for the functions Y representing a material 
coordinate transformation so as not to confuse these functions of 3 variables with 


K : : 
the motion X(#”) represented by functions of 4 variables +, 
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where we put (C2)F =Gpy(C)E™, (C2)*4=G*™ (Cn (C4). The F's are 
functions of I, II, and III, now given by 


TS» 


T=(CYE, IL =>, OBR(CHE (CM, IL = 2 OB RY (CA)R(CHE(CAE. (7-7) 


8. Two-point tensor fields and world invariant kinematics 


We have seen that in continuum mechanics we are interested in functions 
of a motion such as the stress tensor in elasticity theory which are invariant 
under two groups of transformations. The first of these groups involves trans- 
formations of the coordinates of points (events) in a 4-dimensional space and the 
second of these groups involves transformations of the coordinates of points 
(material points) in a 3-dimensional space. Multiple point fields are familiar 
objects in pure geometry and have been used to advantage in continuum me- 
chanics [2, 19, 17]. The concept of a 2-point field is an easy generalization of 
the concept of a 1-point field set forth in §2. Let “ and XY be two spaces, 
not necessarily of the same dimension, with groups ©, and ,, respectively. 
Let X°, Q=1,2,...,.N denote coordinates for the points in “ and let x, 
fu =1,2,...,m denote coordinates for points in &. A 2-point field is a set of 
functions F,(X®, x) of the coordinates of a point in “ and of a point in 4. 
Thus the components F, of a 2-point field will depend in general on n+ N 
independent variables. A 2-point field has a law of transformation for its com- 
ponents under independent transformations of the coordinates of either point. 
The law of transformation may be any law consistent with the property that a 
representation F,(X®, x“) in any system of reference (X®, x“) determines a 
unique representation F,(X®, x”) in every other allowable system of reference 
(X?, x). An invariant property of a 2-point field is a property held in common 
by all of its representations. Joint invariants, differential invariants, efc., are 
defined as in the case of 1-point fields. 

A 2-point tensor field is a set of functions T¥-§-:(X°, x“) with a trans- 
formation law of the general form 

Ty 6. = |(x! x)|-* |(X'X)|-™ (sgn (x’ x)? (sen (X' X))¥ x 
Ot OX Oe Te 
ant ax® oe” OXF 


(8.1) 


fb ses Sonzee 
wR. aks Jes 


An absolute 2-point tensor field is one for which w=W—=y=Y=0. Names of 
2-point fields with other values of the exponents w, W, y, and Y are assigned 
on a basis similar to the case of 1-point fields. 1-point tensor fields are special 
cases of 2-point tensor fields which are constant scalars with respect to trans- 
formations of one of the points. Thus, in all of the work preceding this section, 


the 1-point fields in space-time can be regarded as special cases of 2-point 
fields. 


Let 4 be the 3-dimensional space A; of material points with group ®y4 
and let Y, be one of the 4-dimensional spaces A, SL, or SF, with group G,. 
The invariant theoretic problems of classical continuum mechanics can be phrased 
in terms of 2-point fields F,(X*, x). We call such a 2-point field a world-material 
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field. We have already met with an example of such a world-material field 
(C1)X*(X™, T) that we defined by setting 


(CEE (X™, T) = CE (X™, T)). (8.2) 


One readily verifies that (C)*" is a world-material tensor field having a trans- 
formation law of the general form (8.1). The dependence of this field on the 
space-time coordinates x“ is somewhat special since it depends on the x“ only 
in the combination T =i(«“). Now it is possible also to regard (C4)** as a 


KL 
1-point field in space-time, C~*(x"). But these quantities do not have the con- 


. ‘ KL 
venient tensor law of transformation under the group Gy4. Rather, the C~!(x") 
have the odd transformation law 


KL’ Rly itera Re 
CH (a) = B (a!) B (xt) CA (x), (8.3) 
bea 
where the functions B are determined by 
my = 28 (Ep 
B(x) = 22 (x(a), (8.4) 


K 

and the X («") are the functions representing the motion relative to the general 
KL 

coordinate system x”. Thus the transformation law of the C~1(x) under the 


group ),, depends on the motion x (x“) and cannot be written down without it. 
Of course, we used the motion to define the functions (C1)**(X*, #), but they 
have a more convenient tensor law of transformation under ®y,,4. 

The concept of a 2-point field and a pair of spaces %, %, or A and Ay 
each with its own geometry leads naturally to the interpretation of a motion 


Ge 2 (x“) as a mapping between the points of space-time and the material 
points of A,. The mapping is one-to-one only in the direction x“ X*. In the 
other direction X* > x", a single material point is mapped onto a 1-dimensional 
set of events, the world line of X*. 


The geometry of world-material fields is enriched still further by adding the 
connectors or shifters gy (X™, x), gk (X™, x”) to the list of 1-point fields g*’, ¢,, 
Tye, y=, HEL (X™), ..., Gx, (X™) which characterize the geometry of each of 
our two spaces separately. The connectors are 2-point world-material absolute 
tensor fields providing a linkage or connection between the two spaces. Quantities 
similar to the connectors we now introduce have been considered previously 
in [19, 17]. For our purposes here, we shall define the components of the con- 
nectors in a general system of coordinates (X*, x) as follows: Let there be a 
class of preferred Euclidean, Galilean, or Lorentz frames z/ and rectangular 
Cartesian material frames Z* such that the spatial frames z, and the material 
frames ZX are “at rest’ relative to each other and whose axes coincide. With 
respect to such a system of reference (Z{, 24) we assume that the connectors 
have the joint canonical form 


ek = (54.0), ef = (88,0). (8.5) 
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The components of the connectors ina general system of reference are then given by 


. L LU K 4 

ge (X™, x”) = 0; a we gk (XM x) = OF eS Ss (8.6) 
Stated more simply and directly, we assume that the connectors are absolute 
2-point tensor fields such that by suitable choice of a Euclidean, Galilean, or 
Lorentz frame and a Cartesian material frame they are reducible to the joint 
canonical form (8.5). A pair of coordinate systems (Z*, z*) in which we have 
(8.5) is called a common frame. If the components of the connectors corresponding 
to A and FY, are referred to an arbitrary Cartesian material frame Z* and an 
arbitrary Euclidean or Galilean frame z“ they will have the canonical form 


Sk = (Sx, 0), be hea Vale (8.7) 
aes K yi 4 CK __ st ; 
VE = — SEV", SSF = 8, 

where V‘ is the relative velocity of the origin of the spatial frame z‘ and the 
material frame Z*, and Si. is an orthogonal matrix representing the relative 
orientation of the two systems of axes Z* and z‘. In Euclidean space-time, the 
S* and V' will be general functions of z4, while in Galilean space-time they will 
be constants. Let V“ be a world vector in Y% or “ with the canonical form 


Vites(ViA)s (8.8) 
The Euclidean, Galilean and Lorentzian connectors satisfy the invariant relations 
eee = OR, beer =o —V"i,, 
gh ey ei CR LG gee Seat (8.9) 
PRE hat ORs cela fee pee 


All of these relations can be verified by referring all quantities to the common 
frame. Since they are tensor equations holding in one frame, they will hold in 
a general system of coordinates. 

As an illustration of the kind of world invariants of a motion in A or G 
that one can construct with the connectors, we consider the problem of defining 
a world tensor measure of the finite rotation of a motion relative to the common 
frame. The considerations given here are natural generalizations of those given 
in [17, § 4| to the case of arbitrary moving and deforming coordinate systems in 
Space-time. 


Consider the world tensor defined by 


Cun = Gc (X (x7)) 8, X¥ 3, XE, (8.10) 
The canonical form of this tensor is 
eee heared —C,,v" 
wy — k na |e. COL Oe hy Aaa (8.41) 


ae OF UGS Cr Oe 


The quantities c;; are the spatial measures of finite deformation introduced into 
elasticity theory by CaAucny and GREEN [8]. It follows immediately from the 
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canonical form (8.11) that det c,,—0 so that c,, is a singular tensor. Its null 
eigenvector is the world velocity vector v“. Thus in defining a world spatial 


measure of finite deformation it proves more convenient to use the tensor c*“” 
defined by 
cr = GO x Oy x", (8.12) 


where x (XX, T) are the functions (4.5). The canonical form of c*"” is 


UY (C22 0 
c =| A a (8.13) 


where the (c7)'7 are the components of the inverse of c;;. e**” has the null 


eigenvector ¢, and (c™)'/ is positive definite. Consider the eigenvalue equation 


cern, = cA gh’ n,, NO fez boa (8.14) 


This equation has 3 solutions (m,, c) with positive c and vectors n, satisfying 
bes 22 Q 2 
the two conditions 


gh’ My, % = Soe, ve). (8.15) 


The scalars ( é — 1) are called the principal extension ratios. The canonical form 


of the vectors My is 
=(n,, —v'n,), 8.16 
1, = (i, — 0%) (8.16) 


where the unit vectors ”; determine the principal axes of strain in the deformed 
2 


body [8,17]. All of these results follow from the canonical form of c*”’, and 
we have simply placed them in world invariant form. 


Now consider the world-material vector fields N (X*, T) satisfying the eigen- 
value equation 
(C>)K* N, ar GY N,. (8.17) 


Since (Grays is positive definite, the eigenvalues oo are all positive, and there 


exist three linearly independent eigenvectors N, satisfying 
KL = 
G Nk Ni = 900. (8.18) 


Consider next the world-material vector fields Mx (X™, T) obtained from the 
fields Mn (x") according to the rule 


ty (X™, T) = ge (X™, #) (4) (8.19) 
Since we have the identity (8.9),, it follows from (8.15) and (8.19) that 


GRY tx oe = doe: (8.20) 
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Thus Ny and mx are two sets of orthogonal unit vectors at (X*® T)y therefore 
os & 7 . . 
there exists a unique matrix Rf (X™, T) satisfying 


a VON MN pK RPL _ CKL ae KN f 
ya ii Nx: G Ry Ri G 5 Aes oft (8 21) 


When the material system of coordinates is rectangular Cartesian, Rf vis an 
orthogonal matrix. The sense in which this matrix is a measure of the finite 
rotation of a motion and the relation in which it stands to the classical measure 
of infinitesimal rotation has been explained in [19, 17]. From the canonical 
forms of c*“” and (C+)¥¥ we see that c;; and (C)*” have equal eigenvalues. 
Since these eigenvalues are absolute scalars under general transformations of the 
coordinates (X*, x"), we shall have in general 


Ne fel 
e(%) =aiG ; (8.22) 


provided we order the two sets of eigenvalues appropriately. With this result 
we can show that the vector fields Mn given by 


K a= 
n, = 0,X* Ng ({X, 1) VC (8.23) 


satisfy all of the equations (8.14) and (8.15). If the eigenvalues of Cx; are 
distinct so are the eigenvalues of c*”” and the eigenvectors Nk and um, are uniquely 


determined. This is not true if two or more of the eigenvalues are equal. However, 
in the case of distinct eigenvalues, the motion determines a unique matrix R*¥ 
provided we order each set of eigenvectors in some definite way such as that 
corresponding to the equalities (8.22), (8.23). 


Multiplying (8.23) by Sir Ni and summing on uw and Q we get 


Gur RE = Dim, Nghe =D tu Ne = 0,X¥ the SVC NeM, 


Eyl ath 8.24 

Ryt= 0, x gir (C*) x. 

Multiplying this last equation through by (C~3)4% ge” and using (8.9), we get 
finally 

O20 ee 0,X* + (C~4)¥" Ru e™. (8.25) 


The canonical form of this world-material tensor equation is 


O;X® = (CEY Ry Si, (8.26) 


axk “4 = 
GE AICHE Rigs VEL ViO.XE 4 Oe (8.27) 


Equation (8.26) corresponds to the result (4.19) of [17]. The last equation is 
an identity satisfied as a consequence of (8.26) and (8.7)s. Equation (8.26) is 
the familiar decomposition of a deformation 0,;X* into a pure stretching without 
rotation followed by a rigid rotation. In Euclidean space-time, the orthogonal 
matrix SM will in general depend on the time. As we have said, it represents 
the time dependent relation between the material coordinate axes and the 
moving, rotating axes of the Euclidean frame. 
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Die Stabilitat der Strimung in einem gekriimmten Kanal 


GUNTHER HAMMERLIN 
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1. Einleitung 


Seit G.I. TAyLor im Jahre 1923 seine Untersuchungen iiber die Stabilitat 
der Strémung einer zahen Fliissigkeit verdffentlichte [1], die sich zwischen zwei 
rotierenden achsengleichen Zylindern befindet, gewann die Betrachtung der dort 
zum ersten Male aufgetretenen stehenden Wirbel mehr und mehr an Interesse. 

Das Wesentliche an dem Taylorschen Modell ist es, daB eine Str6mung langs 
einer gekriimmten Wand vorhanden ist, die sich in einem labilen Zustand be- 
findet. W. R. DEAN [2] tibertrug im Jahr 1928 den Taylorschen Grundgedanken 
der Stérung einer laminaren Strémung durch Uberlagerung einer stehenden 
wirbelartigen Instabilitat auf die Strémung in einem gekriimmten, unendlich 
hohen Kanal; die Wande des Kanals werden dabei von zwei achsengleichen 
Zylindern gebildet. Kiirzlich wurde dieselbe Aufgabenstellung von C. S. YrH und 
W. M. SancsTER [3] aufgegriffen, wobei diese Autoren die Arbeit [2] offenbar 
nicht kannten. 

Der Ansatz einer Taylorschen Instabilitat in die vollen Navier-Stokesschen 
Gleichungen der rotationssymmetrischen Kanalstrémung fiihrt auf ein Eigen- 
wertproblem der Art, wie es schon bei TAYLor auftritt. Als Parameter treten 
darin die Wirbeldicke der angesetzten Instabilitat und deren Anfachungskonstante 
auf, ee wir als den davon abhangigen Eigenwertparameter eine GréBe 


Re? z - betrachten wollen (d= Kanalbreite, R, = Kriimmungsradius der inneren 


Kanalberandung und Re = Reynoldssche Zahl der Grundstrémung, gebildet 
mit d). 

Unter den Eigenwerten interessieren dabei vor allem diejenigen neutraler 
Wirbel und unter ihnen wiederum derjenige kleinsten Betrages, den man gewoéhn- 
lich als ,,kritischen Wert“ bezeichnet. Seine Bedeutung liegt darin, daB eine 
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a 


Strémung sicher dann nicht gegen Instabilitaten der betrachteten Art anfallig 
ist, wenn der ihr zukommende Wert des Parameters Re2 unter dem kriti- 


schen liegt. Zu diesem kritischen Wert gehért auch eine bestimmte Wirbeldicke, 
die also den gefahrlichsten der neutralen Wirbel eigentiimlich ist. Das Aussehen 
eines solchen Wirbels kann man beschreiben, wenn man die zu dem kleinsten 
Eigenwert gehérenden Eigenfunktionen kennt, die gerade die Geschwindigkeits- 
komponenten der Wirbel darstellen. 


In der Arbeit von W. R. DEAN wird der kritische Wert angegeben, und es 
wird gezeigt, daB auch gedampfte und angefachte Stdrungen existieren. Uber die 
Eigenfunktionen und damit das Aussehen der Wirbel werden keine Aussagen 
gemacht. YIH und SANGSTER behandeln nur neutrale Stérungen und geben 
einen kritischen Wert an, der zu dem Deanschen im Verhiltnis 1:150 steht. 


Die vorliegende Arbeit verfolgt zwei Ziele: Einmal wird eine Neubehandlung 
des Eigenwertproblems fiir neutrale Stérungen mitgeteilt, die es erméglicht, die 
Eigenwerte kleinsten Betrages in Abhangigkeit von der Dicke der angesetzten 
Wirbel genau zu bestimmen. Es wird sich zeigen, daB die so gefundenen Werte 
mit den von DEAN angegebenen iibereinstimmen. Des weiteren werden wir auch 
die Eigenfunktionen finden, die zu diesen Eigenwerten gehdren, so daB damit 
die Gestalt der auftretenden Wirbel beschrieben werden kann. Die Genauigkeit 
unseres Vorgehens erméglicht es, die Resultate numerisch zu sichern. 


2. Die Stérungsdifferentialgleichungen 


a) Zur Untersuchung der Strémung in einem gekriimmten Kanal hat man 
von den vollen Navier-Stokesschen Bewegungsgleichungen dieser Str6mung aus- 
zugehen. Wir wahlen die tibliche Darstellung in 
Zylinderkoordinaten ; dabei seien 7 die radiale Rich- 
tung, g das WinkelmaB und z die Koordinate in 
Richtung der Zylinderachse. R, sei der Radius des 
inneren, R, der des 4uBeren begrenzenden Zylinders 
(Fig. 1). Die Geschwindigkeitskomponenten werden 
folgendermaBen benannt: 


g-Richtung: u 
y-Richtung: v 
z-Richtung: w. 


Man kennt eine exakte, stationdre Lésung der Fig. 1. Das Koordinatensystem 
infolge der Rotationssymmetrie des Problems ver- 
einfachten Bewegungsgleichungen, die die Randbedingungen u(R,)=v(R)) = 
w(R,) =0, “(R,) =v(R,) =w(R,) =0 befriedigt: 


R, Ofo { RzIn(R,/R;) ( y =) net Int , 
2 


lige cM Rite Rial. Rak ba tole 
V,=0, 
Wy = 0 (GoLDSsTEIN [4], S. 315). 
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Fiir den Druck f, ergibt sich dabei 
Po=ko +t), 


also 

oPo =k,  konstant. 

op 

b) Um die Stabilitat der Grundstrémung wg, v9, W) zu untersuchen, tiber- 

lagern wir dieser kleine Stérungen u*, v*, w*; d.h. wir gehen mit dem Ansatz 
U=Up(r) + u(r, z, t); v=v*(r, z, 1); w=w*(r, 2, 1); P=Pol”, v) + P*(”, z, 2) in die 
vollen Navier-Stokesschen Differentialgleichungen. Diese werden unter Beriick- 
sichtigung der Tatsache, daB wu) und #9 stationdare Lésungen sind, hinsichtlich 
der gesternten Stérgeschwindigkeiten linearisiert. Im Hinblick auf die geometri- 
schen Verhdltnisse unserer Aufgabe machen wir dann Gebrauch von der An- 
nahme, daB die Breite des durchstrémten Kanals klein sei gegeniiber den Radien 
seiner Wande, also R,—R,=d<R,. Da die Stérungen g* naturgemaB auf das 


Innere des Kanals beschrankt sind, folgt in bekannter Weise 4 gt<— ge, und 


Br K Bry. 
Mit diesen Annahmen, die eine erste Approximation hinsichtlich der geo- 


metrischen Verhaltnisse darstellen, reduzieren sich die Stérungsdifferentialglei- 
chungen schlieBlich auf 


(4.1) u;* 4 ve ( Se Mo) = (eh + Uss) » 

(1.2) Pty =p toot tot), 
(1.3) wit == Gb: trent wh), 
(1.4) Uae 1 —4iy 


Als Randbedingungen hat man das Verschwinden der Stérgeschwindigkeiten 
an den Berandungen zu fordern, also 


u*(Ry) = v*(Ry) = w*(Ry) = 0 
u*(R,) = v*(R,) = w*(R,) = 0. 
Fiir das Folgende fiihren wir anstelle von 7 eine neue Koordinate nein: 


free aE ok r= R,>y=0 
R,— Rk r7=Ri >= 1. 


Fiir die gesternten Stérglieder machen wir nun nach dem Vorbild von TAYLOR 
den speziellen Ansatz 


u* = u,(n) cos az eF* p* = p,(n) cos az e** 
v* = v,(y) cos a ze = 270A 


w* = w,(y)sinaze®* 4/2 = Wirbeldicke. 
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Unter der schon oben benutzten Annahme d< R, wird die Grundstrémung 
u(r) in demselben Grad der Approximation, in dem die Gleichungen (1) gelten, zu 


Mon) = KS (0), 
Ry 
1 Rk, OP 1 
wobei der konstante Wert ii ie <0 gesetzt ist. 


Um die Stérungsgleichungen dimensionslos zu machen, wahlen wir als Bezugs- 
geschwindigkeit 


2 
c—«ad, Gates G7) Rolce, pe ie U4=—Uy, v=6Rer,.: 


Nach Elimination des Stérdrucks #, und der Geschwindigkeitskomponenten 
w, fiihrt dies auf die Gleichungen 


(2.1) here u(0) = v(0) ms) 
(Ave OA Ven (1) 0) 
(2.2) L,fv=—s0 U4, 
mit 
= int a a eC i oe — 0". 


Die Randbedingungen fiir v’ ergeben sich aus denen fiir w nach der Gleichung 
(2.3) DV =O =, 


die aus der Kontinuitatsgleichung (1.4) flieBt. 


Damit ist unser Problem in diejenige Form gebracht, an die sich die weitere 
mathematische Behandlung anschlieBen soll. 


Die Gleichungen (2) finden sich im wesentlichen schon bei DEAN [2], und 
Yiu und SANGSTER behandeln ebenfalls dieses System fiir den neutralen Fall 
t=o. Dort wird die Variable =1—y eingefiihrt, und einige Bezeichnungen 
sind anders gewahlt. Ein Zusammenhang mit dem friiher behandelten Problem 
der Stabilitat einer Grenzschichtstrémung langs einer konkaven Wand ([5], [6]) 
14Bt sich nun leicht herstellen: Die Storungsdifferentialgleichungen sind hier 
dieselben wie dort, jedoch haben sich die Randbedingungen insofern gedndert, 
als an die Stelle des 4uBeren Randes der Grenzschicht (y= co) nun der innere 
Zylinder (7 =1) tritt, und entsprechend hat man hier natiirlich ein anderes Grund- 
profil U(y). 


ATS 
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3. Die Vorbereitung des Eigenwertproblems 

Die Gleichungen (2.1) und (2.2) stellen ein Eigenwertproblem dar, dessen 
Parameter o, t und S sind. Diejenigen Wirbel, die den Grenzzustand zwischen 
gedampften und angefachten Stérungen darstellen, bezeichnen wir als neutrale 
Stérungen; sie sind durch B =0 (t =o) gekennzeichnet. Sie interessieren beson- 
ders, da sie die kritische Kurve des Parameters S liefern, wenn wir von der 
iiblichen Auffassung ausgehen, daB zu einer vorgegebenen Wirbeldicke (beschrie- 
ben durch o) der zugehérige Wert S zu berechnen ist. Wir werden uns 1m folgenden 
auf diesen neutralen Fall beschranken. 


Fiir t=o lauten die Gleichungen (2) 


(3.1) Lu=(1—2n)v UO) == 210) === 90) == 0, 
(3.2) L?u = — So? (yn — 7°) u aly Sv = te 0, 
a 2 
See 


Die Aufgabe besteht darin, den kleinsten Eigenwert S in Abhangigkeit von o 
zu finden. Dazu ist zu bemerken: Fiihren wir an Stelle von 7 die Koordinate 
7=1—7 ein, so gelangen wir zu dem System (4#=u, v= —v) 


mit den Randbedingungen 


Dieses System unterscheidet sich von (3) lediglich durch das Vorzeichen von S. 
Hat also (3) fiir den Eigenwert S=-+ A eine Lésung, so hat das transformierte 
System eine Lésung fiir S = — A, in anderen Worten: Mit S=-+ A ist stets auch 
S=—A ein Eigenwert, oder es treten stets zwei Eigenwerte gleichen Betrags 
auf, die sich nur durch das Vorzeichen unterscheiden. Die Ursache dieser Er- 
scheinung ist in der Willkiir zu suchen, die in der Wahl der Variablen 7 bzw. 7 
liegt. Es wird sich zeigen, daB die Eigenwerte reell sind, und fiir unser Problem 
ist immer der positive maBgebend. Jedoch hat das Auftreten eines zweiten 
Eigenwerts gleichen Betrags eine noch zu erlauternde Bedeutung fiir die Berech- 
nung der unteren Schranken und fiir das Iterationsverfahren, welches ange- 
wendet werden soll. 

Zur weiteren Behandlung verwandeln wir das System (3) in ein System von 
zwei Integralgleichungen. Das geschieht mit Hilfe der bereits von H. WittinG [7] 
berechneten Greenschen Funktionen MW (yn, €) und N(y, €) der linken Seiten von 
(3). Sie haben die Gestalt 


M(n, £) = — a sinhoysinho(—1) in 7<é 
und M(é,7) = M(n, é), 


N(n,§) = {IM (0) + N2(0)] fe(o (€ — 1)) —o N (0) A(olE — 1))} fa(on) + 
+ {LN (0) = 2 (0)] Alo € = 1)) +6 Nj (0) flo (E —1))$ (on) fiir 4 SE 
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und N(n,&) = N(é,7). Dabei sind 


Nyojeels aa Sis fi(y) = ysinh y 


20 o*%sinh?o — o4 


Ni, (seer base f(y) = sinh y — ycoshy. 


2 o*sinh?o — ot 


Mit Hilfe dieser Greenschen Funktionen entsteht aus (3) das System 
1 
(4.1) u(n) = — J M(n, €) (1 — 2) o(€) dé, 
0 


(4.2) o() = Sot Neq.8) (€— 8) ule) a8. 


Das System (4) erlaubt einige Abschatzungen fiir S (o?), die der tatsadchlichen 
Berechnung vorausgehen sollen. 


4. Untere Schranken fiir die kritische Kurve 


Aus den Integralgleichungen (4) 1aBt sich durch Einsetzen von (4.1) in (4.2) 


eine einzige Gleichung herstellen: 
aL 


(5.1) On) SIE) vO) dt. 


Dabei ist der Kern K(n, t) dieser Gleichung von der Form 


1 
(5.2) K(y, t) = — (4 — 22) J o® N(n, £) (& — &) ME, 2) dé. 
0 
Der Kern K(y, #) ist stetig und nicht symmetrisch. 
Nach J. ScHUR (HELLINGER-TOEPLITZ [8], S. 1550) gentigen die Eigenwerte 
A, einer Integralgleichung 


f(x) = AS K(x, 9) f(y) dy 


mit beliebigem quadratisch integrablem Kern K(x,y) der Beziehung 
& a Fee 

Pelee ree age net ae 
» a <f [\K@»| dx dy; 


dabei tritt jeder Eigenwert so oft auf, wie seine Vielfachheit angibt. Entsprechend 
gilt fiir die iterierten Kerne 


ial 
(6) De ae <f [|KO )Paxdy. 
al t ene 


Diese Beziehungen liefern die in unserem Falle besonders wiinschenswerten 
unteren Schranken fiir den kleinsten Eigenwert der Gleichung (5). 

Im 3. Abschnitt wurde gezeigt, daB je zwei Eigenwerte gleichen Betrags auf- 
treten. Damit ist S,—-— S,, und die aus (6) folgende Abschatzung fiir S, ver- 
scharft sich zu 14 
ae =| |K™ (, |? dy dt. 


00 


218 GUNTHER HAMMERLIN: 


In dieser Abschitzung bereitet die Berechnung der iterierten Kerne gewisse 
numerische Schwierigkeiten. Zundchst ist nach (5.2) der Kern A(y, 2) = Kk (n, t) 
zu gewinnen; diese Integration iiber & wurde mit den Schrittweiten /=0,2 und 


h—0,1 durchgefithrt. Daraus ergibt sich, wiederum durch Anwenden der auf 
a iat b 


Gebiete iibertragenen Trapezregel, {fK?dydt, und damit die Abschatzung 


00 ee 
y2 


{ f fK2(n,t) dnat\* 
00 


Da wir numerisch integrieren, ist die Angabe des Doppelintegrals mit einem 
gewissen Fehler behaftet. Eine iiberschlagige Fehlerbetrachtung zeigt, daB die 
Angabe von E,(c) in Fig. 2 (gewonnen aus den Werten fiir h=0,2 und h=0,14 
durch Extrapolation mit der bekannten  tiber- 


] 3 : 
§-19°§ Saas schlagigen Uberlegung) sicher auf 1% genau ist. 


| Sy 2 hy= 


E, (0) Ebenso ergibt sich mit dem iterierten Kern 


1 


| K® (n,t) =f K(n, w) K(u,t) du 
mer 2 


die Abschatzung 


ba 
Spe foe Fe = : les 
: ; LKGA(n, t tt 
Gin tae) oe whe ORO (Jf* (9, t)dn at, 
Fig. 2. Untere Schranken fiir |S,(o)|. Peres ; : 
Be pi ones 2is ale ‘l- Die in Fig. 2 angegebenen Werte fiir E, sind auf 
3% genau. 


Eine weitere Verbesserung dieser Abschatzungen mit Hilfe der hdheren 
iterierten Kerne ist méglich. Dabei ist jedoch zu bedenken, daB der Integrations- 
fehler die Angabe héherer iterierter Kerne in steigendem Mae unsicher macht 
und daB auBerdem die zu erwartende Verbesserung der unteren Schranke von 
Schritt zu Schritt geringer wird. Wir begniigen uns deshalb mit der Angabe der 
Schranke E,(o0) — die schon sehr brauchbar ist — und gehen zur Berechnung der 
Eigenwerte iiber. An und fiir sich macht das genaue Verfahren, das im folgenden 
verwendet wird, die Angabe unterer Schranken iiberfliissig. Sie sollen in diesem 
Fall auch nichts mehr bringen als eine vorlaéufige Aussage itiber die GréBe des zu 
erwartenden Eigenwertes S,. 


5. Berechnung der Eigenwerte kleinsten Betrags 
Die Berechnung des kleinsten Eigenwerts S, soll ihren Ausgang vom System 
(7) nehmen. Wir verwenden dazu ein Iterationsverfahren, das bereits in [6] und 


spater in [7] erfolgreich durchgefiihrt wurde. Die exakte Fundierung des Verfah- 
rens findet sich bei H. WieLanpt [9]. 


Dabei gehen wir von der Anfangsfunktion v© (&)=16(&—£)? aus, die die 


Randbedingungen v (0) =v'(0) =0 und v(1) =v'(1) =0 erfiillt, und berechnen die 
weiteren Iterierten nach der Vorschrift 


(7.1) u(y) = —f M(y,8) (1 — 2€) (6) dé, 


(7.2) v(m) i 0? N(, &) (E — &) ul(€) dé. 
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In der Arbeit [6] waren die iterierten Funktionen durch elementare Inte- 
gration ermittelt worden; das ware auch hier mdglich, jedoch ziehen wir jetzt 
der Einfachheit halber eine numerische Integration vor. Dazu wird durchgangig 
die Trapezregel verwendet, die in (7.1) mindestens teilweise unumganglich ist, 
da M(n, €) als Greensche Funktion einer Differentialgleichung 2. Ordnung einen 
Sprung in der ersten Ableitung besitzt. Aus Griinden der einheitlichen rechneri- 
schen Durchfiihrung wird auch in (7.2) mit ihr gearbeitet. 

Wahrend eine Iteration durch elementares Integrieren den exakten kleinsten 
Eigenwert mit beliebiger Genauigkeit liefert, diirfen wir das jedoch hier nicht 
mehr erwarten. Die Ersetzung der Inte- 


grale durch endliche Summen bedeutet ohe? d 5,108 ; 

ja eine Ersetzung des Systems (4) durch | 3 450} t 
ein System linearer Gleichungen, das ik 

dann iterativ gelést wird. Durch die Ite- 42. 


ration wird nun die Lésung dieses letz- 
teren Systems beliebig genau geliefert, | 72 
und es hangt von der gewahlten Schritt- 4°” |— 9% 


weite ab, wieweit diese Lésung mit der a7. GE 
des Systems (4) bzw. (7) tibereinstimmt. 
Wegen des Auftretens von je zwei G50 T 
Eigenwerten gleichen Betrags haben wir | 
mit dem Eintreten des Oszillationsfalls O25 
der Iteration zu rechnen ([9], S. 126). | 
Dann approximieren zwar die Iterierten ar Oe A ca aS ET 
noch nicht die gesuchten Eigenfunk- Fig. 3. Kleinster Eigenwert S,(c); kritische Kurve: 
tionen, aber diese stellen sich als Linear- gine Bis. 


R, 36 
kombination zweier aufeinanderfolgen- 


der Iterierter dar. Betrachtet man etwa die Iterierten wu” und ut, so werden 
die zum kleinsten Eigenwert S, gehérenden Eigenfunktionen durch 


(8.1) ure uM +S yerh) 
und 
(8.2) vw vl) 4S ytd) 
approximiert. 
Aus demselben Grund wird hier nicht der Quotient gigi linearer Funk- 
tionale gegen S, konvergieren, sondern es geht wire —> Si. 


Die Rechnung wurde mit den Schrittweiten h=0,2; h=0,1 und h=0,05 fiir 
o=1, 2, 3, 3.4,4,5,6,7 durchgefiihrt. Die in Fig.3 angegebene Kurve S,(o) 
ergibt sich daraus durch quadratische Extrapolation der sehr schnell konver- 
gierenden Werte fiir 0. 


6. Die Eigenfunktionen und das Aussehen der gestérten Stromung 
Wie oben schon erwahnt wurde, liefert das Iterationsverfahren nicht nur den 
kleinsten Eigenwert, sondern auch die zugehérigen Eigenfunktionen. Die Kennt- 
nis dieser Eigenfunktionen gestattet es, das Aussehen der entstehenden Wirbel 
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zu beschreiben. Die gefahrlichsten unter den neutralen Wirbeln sind diejenigen, 
zu deren Wirbeldicke der kleinste der im 5. Abschnitt errechneten Eigenwerte 
gehort. Das Minimum der Kurve Sj(o) liegt beio=4. o=4 bedeutet, daB die 


Dicke eines solchen Wirbels a = ae =0,79d ist. In Fig. 4 sind die Eigenfunk- 


\ 
G2 0% G6 G8n 40 


Fig. 4. Die Eigenfunktionen fiir die 

kritische Wirbeldicke A/2=0,79 d 

(o=4). n=0 entspricht r= R, (auBe- 

rer Zylinder), 7 =1 entspricht ry = R, 
(innerer Zylinder) 


tionen fiir o = 4 angegeben. u(y) und u(y) ergeben sich 
nach (8), wahrend sich w(y) aus der Beziehung (2.3) 


OY, 
w(y) = '(n) 
errechnet. 


Die Komponenten der gestérten Strémung sind 


u(y) + u(y) cosaz, v,(y)cosaz und w,(7)sin«z. 


Uber die tatsachliche GréBe der Wirbelkomponenten 
kann nichts gesagt werden, da die Eigenlésungen u,v, w, 
die den Stérgeschwindigkeitsamplituden w,, v,, wy, 
proportional sind, als Lésungen eines homogenen 
Systems nur bis auf einen Faktor zu bestimmen sind; 
auch die relativen Groé8enverhaltnisse von uv, v und w - 
entsprechen nicht denen von u,, v, und w,, da diese in 
der dimensionslosen Darstellung mit verschiedenen 
Faktoren versehen wurden. 


Dennoch erlauben es uns die Eigenfunktionen, die 
Gestalt der entstehenden Wirbel anzugeben. Fig.5a 


zeigt einen Querschnitt senkrecht zu den Wirbelachsen, wobei der Betrachter 
in Richtung der Strémung blickt. Fig.5b zeigt die Profilformen in Strémungs- 
richtung, also (7) =%)(y) +-%() cos az. Die ebenen Schnitte senkrecht zu 
der Achse der Kanalberandungen sind dabei in den folgenden Héhen gelegt: 


8 = Uy + U, COS aa z Uv =v, Cos as | o= Ww, sin = by 

RA Up + % VY; ) 
shit, ; : 

4 0 Wy 
gist A Uy — U = 

4 0 1 Vy O 
4h+3 | , : 

4 0 ap a 


B= Oy aA eee 


_ Zur Zeichnung der Profile w(y) ist angenommen, daB Max |1(n)|= fy Max |u(7)| 
sei. Die angegebenen Wirbel sind diej enigen der Dicke 4/2 =0,79d. 


Die Kerne dieser Wirbel sind etwas nach der 4uBeren Kanalwand hin ver- 
schoben; diese Erscheinung riihrt von der stabilisierenden Wirkung der konvexen 
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inneren Kanalberandung her. Die Grundprofile werden abwechselnd nach innen 
und nach auBen ausgebeult, wahrend die Stérkomponente uw, in der Héhe der 
Wirbelkerne verschwindet. 


rahe 


Fig. 5a u. b. Das Aussehen neutraler Wirbel. a Schnitt senkrecht zu den Wirbelachsen; 
b Strémungsprofile in Strémungsrichtung (g-Richtung) 


7. Vergleichende Betrachtungen 


Es besteht die Notwendigkeit, noch ein Wort zu den numerischen Ergebnissen 
im 5. Abschnitt zu sagen. Der Eigenwert S(c) hatte die Bedeutung 


_ 79(4\? a 
Ss Pelee 


Bei den Taylorschen rotierenden Zylindern erscheint ein entsprechender Para- 

meter p= 2Retat und bei der Untersuchung von Grenzschichten langs kon- 
1 

kaven Wanden tritt w= 2RePS auf. Die Reynoldsschen Zahlen sind dabei 


jeweils auf die maximale Grundgeschwindigkeit und auf Spaltbreite d bzw. 
Grenzschichtdicke 6 bezogen. Wegen |#)| = 3 Max |u| = $ uy, ist 


Uy d 2d 1 
( v Ey 16 


Mit diesem Wert wird ein Vergleich moéglich. In der Umgebung des Minimums 
hatten sich folgende Werte ergeben (die angegebenen Ziffern sind als sicher zu 
betrachten) : 


Uy da 24d pe ad (sac) a 
o 7) R, 2 oan: 2 is R, 
| 
3 1,014-10° 2,82 - 108 | 6,34 - 10° 
0,918 - 10° eeu OLe | Weukeds ee 
5 0,990 10° Neee2e75 3405" |) oe 0,10.10° 
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Das Minimum von S,(o) liegt bei o=4 und hat den Wert 


2d 
Fie 5740. 


2 . ,(umd 

Simin = 91800 oder 2,{ 2m =) 
W. R. DEAN [2] gibt hierfiir den Wert N=36(4) an; das ware in unserer 
d 
oe 
dieser Wert nur mit einer Sicherheit von etwa 2,5% gilt. An und fiir sich ist 
diese Genauigkeit ausreichend, und wir stellen fest, daB unser — genauerer — 
Wert 5740 innerhalb der Genauigkeitsgrenzen des Deanschen liegt. In diesem 
Sinn haben wir eine Ubereinstimmung des kritischen Werts mit der friiheren 
Angabe von W.R. Dean. Dasselbe gilt fiir die kritische Wirbeldicke, die wir 
mit A/2=0,79d (c=4) angeben. Der von Yiu und SANGSTER [3] angegebene 
kritische Wert S;yjn—=600 wird als falsch erkannt. 


Zu den Methoden, die in den erwahnten Arbeiten verwandt werden, ist 
folgendes zu sagen: YIH und SANGSTER setzen die Eigenfunktion w als Fourier- 
Reihe an, lésen dann (3.2) und setzen das so gewonnene v in (3.1) ein; homogene 
Bestimmungsgleichungen fiir die noch unbekannten Koeffizienten des Ansatzes — 
ergibt der Vergleich der Koeffizienten der nach trigonometrischen Funktionen 
entwickelten Gleichung (3.1). Die Determinante dieses Systems liefert, gleich 
Null gesetzt, eine Bedingungsgleichung fiir den gesuchten Eigenwert. Die Ver- 
fasser beriicksichtigen von dieser Determinante unendlich hoher Ordnung nur 
ein Glied und berufen sich dabei auf CHANDRASEKHAR [10]. Tatsachlich fihrt 
jedoch CHANDRASEKHAR, der die Taylorschen rotierenden Zylinder so untersucht 
hat und zu brauchbaren Approximationen gelangte, den Vergleich in anderer 
Weise durch. Hinzu kommt, da8 die Beschrankung auf ein Glied bei Yru und 
SANGSTER letztlich darauf hinauslauft, die Eigenfunktion «(7) durch C, sin 77 
zu ersetzen. Wie Fig. 4 zeigt, besteht damit kein Grund zu der Hoffnung, so zu 
einem brauchbaren Resultat zu gelangen. 


Bezeichnungsweise 2 (Huey 5832; zu beriicksichtigen ist allerdings, daB 


Es ist noch der Bemerkung wert, daB bereits DEAN [2] mit Fourier-Ansatzen 
fiir die gesuchten Eigenfunktionen arbeitet, die ebenfalls auf eine unendliche 
Determinante fiithren. DEAN treibt die Annaherung dieser Determinante durch 
solche endlicher Ordnung jedoch weit genug (bis zu 7reihigen), um seine numeri- 
schen Angaben sichern zu kénnen. 


Wir kénnen unser Resultat noch mit dem kritischen Wert einer laminaren 
Grenzschichtstrémung vergleichen. Fassen wir die halbe Kanalbreite als Grenz- 
schichtdicke 6 =d/2 auf, so rechnet sich der kritische Wert der Kanalstromung 
Un } 


um zu 2( sep ae! ha Fiir laminare Grenzschichten ergab sich in ent- 


Min 

sprechender Approximation bei parabolischem Grundprofil der Wert 80, der 
betrachtlich niedriger liegt. In der Grenzschicht waren eben gerade diejenigen 
Wirbel die gefahrlichsten, die eine groBe Dicke hatten und auch weit iiber die 
Grenzschicht hinausreichten. Das Vorhandensein der inneren Wand verbietet 
bei der Kanalstroémung solche Wirbel; wir haben keine Grenzschicht im obigen 
Sinne, sondern diese wird von der konvexen Wand her stabilisierend beeinfluBt. 
Der kritischen Wirbeldicke 2/2 =0,79d entsprach im Fall der Grenzschicht ein 


v 
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kritischer Wert 2Re?6/R=275, der, verglichen mit 717, ebenfalls die stabili- 
sierende Wirkung der inneren Kanalwand zum Ausdruck bringt. In unveréffent- 
lichten Rechnungen fand H. Wirtinc fiir die Sekundarinstabilitat einer ebenen 
Kanalstr6émung den kritischen Wert 2 Re?6/R = 822, wobei mit 6 die halbe Kanal- 
breite bezeichnet ist, also einen Wert, der in der Nahe von 717 hegt. 


Zusammenfassung 


Die Wirbelinstabilitat einer laminaren Strémung durch einen gekriimmten, 
unendlich hohen Kanal, dessen Breite klein ist gegeniiber den Kriimmungs- 
radien der Kanalwande, wird durch die beiden Differentialgleichungen (2) be- 
schrieben. Das System (2) stellt ein Eigenwertproblem dar, bei dem vor allem 
die kleinsten Eigenwerte S, in Abhangigkeit von dem Parameter o mit den zu- 
gehorigen Eigenfunktionen interessieren; dabei beziehen wir uns bereits auf (3), 
wo neutrale Stérungen (t=o) betrachtet werden, und S,(o) ist bis auf einen 
Zahlenfaktor die kritische Kurve des Parameters 2 Re? d/R, in Abhangigkeit von 
der Dicke 4/2 der angesetzten Wirbel. Mit Hilfe Greenscher Funktionen werden 
die Differentialgleichungen (3) in die Integralgleichungen (4) verwandelt, die es 
zunachst erlauben, untere Schranken fiir S,(o) anzugeben (Fig. 2). Das Itera- 
tionsverfahren nach [9] liefert S,(o) mit hinreichender Genauigkeit nach wenigen 
Iterationsschritten (Fig.3). Ebenso ergeben sich damit die Eigenfunktionen 
(Fig. 4), also die Wirbelkomponenten, so daB in Fig. 5 das Aussehen der neutralen 
Wirbel angegeben werden kann. Der kritische Wert — aus dem kleinsten unter 
den Eigenwerten S,(o) — ist 2Re?d/R,=5740, wahrend DEAN [2] 1928 den 
ebenfalls richtigen Wert 5832+2,5% gefunden hatte. Die zugehorige kritische 
Wirbeldicke ist 4/2 =0,79d. Die Ergebnisse von Yiu und SANGSTER [3] erweisen 
sich als falsch. Die Arbeit verfolgte zwei Ziele: Eine Neubehandlung der Diffe- 
rentialgleichungen (5) mit exakt begriindeten Methoden, um eine Beurteilung 
der Ergebnisse von DEAN und von YIH und SANGSTER méglich zu machen. 
Dann die Angabe der Eigenfunktionen, die es gestatten, das Aussehen der ent- 
stehenden Wirbel zu beschreiben. 


Diese Arbeit wurde durch das Wirtschaftsministerium des Landes Baden- 
Wiirttemberg gefordert. 
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Determination of the Vorticity 
and the Gradients of Flow Parameters behind a 
Three-Dimensional Unsteady Curved Shock Wave 


R. P. KANWAL 


Communicated by C. TRUESDELL 


1. Statement of the problem 

The differential effects of shock fronts have been the subject of several 
researches of gradually increasing generality. My recent papers on three-dimen- 
sional flows of perfect gases calculated the jumps in the pressure, density and 
velocity gradients across shocks that are stationary [J, 2] or pseudo-stationary [3]. 
The determination of vorticity in these two cases led to characterization of shock 
surfaces which may join two regions of irrotational flow and, in the latter case, 
to a theorem regarding the mode of propagation of such shocks. 

By a different method, LIGHTHILL |4] derived an expression for the vorticity 
behind a steady shock in a gas obeying an arbitrary equation of state. HAYEs [5] 
remarked that for plane shocks this result is contained in an earlier analysis of 
TRUESDELL [6], who concluded that the magnitude of the vorticity generated by a 
shock of given strength and curvature depends only on the magnitude of the tangential 
component of the velocity and is independent of the form of the equation of state. 
The strength of the shock is defined as the density jump divided by the density 
in front. As Hayes observed, TRUESDELL’s statement suggests that the result 
is purely dynamical and thus should be derivable without recourse to any thermo- 
dynamic law, and Hayes has provided a proof based on the momentum equation 
alone. He did not mention, however, that this same independence may be inferred 
also from TRUESDELL’S analysis. The reason the energy equation appears there 
is that, as in earlier work of THomaAs [7] and in my papers cited above, it is 
used to calculate the jumps of the full set of velocity gradients. By following 
the details in all these derivations of more complete results [J, 2, 3, 6, 7], it is 
easy to verify that in them, too, the energy equation is not used to calculate the 
vorticity Jump. 

The purposes of this paper are 

1. To calculate the jumps of pressure, density, and velocity gradients in full 
generality, that is, for possibly unsteady shocks in fluids obeying an arbitrary 
equation of state. As usual, dissipative mechanisms such as viscosity and heat 
conduction are assumed absent. 
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2. To arrange the analysis, which follows and extends the methods used in 
several earlier papers [7, 2, 3, 6, 7], in such a way as to make plain the independence 
of the jump in vorticity from the equation of energy and from the form of the equation 
of state. In particular, the following results show that if a shock wave were 
possible in an inviscid but thermally conducting gas, the vorticity it generated 
would be independent of the thermal conductivity. 

For the case when a curved shock is propagating into a quiet region, the 
results reduce to a simpler form. In parallel propagation, when the normal to 
the shock is steady, as is the case for example for spherical waves having a 
common center, the flow behind the shock remains irrotational. 

The paper ends with a simpler derivation of the extended geometrical and 
kinematical conditions of compatibility recently discussed by THomas [8]. 


2. Notation and assumptions 


As in previous works [/, 2, 3], we assume the shock surface 2(¢) represented 
by a continuously differentiable function x;=%;(y1, y?, #), where the x; are 
rectangular Cartesian coordinates in an inertial frame and where the y% are 
coordinates on L(t). The range of Latin indices, referring to spatial tensors, is 
1, 2,3 that. of Gréek-indices,, referring to surface tensors, is.f, 2. .The unit 
normal €;, assumed continuously differentiable, is directed downstream. 

A quantity f if evaluated on the upstream side of the shock will be denoted 
by /,; if on the downstream side, by /; and [f/f] =f—f. 


3. Basic equations expressing balance of mass and momentum 


The differential equation expressing the balance of momentum in an inviscid 
fluid subject to no extrinsic force is 


Ou; 
(1) lary + 0Uu, 4, + p,;=0. 


Let G denote any speed and &; any direction, and let 5/d¢ denote the time 
derivative as apparent to an observer moving with velocity Gé,, viz 


ley 


Py en) A 

2 oi hae iS 

(2) Oe. ot + Ghiré;. 

Then the momentum equation (1) may be put into the form 
Ou; 

(3) 0 ra THO Uy (Ui GE ots Peete th, 


In what follows, we use this result with G &, taken as the velocity of propagation 
of the shock. 


The strength 6 of the shock is defined by 


(4) jel 
Q1 
Write 1¢,==1, 655,21, 6. Jump conditions expressing conservation of mass 
and momentum are 
m) 
(5) [4] =— aa) (1, — G) &;, 


(6) [1 = 75 a1 (1n— 6) 
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Quantities expressed as functions of 0,, £,, “1, &;, G, 6, the metric tensor of 
the shock surface, and their derivatives along the shock surface we shall call 
effectively calculated. Such quantities depend only upon the shape, speed, and 
strength of the shock and upon flow conditions on its upstream side. The objective 


of differential shock analysis is effective calculation of the jumps in various 
gradients. 


4. Properties of shocks following 
from the balance of mass and momentum alone 


Differentiating both sides of (5) and (6) with respect to the coordinates on 
the shock surface, we have 


m) 

(7) Uj; Xiu = Uy, Xia {3 (u,, — G) = AS 
a) 

(8) Pri Xin = Pr i%jat [apy an 6)? = BE, 


where %;,=0%x,/0y%. The quantities Aj, and Bx are effectively calculated. 


From (3) and (8) we eliminate #, ; and obtain 


5 Be Ou; 
(9) Uj, 5 (Me; i G §) Nore a aS ve Nic 
We set 
u;—G&; 
(10) Oe ze en ae : Go % unsummed, 
w 


where g, = (dxp)~%, %g being the surface metric. From (5) it follows that C;; 
is effectively calculated. The spatial unit tangent vectors to the surface coordinate 
curves are %;, 8 q, « unsummed. Thus det C;;=1; therefore the matrix C;; has 
an inverse D; ;: 


(11) D;; = cofactor of C;; in Coe i; D;, Cy; = 9; Dj; Cj, = 943. 
Now set 
(12) rere A gare 


By (11) we may invert this relation and obtain 


(13) Chips B;;D;,D; 


jm 
From (11), (10) and (5) we see that D;; is effectively calculated. Thus if B; 
is effectively calculated, by (13) it follows that u,,, is effectively calculated. 
Now from (12), (7) and (9) we find that 


7 


Byg = Alp %ta8aakep  % B unsummed, 


(14) Day Y [pees % unsummed, 


\ 


ze Be Ou; 
Pas aE on boar) 
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Thus the components B,,, Byzg and B;, are effectively calculated. For B33, 
however, we obtain only 


(1 5) B33 = Um — (Un =Cr , 


and this does not provide effective calculation of B35. 


5. The vorticity jump 
By (13) and (11) we calculate the vorticity behind the shock: 


(16) WO, = C7 n Un, 7 = Fgh Bim Dir Dajg= —Cx Cemp Bee 


Thus the vorticity is determined by B,,,—B,,,, and since these quantities are 
effectively calculated by (14), it follows that the vorticity behind the shock 1s effec- 
tively calculated. For steady or pseudo-steady shocks, both the method and the 
result reduce to those obtained in the earlier investigations [1, 2, 3, 6] described 
in §1, and thus it is clear that the energy equation and the equation of state are 
not required for determining the vorticity jump across a shock of a given strength. 


6. Use of the balance of energy to determine B,, and 6 


In order to calculate effectively the full set of gradients w, ; behind the shock, 
we need to determine B,,, and to this end we introduce an equation of state 
f(€, @, n) =0, where ¢ is the specific internal energy and 7 is the specific entropy. 
This equation is assumed soluble for any of its variables as a function of the 
other two, and the resulting equation is assumed continuously differentiable. 
The pressure # occuring in (1) and (6) is assumed to be related to the equation 
of state as follows: = — de/0(1/0), where 7=const. The speed of sound, c, 
is defined formally as c?=(0f/0e),. We assume further that the fluid does not 
conduct heat. When (14) and the equation of continuity are satisfied, the equation 
of energy is equivalent to 
(17) oO gu; 98 — uu; mi, + 20 uy, = 0. 


By use of (2) we may put this condition into the equivalent form 


(18) em; s(4j — GE) (4% — GE) =P — 9 U8 (4, — GE) + oman. 
The quantity u, , may be eliminated by (13): 
(19) i a B;; Dj, D5, 
= By (1+ 9a) + Boa (1-+48) + 2 Boe x1 %2— 2Basym1 — 2 Bees Xe +Bes, 
where B;;=3(B;;+ B,,), and 
(20) —— nee % unsummed, v 


From (15), (18) and (19) we observe that Bg, is effectively calculated. We | 
have thus determined all the components of the tensor B;;. | 
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Moreover, conservation of energy at the shock implies 


d{1+36 
(21) [A] = I (in — 6 
where h is the enthalpy, h=e+ 4/0. Since h=h(p, 0), a formula for 6 as a 
function of (“,,,—G), p;, @,, may be obtained in principle by elimination among 
equations (21), (4) and (6). For the special case of a perfect gas with constant 
specific heats, there thus results 


22 es 2{01(Uin — G)?— yPy} 
( ) 2¥ Pi + (y—1) @1(%1n — G)?’ 
where y =¢,/c,. 


7. Evaluation of the pressure and density gradients 


So as to calculate the pressure gradient behind the shock wave, we substitute 


for u;; in equation (3) and obtain 


Ou; 
(23) P.i = — @ (4, — G) B;3D;;—e—4. 


To evaluate the density gradient we differentiate both sides of the equation 
(4) with respect to the coordinates y* and get 


(24) 0: Xia = 01,4 %ia + (60)).0= CF. 


The equation of continuity gives 
(25) =~ + 054 (u% — G &;) + OU, ,=0. 
The equations (24) and (25) can be written as 


(26) 0 


where 


Hence by inversion we have 


8. Calculation of the quantities B;; and d; 

As in the case of steady and pseudo-steady flows, it is possible to calculate 
explicitly the quantities B;; and d; in the case of unsteady flow of a perfect gas 
with constant specific heats, if we assume that the flow upstream of the shock 
is uniform. For simplicity we take the lines of curvature as the coordinate curves 
on the shock surface. Thus 


(28) Ae=Oy,.=0; Ky= dy Jay, Ky = dy9/420, 


where b,, are the components of the second fundamental form and ky and K, 
are the principal normal curvatures of the shock surface. WEINGARTEN’Ss formulae 
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assert that 


a vB eae 
(29) Ey = — Ob Aggy a a yas % unsummed, 


where a*? are the components of the tensor reciprocal to a,g. We further have 
the relation [5, 8] 
0; 
(30) ar a a“’G,, Nip 
Keeping these results and the assumption of uniform flow upstream of the 
shock in mind, we get 


a, ore BE gi : Lo] ra 
(31) Ain= Loti Gh me 0; (u,,— G)? -( (Uy » G) K a Kia + &; (Vy K me) 
32) Bt SON AEG) Ws Deal 
y+1 
(33) Gia 4otyv(y +1) (“in — G) (Va Ka + G, «) Pr 
: v {2y p+ + (y —1) 01 (41 — G)?}? : 
ou; Op E; n) a c | 
ee ee eee eee ee), 
. dG 
(35) Cn 401 (4 G) (v Gat} 
ot y+1 : 
F ] \G 
6 Oo 4 ot y(y + 1) (4— G) (Gato), 
aes Brkt O=1 alan O 


The quantities B;; and d; can now be readily calculated and are given as 
By, = 0(u, —G) K,, 
By, = By, =0, 


Byg= 0(u, — G) Kg, 


(4 Ay + G1) 
ES 


ay , 


By, = 


if 
(37) Bey= 4a, ve “ait Ga) _ 3G,,| 
mes es 
eI) 
ae 


\ 
J 
iT + 6K, — 25, (Ke+G,)h, 
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and : 
d,=— 44, 8 V(Y +1) (Hin —G) (Ky +EG; 1) Pr 
{2y Pi + (vy — 1) @ (Uin — G)?}? é 
(38) Ti 4ay0 01 (y +1) (Uen—G) (v2. +G, ») Pr 
" {2y Pi + (y — 1) Q1 (Min — G)?}? : 


0G 
‘ / 40 y(y + 1) (41, — G) (*6,.+ =} Pr 
Paaey. Ug — G e “hy k be {2y Pit Op 1) 01 (Mn — G)?}2 : 


where the quantity «, , entering the expressions for B,, and d, is given by the 
equation (19). 

Thus when the flow upstream of the shock is uniform the expression for the 
vorticity vector behind the shock is 


0? 2 5 
(39) C= 411422 {— (v, Ky + G, ) %j2 + (Up Ke + G, >) He}, 


which agrees with HAYEs’ result. 


Now if the shock propagates into still air, the values of B;, and d; take a 
very simple form. In that case we notice from the expressions (37) that im the 
parallel propagation of shock waves, i.e., the normal directions to the shock surface 
d'(t) at t=ty remain normal to X(t) for t>ty, e.g. in the case of plane parallel 
shocks or spherical shocks having common center, the flow behind them remains 
irrotational., 


Having determined the gradients of flow parameters we can determine their 
variation along the shock wave as was done in the stationary case [2]. 


9. Derivation of Thomas’ 
“geometrical and kinematical conditions of compatibility” 


Finally we give a simpler derivation of the “first extended compatibility 
conditions” defined by THomas [8]. In fact, if we are given a jump relation 
(40) [Zl =A 
and we differentiate both sides of this relation with respect to y%, we get 
(41) SEN WS cee 
Multiplying both sides by a%’x,, and using the relation a%?x;,%;,—=0;; —§;§; 
yields 
(42) [Z,,] —=41[Z, 16) é-- a A, Xe: 

This is THomas’ “‘first geometrical condition of compatibility”’. 

Applying the operator 6/6¢ to both sides of (40), we get 

OA 


(43) Se) —iZdahe+ Se. 


which is THomas’ “‘first kinematical condition of compatibility”. 
Tuomas’ “‘second conditions of compatibility”’ can be derived from the first 
ones, as he has remarked. 
10* 
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On the Convergence 
of the Rayleigh Quotient Iteration for the Computation 
of the Characteristic Roots and Vectors. I 


A. M. OSTROWSKI * 


1. Let A=(a,,) be an (n xm) real symmetric matrix. Then for a vector = 
eet ic,) we: put 


(1) Qu (€) =EAE’ r= dias Xu Xy. 
The quotient ms 
Qa (é) 
2) ANS) 
os [EP 


is called the Rayleigh quotient corresponding to &. If & is a characteristic vector 
belonging to a characteristic root A, then the corresponding Rayleigh quotient is . 
Therefore the following procedure has been devised for obtaining a sequence of 
numbers 4, (x=0, 1, ...) converging to a characteristic root: 

For any A,, of the sequence (%=0,1,...) find an approximate solution &, of the 
homogeneous system 


(3) At Ney 
and put 
(4) ee oe (Oe 


My attention was drawn to this method by Joun Topp, who used it in his 
lectures as long ago as 1945. It appears to converge fairly well in numerical 
practice. In what follows I give some theoretical results on the convergence 
of this method. 


2. The crucial point in the discussion of the above method is of course a 
suitable rule for the computation of the “approximate solution” &, of (3). 
The rule I shall use in the first part of this discussion consists in taking an arbi- 
trary vector 7 + 0 and in putting 


(5) &,= (4 -4,£) 11’. 


The theoretical arguments in support of this rule are given in another paper**. 


* In writing this paper I had very valuable discussions with Mr. Cur. BLATTER. 

xx A OsTROWSKI, ,,Uber naherungsweise Auflésung von Systemen homogener 

linearer Gleichungen“, Journal of Applied Mathematics and Physics (ZAMP), Basle, 
Vol. 8 (1957), pp. 280—285. 


234 A. M. OsTROWSKI: 


Since the formulae (4), (5) are invariant, for our discussion we can introduce 
normal coordinates from the beginning and therefore without loss of generality 
put 


n 


(6) Q4(é) = Lt ap 
where 
(7) My She So Sha 


are the characteristic roots of A, ordered increasingly. Then, if 


(8) TMV ey DS Os 
we have from (5) 
Se) ee ee a0 
E, = (% ) 7) i = on 
7s “ My Vv fe == a Bf 
Qa (Ex) ats (My — A)? : IE. ya (Ht, 4)" . 


and finally 


(9) 


ey 
x 
+ 
f 
| 
Ms |iM= 
| 
ze 
x 
! 
S 
— 


(4, —A,)? 


= 
I 
pany 


3. The expression on the right side in (9), if all products w, y; in the numerator 
are replaced by p,y? or ,, V2, reduces to 4, and w,,, respectively. We see that in 
any case 
(10) ieee Net ae 


In the expression on the right side of (9) a uw, drops out if the corresponding y, 
vanishes. Denote the remaining distinct wu, in increasing order by 


(11) 0, S Og S++ Oy. 
Then the formula (9) becomes 
™m On 
= (Git 
(12) Aye = Ge 0A a, 
py tae, 
(a, —h,)" 


w=1 


where the coefficients p, are all positive. 


4. Denote one of the o,, by o and the corresponding , by ~. By subtracting o 
from both sides of (12) we obtain 


m m 
2 eg ea Uy | kee ee 
Py Pu (6, in Ah,\e 2 Pr ( be rsA,\2 
(13) Agia iby teu at Pa =\(0 Ay) Mm guine2 
Bib ht Ly 6 (<=> 
pau Gia im p 2 < ot, 


where in the last sums the terms with the index w for which o,, =o are to be 
omitted. From this formula we see that if one A, gets sufficiently near to o, 
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the whole sequence 4, tends to o. Then dividing the first and the last term 
of (13) by (A, —o)?, we obtain 


(14) ee ee 
ret 


(A,,— 9)? Pema 


We see that in this case the convergence is at least quadratic. The convergence 
could even be faster than that, if the limit in (14) were 0. However, for C=O; 
or 6=o,, the convergence is exactly quadratic *. 


5. In order to characterize the convergence neighbourhood of o put 


(45) d= Min |o, —o|, Pipe 
o o =1 
If we then assume ; i 


(16) |4,-o] < 
for o,,==o it follows that 


lon 9] |o,—9| |o.—9| 4 Gs 
(o,=4,)7 (o,-—o-- o—/,)" a 


Introducing this in the numerator of the last term in (13) and replacing the 
denominator by #, we get 


pat Pu 
| Agel eA rae din Paep 
(Qi 0)" ede A Pp’ 
7 [Are —o| ZVI 
(17) ee ae fe ue 
AS 4 P—ph 
(18) jie = p | A, o| 


Therefore it follows that if for one index x we have 


(19) |A,=o| <$Min( ph 1), 
then the distances |A,—o| from that index on are strictly diminishing and 
therefore converge to 0. We see that a convergence neighborhood of @ is given 
by (19). 

Of course in this way we obtain only sequences /, converging to those charac- 
teristic values uw, of A which remain among the o,. If the constant vector 7 is 
orthogonal to all characteristic vectors corresponding to a characteristic root j,, 


* This result agrees with the note of WALTER Kouwn, “A Variational Iteration 
Method for Solving Secular Equations’, Journal of Chemical Physics, 17, 670 
(1949). In this note W. Koun discusses the application of the Rayleigh quotient 
method in taking one of the coordinate unit vectors for 7. He says then (in our 
notation): “‘A more careful analysis shows that 


61g = K,(04,)2 


where K, is in general of the order 1.” 
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then this characteristic root drops out, since all corresponding y? vanish. How- 
ever, we obtain for each choice of 7 at least one characteristic root, if the start- 
ing A, is chosen near enough to such a root. In particular cases it may be neces- 
sary to try out several choices of 7, for instance to try each of the coordinate 
unit vectors. 


6. In order to discuss the global convergence situation we have to consider all 
fixed points of the iteration A,,;—(A,), where 


AG 2 


(20) g(a) =* 


The corresponding algebraic equation for the fixed points, A= (A), becomes an 
equation of degree <2m—41. We know already m different roots of this equa- 
tion, given by the o,,, and we have seen that all these fixed points are points of 
attraction. By a theorem which we proved in another communication*, between 
two consecutive fixed points of attraction there is always at least one fixed point 
of repulsion. We see that besides the m fixed points (11) the iteration by @(A) has 
at least m—41 further different fixed points. Therefore the iteration by (A) has 
exactly the 2m—1 fixed points 


(21) Op U 6p Oe De 


where the #, are points of repulsion. 


In order to obtain an algebraic equation of degree m—1 satisfied by the #,,, 
subtract A from the expression (20) and multiply by the denominator. Then we 
obtain 


(22) - ie : 


Therefore the polynomial equation for all fixed points is given by 


On Pu APu 
(A=6;,)" (k= on)" 


m mt 


(23) in =0, 


while the polynomial equation satisfied by the % is obtained in the form 


m m 


(24) [TGs a 
w=1 6 


w=1 
The equation for the repulsive fixed points #.,, given above in normal coordi- 
nates, can be written in the invariant form 


(25) Quaz—a>(y) = 0. 


7, Since the #,, are points of repulsion, we have | p’(#,)| = 1. It can easily be 
shown that at each of the points #, we have 


(26) OW) 24. 


* A. OSTROWSKI, Mathematische Miszellen XXV, ,,Uber das Verhalten von 
Iterationsfolgen im Divergenzfall‘‘, Jahresber. d. DMV, Bd. 59 (1956), pp. 69—79. 
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Indeed, if we take A> o,, In a sufficiently small neighborhood of 6, we have 
g(a) <A. 


If 4 increases, this inequality remains valid until we get the first point for which 
p(A)=A. This is #,. We have therefore for sufficiently small positive e: 
p(3,,—€) <8,—«. On the right side, replace 0, by p(G,,), subtract p(9,,) on both 
sides and divide by —e. We obtain 


p (3, —€) —9(h,) SS 4 : 


= € 
and from this, since €|0, (26) follows immediately. 


8. The further discussion of the global convergence problem appears to present 
considerable difficulties if treated by the method of conjugate couples of points 
(cf. the paper cited in § 6). Indeed the determination of such couples depends on 
the solution of the equation 

A= 99 (A), 


which reduces to an algebraic equation of degree 4m—1; we may expect as 
many as m—1 couples of conjugate points. 

Only when m= 2 is the problem solved immediately. Indeed it follows then 
from the analysis given in Section 7 that as soon as a A, lies in one of the open 
intervals (o,, 7,), (9, 02), the sequence A, converges to o, in the first case and to 
6 in the second. On the other hand it follows from (10) that A, lies in <o,,0,>. 
Therefore the decision in this case is possible after the first iterative step, as 
soon as we have determined #,. But here the equation (25) gives immediately 
for #, the expression 

1, V3 24421 Vo + 420 Vi 
Ce a a vi ye SOF 


9. In the foregoing discussion, in order to obtain &, we used an arbitrary /rxed 
vector 7 in (5). On the other hand, in the theory of direct iteration a variant due 
to H. WIELANDT*® and called broken iteration is often used and consists in forming 
recursively the vectors &, given by 

3, == As ae , 
starting with an arbitrary vector €,—7. In combining this idea of broken 
iteration with our rule (4) we obtain the following modification of our rule: 


For any i,, define a vector &, by 
(28) 6 = (AS Eye ee, GA Oe 
where &_, is an arbitrary vector y+ 0, and then put 


(29) aye ae eS Og eae 


It will turn out in this case that the convergence is indeed considerably 
hastened, becoming cubic instead of quadratic. 


* H. WIELANDT, ,,Beitrage zur mathematischen Behandlung komplexer Eigen- 
wertprobleme, V: Bestimmung hGdherer Eigenwerte durch gebrochene Iteration", 
Bericht B 44/J/37 der aerodynamischen Versuchsanstalt Gottingen, 1944. 
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10. In order to discuss this procedure we assume again without loss of gener- 
ality that the coordinates are normal and that (6), (7), (8) hold. If we again put 


(x—1) 
E,= (x), it follows from (28) that «” = ie -, and therefore 


dy — hy 
(0) Oa Sls BATE OIL 
LM i,) 
If we put k 
Ne = LT (ty — 4.) (v= 1,--.05 B= 0,1,--), 
we obtain sl 
=a Qale) = DN, 
(31) A131 = cara 


Here we again disregard the , corresponding to the vanishing y, and denote the 
remaining distinct w, by (11). Then (31) becomes, with appropriate positive ,,: 


m Oo 

ae 
(32) pet se ae Same ag 13) A 

Pe 

w=1 Mii; 
where 

k 

(33) M,,.= TI] (¢, —4,) (Ris= 43 nte Rie= OAs 


11. Again denote by o one of the o, and by p, M, the corresponding by, Me, 
Then from (32) follows 


Se (n=O) Bu 
M? ;, 
(34) Ai -o= a 2 


2 Soe ea 
Sarge ae 
M; 2 Mik 


where in the sums >)’ the terms with the index u for which o,, =o are to be omitted. 
From (34) putting 


$i (u=0) by 
as a 
m : MP ? 
-}- eee S 
p+ 2 Pup, 


w= 


(35) dD, = 


we have again 


(36) y+, — 0 =D, My. 

Now put 

(37) @=Min|o,—o|, P=)ip,. K=2 (6,9), 
Oat jal 
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and take a 6 > 0 such that 


(38) 6 < 6’= Min (¢ x): 
By LATS 
Then we have obviously 
Ll 

(39) é= 5 <1, 

ee d*/4 

Ke ar a ee ; => 
9 

(40) Ts? 40. 


12. Suppose now that we have forx=0,1,...,k 
(41) |4,—0| <6 (Gee 10, Aitete ih) 
Then we have from (33), since o,, =o, 

(42) |M,,2|2@—6)" 


and therefore by (37) and (35), since all terms in the denominator of (35) are 
positive, 


epee Men el. 
2) Geo 


|D,| <K (a — 6)-2*-2, 
It follows then from (39) and (40) 
(43) UD OR eo ee eho (ega= Oye lar a), 
We have now from (36), since Mj? <6?**? by (41), 
(44) [Arty —9| Se?" 6 (ies Ope leeen) = 


We see that the sequence A, is convergent to o and we have for each k: 
| A, —o| <0, provided only that 


(45) |Ag—o| SO. 
13. We now prove that 7f A, ->o and if none of the A, 1s equal to o, then 
Arti 9% ; 
(46) Termes (% — oo), 


where y is a positive constant equal to one of the quotients = —. 


We assume first that (45) is satified. Observe that from (44) and (45) by 
definition of MM, we have 
k 
(47) M? ra 62 i bi (Set) see O2%+2 eek (k-1) 
co 
On the other hand, if we divide both sides of (36) by A,—o and use (47), it 
follows that 
Anti Oo 
ea 


k-1 
= |D,M,M,-1|<!D,| (TT (6 at) Oe a 


“=1 


rm | D,| sane ei 


2 
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and therefore by (43) 
(48) 


On the other hand, if we write (36) for k and k—1 and divide, we obtain 


Anti— o ee! 
(49) (A,—o)® — Dray’ 


therefore we have to discuss D, as given by (35). 


14. In the formula (33) for M,, the general factor o,,—A, can be written as 


(G79) (1 | oan) Therefore, putting 
ane k 

o—h, 
(50) Tn = (! ey 
we have 
(51) M,. = Gee Tie 
Since >) |o—A,] is convergent by (48), we see that 

x=0 
(52) Lg ae bie (k rac oS) , 
where ¢, is finite and positive, since = ES a ra 5 (4=0>4, os.) ones 
70 


other hand we have 


— dys 
Ee ae el rhe itt = O(o — A, +1) 


uk KM, 


and therefore further 


A,—0 
dG fy cece Fue) =Ol0— dns) Do uae 
But, by (48), the sum on the right is convergent and Sate ; therefore, since the 
sum on the left is ¢,—T),,, we obtain 
(53) IRs = t, + O(a — dy). 


15. From this it follows further by (51) and (37) that 


22a": & Pu o( Ane). 


Min laud) Ply 
Therefore, if we put 


, 
ae 


we obtain finally for the numerator of D, in (35) the expression 


Dt +(e). 


B=1 |o,,— a” 


On the other hand the denominator in (35) can be written as p+n,, where 


1h Pv TT (z=! = -J'>0 (eles), 
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since 6 —/A,->0 and o,,—4,->0,,—o. Therefore we now obtain 
(54) en cD) oe heey? | o(Aea =?) 
(6 + nr) Dr Os (o,—o)2* pk 


16. In the right hand sum in (54) it could happen that some values of |o,—o| 
occur twice, if there are two o, symmetric with respect to o, and it could even 
happen that two such terms od each other, if the corresponding ;, have the 


sum 0. Denote the distinct quotients ae ae which are not cancelled out by 
10 
(55) VV Va Pt yy 
Then we can write 
Anti 
(56) (o +m) D, = 8 Ol eer 


where the s, are non- vanishing constants, as long as there are any terms left, 
that is if 7 S Ais 


17. But if we had y=0, it would follow from (56) that 


— Anti—o\., 
D, =0/ Soe 
introducing this into (36) yields 
Aner — 0 =O(=* Mp), a? — O(MR), 


q2k 


and therefore by (4 
ao (5)" = 0 (e24(t-) 
6 32 
which is impossible, since 0 <<¢<1. Therefore we have 7 => 1, and it follows from 
(56) and (48) that 
(57) Deh a ee 
and (46) now follows from (49). 

Thus far we have proved (46) only under the assumption that (45) holds. 
However, if we assume more generally that 4,—>o, for a certain x, we have 
|2,,, —6| <6, and our result above applies if we put /,,,,=A,. The theorem 
stated in Section 13 is now completely proved. 

It is hardly necessary to add that our results hold also for Hermitian ma- 
trices, for which the discussion above remains valid with some slight and ob- 
vious modifications. 

Note added Octobey 1957. Professor G. ForsyTHE has directed my attention to a 
paper by S. H. CranpaLt, “Iterative procedures related to relaxation methods for 
eigenvalue problems” [Proc. Royal Soc. London, 207, 416—423 (1951) ], in which the 
iteration rules (3), (4) and (28), (29) are discussed. In particular, Professor CRANDALL 


establishes the cubic character of convergence of &, in the rule (28), (29). However he 
does not arrive at our asymptotic formula (46), which is the principal result of our 


paper. 
Mathematische Anstalt der Universitat, Basel 
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On a generalization of Airy's function 


AUREL WINTNER 


Communicated by C. TRUESDELL 


1. The standard normal form (c}., e.g., [7], p. 97) 
Day = ay 
(D =d/dx) of the differential equation 
PZ’ 4+ tZ'’+ (#—A)Z=0 
of the cylindrical functions Z = Z, (t) of index A, where 4 =y (A), can be reduced, 
it A=4, A——F of A=%2t0 
(1) DEY RY FQ; 
cf. [7], p. 189, where (1) is written in the form 
(2) Eee 10 


after a change of the unit of length on the x-axis. Since the work of WEYL [8], 
this particular differential equation has played a considerable part in the theory 
of singular Sturm-Liouville problems on a half-line 0<x<oo and, correspond- 
ingly, KRAMERS [4] observed that, to a good approximation, certain asymptotic 
problems concerning SCHRODINGER’S equation (for an “‘arbitrary’’ potential) are 
determined by the asymptotic behavior of the same problem for (1) (concerning 
(1) itself, cf. the results of ZERNIKE [10}). 


What is relevant in these applications is a particular solution, y(x) = A(x), 
of (2), this y(x) being Arry’s function 
(3) y = A(x), where A(x) =f cos(u3— xu) du; 


0 
cf. [7], pp. 96—97. Since (3) is a solution of (2), it is clear from D—=d/dx that 
(4) Jaa) 
is a solution of (1). 
2. The purpose of this note is to point out how, by an appropriate application 


of an explicit rule, the relevant particular solution (3) of (1) can be extended so 
as to satisfy 


(5) D"""y —xy=0, 
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, 


if m = 2, then (5) becomes of first order and its general solution y (x) is a constant 
multiple of e*, where X = — }2?). 


the generalization of (1) to an arbitrary m=3 (if m=3, then (5) reduces to (1); 


The explicit result in question is very old and seems quite forgotten; it 
supplies the general solution of (5). The relevant particular solution of (5), to 
be obtained for any m — 1 = 2, will lead, if m — 1=2, to an integral representation 
of A(x) which is distinct from AIRy’s integral (3). 


3. JacoBr [3] wrote a note on an observation of SCHERK [6], which, with 
notations simplified so as to adapt them to the question at hand, can be for- 
mulated as follows: For every m (>1), the general solution (containing m—1 
constants of integration c;:c,) of (5) is 


(6) Vile os i if exp(e*xt—t"/m) dt, where ))c,=0, 
k=1 k=1 
th 


if ¢ is a primitive m™ root of unity. The proof consists simply in successive 
partial integrations of the m definite integrals occurring in (6). 

Disregard the trivial case m= 2 of a first order equation (5), choose the two 
c,’s belonging to e*=exp(2z21i/m) and e*=1 as c,=1 and c,—=—1, respectively, 
and for the m—2 remaining c, set c,=0. Then the assumption c,+---+c¢,,=0 
of (6) is satisfied, and (6) reduces to 


Ble SR O2 eee —exp(xt)|exp(—?"/m)dt, where ¢=exp(2zi/m). 


But (5) is real for real x, and so both the real and imaginary parts of (7) must 
be solutions of (5). Since the real part of (7) differs from 


(8) y(x) = f [exp (xt) — exp (a,, xt) cos (b,, xt)] exp(—t"/m) dt, 
0 

where 

(9) Ay = cos(2nz/m), 6, =sin(22/m), 


only by a constant factor if x is real, it follows that (8) is a solution of (5) even 
if x is complex. 

Let x be confined to the half-line 0< x< oe, and let the integration variable ¢ 
be replaced by xt with x fixed. Then the particular solution (8) of (5) appears 
in the form 


(10)  y(x) =x f [exp (x*é) — exp (a,, x? #) cos (b,, x*£)] exp(— «”t"]/m) dt. 
0 


Since the factor [ ] in (10) is majorized by a constant multiple of exp(x*?), it 
is clear from (10) that 
(11) y(x) +0 as x—>0o 


(and that, as a matter of fact, y(x) tends to 0 very fast indeed) if m>2. If, on 
the other hand, m= 2, then (11) is clear from the explicit form of the solution 
of the differential equation (5) of order m — 1=1, this solution being exponentially 
small for large x. 
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Needless to say, standard methods dealing with the asymptotic behavior of 
such integrals as (8) could readily replace (11) by explicit asymptotic formulae. 
But only (44) (or, for that matter, only the boundedness of (8) for large x) will 
be used later on. 


4. The integral (8) is a particular solution of (5) whether m is even or odd. 
If m is restricted to be even, quite another real solution of (5) can be obtained. 
It turns out that what then results is relevant for the cases «=2, 4, 6,... 


of the function on 


(12) F, (x) = f cos(x#) exp (— 2%) dt 
0 

(cf. [9], p. 834). The family (0<a«<ce) of these even functions F,(x) was intro- 
duced by Caucuy [2] as a generalization (from «= 2 to an arbitrary, not neces- 
sarily integral, value of the positive parameter «) of LAPLACE’S symmetric normal 
distribution to an arbitrary symmetric “‘stable’’ distribution, with F(x) as the 
density of probability. C/., however, the discovery of F. BERNSTEIN [/], and the 
subsequent final results of P. Lévy (e.g., in [5], chap. V), concerning the existence 
of these “‘stable distributions”’. 


If m is even, so that (5) becomes 


(13) DPE iY = Os 
where n=3m, then both 1 and —1 are m-th roots of unity. In (6) choose 
as 1 and —, respectively, the values of the c, belonging to e*=1 and e*= —1, 


and for the m— 2 remaining c, set cg=0. Then the case m= 2n of (6) reduces to 
(14) y (x) =f cosh(x#) exp(—1#"/n) dt. 
0 


Hence (14) is a solution of (13). But it is clear that there exists a (unique) positive 
constant A= A,,, for which the case « = 2 of CAUCHY’s transcendent (12) becomes 
identical with a constant multiple of (14) if « is replaced by 74x in (12). 


Consequently, y =F, (Ag, x) ts a solution of 
(13 bis) (— 1)" Dy 4 xy =0. 


5. Consider now (1), the case m= 3 of (5), so that (14), where m=2n, is not 
applicable. 


Since (9) shows that (8) now reduces to 
(15) y (x) ul Lexp (xt) — exp (— x#/2) cos (34x #/2)] exp (— #/3) dt, 


(15) is a solution of (1). It will be shown that this solution of (1) is identical 
with a (positive) constant multiple of A1ry’s solution (4). In other words, if x 
is replaced by 3*x in the integral (15), then the resulting function of x must 
be identical with y A(x), where A(x) is Arry’s integral (3) and y a numerical 
constant (the value of which can be determined by placing x = 0 in both integrals). 

This identification could be carried out by a deformation of the integration 
path OSu<oo of (3) in the complex u-plane. But this explicit work can be 
avoided if it is ascertained that (7) the solution (4) of (1) tends to 0 as x00, 
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and that (iz) some solution y(x) of (1) must tend to co as x->oo. In fact, (ii) 
imphes that not all solutions y(x) of (1) will satisfy (11), whereas (15) is a solution 
of (1) satisfying (11), as is, by (2), the solution (4) of (1). Thus the identification 
of (4) with a constant multiple of (15) can be concluded from (7) and (iz), since 
a homogeneous linear differential equation of second order does not have more 
than two solutions which are linearly independent. 

Since the truth of (7) is exhibited by a well-known asymptotic formula of 
standard type (cf. [7], pp. 189—190 and pp. 202—203), only (27) remains to be 
ascertained. But (2) is quite on the surface; it follows from the fact that if (4) 
is written in the form 


(16) DEV (2); 
thensince)/(%)i== x, 
(17) Eoyiz208 for 0S 7% jcc 


In fact, if f(*) is any continuous function satisfying (17), and if y(x) is that 
solution of (16) which is determined by the initial conditions y(0)=0 and 
Dy(0) =1, then it is clear from (16) that, for reasons of convexity, y(«)=2x 
will hold for 0S %<oo, and so y(x)—> 00 as %¥> 00. 
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Boundedness of Solutions 
of Linear Differential Systems with Periodic C vefficients 


H.R. BAILEY & LAMBERTO CESARI 


In previous papers L. Cesart, J. K. Hare, and R. A. GamBitt [J, 3, 7, 8, 13, 
15, 17] have studied classes of systems of linear differential equations, all Batis. 
to the form 
(1) x’ =Ax+ e@(i)x 


where x =col(x,,..., X,), € is a small parameter, A =[a;,] an Xm constant 
matrix, ®(t)=[9;,(¢)] an mx matrix whose elements are periodic functions 
of ¢t of given period T =22z/w, L-integrable in [0, T]. Theorems of boundedness 
and criteria for unboundedness of the solutions of (1) have been given in [J, 3, 
7, 8, 13, 15] in case the characteristic roots of A are all +0, distinct, and purely 
imaginary. In [17] it is assumed that A has also a number of zero characteristic 
roots. In all cases it was shown that convenient conditions of symmetry, or of 
evenness and oddness of the periodic coefficients, may assure boundedness. 


In the present paper it is assumed that some of the characteristic roots of A 
are distinct and purely imaginary and one possibly zero, and the remaining 
characteristic roots are real or complex with real parts negative. The coefficients 
in ¢? of the Floquet characteristic exponents of system (1) are given in explicit 
and easily computable form. Consequently simple criteria of boundedness and 
unboundedness follow (§ 3). Numerous examples are given showing the critical 
part that such coefficients have in the present essentially more general situation. 
The same method of successive approximations (§ 2) is used in the present paper 
as in the ones quoted above, namely, the method introduced by L. CEsart in [7] 
and successively developed by L. CEsari, J. K. Hate, R. A. GAMBILL, W. R. 
FULLER for proving theorems of boundedness of linear differential systems with 
periodic coefficients (Joc. cit.) and existence theorems for cycles of weakly non- 
linear differential systems (periodic or autonomous) [2, 4, 9, 14, 16]. Nevertheless 
some of the theorems of §3 could be obtained also by a different method. 


§ 1. Preliminary remarks 


1.1. Characteristic exponents. We shall first recall some properties of linear 
differential systems with periodic coefficients, due essentially to G. FLOQUET 
[5, 6, 10|. Let us consider the system 


(A434) | Pages yltaye's. 
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where x=col(%,,..., %,), P(t) is an Xn matrix whose elements pjn(t) are 
complex-valued functions of the real variable ¢, periodic of period T=2z2/a, 
L-integrable in [0, T]. We shall consider only those solutions x(t) =[%,(#),..., x, (d)] 
of (1.1.1) whose elements x;(f) are absolutely continuous (AC) functions of ¢ 
and which satisfy (1.1.1) almost everywhere (a.e.). Then there exists a fundamental 
system of n AC solutions of (1.1.1) 


NOES AC eee enn 1a ee ac? 


which can be distributed in disjoint groups (HAMBURGER’s groups) x”) (2), 
h=k, k+1,...,kR+A—1, 15kSk+A—1Sn, each of which satisfies the 
following relations: 


CGI (A Wee ic (4 Pe eee (4) aes sa (sa 
KGET Visas 4) (g)i ey taDipy 
a) (t) = e%* [29 (2) + ee VQ) +--+ oP tet MOT A=k,....RtA—-1, 


where s, « are complex numbers (a pair for each group) and all vectors z have 
elements which are AC periodic functions of period T [z"~” (é) non-identically 
zero, h=hk, ee aid 

The numbers s above are the characteristic multipliers and the numbers « 
(defined mod wz) are the characteristic exponents. If s;, /=1,...,m, are all the 
distinct characteristic multipliers, and if « are the corresponding characteristic ex- 
ponents, if ju, is the sum of all the numbers / above relative to all the Hamburger 
groups having the same s=s,, then #,+-:-+m,—=n. The following further 
information on the characteristic multipliers is needed. Consider the fundamental 
system X™ (t) = [X,,(),7=1,..., 2], 4=1, ..., , of solutions of (1.1.1) verifying 
the initial conditions X;;,(0) =0;,=1if7=h, =Oif7--h. Then the characteristic 
roots of the matrix [X,,(7)] are the numbers s,,...,s,, with multiplicities 


j 
ie Une c sual denote Dy, .5.. >S, and joaaens te characteristic mul- 


>” 

tipliers and exponents, each repeated as many times as the multiplicity yw of s. 

From the considerations above it follows that the solutions of (1.1.1) are 
bounded in [0, + cc) if and only if all characteristic exponents « have non- 
positive real parts, and those with zero real parts correspond to HAMBURGER’S 
groups all made up of only one element (A=1). A sufficient condition (though 
not necessary) is, therefore, that all characteristic exponents « have non-positive 
real parts and those with zero real parts correspond to characteristic multipliers 
s, which are simple (u;=1). The case where P is a constant matrix in / is no 
exception. Then T is arbitrary, and we can assume that the characteristic 
exponents «%, J=1,2,...,, are the characteristic roots of the matrix P. 


Let us consider now a system 
(1.4.2) x’ = Olt, e) x, 


where x=col(%,,..., %,), Q(t, €) =[Gal(é, €)] is an mx matrix such that (a) 

all g,;,(t, €) are periodic functions of ¢ of period T, L-integrable in [0, 7} for every 

| e]| <e, and some eg >0; (0) all g;,,(¢, €) are holomorphic in ¢ for | ¢| < éy for every 1; 

(c) |gn, €)| Spd for all 7, h, t, |e|<é, where y(t) is a non-negative periodic 

function of ¢ of period 7, L-integrable in [0, 7]. Then the solutions X;,(é, ), 
17* 
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considered above, are holomorphic functions of ¢ for |e|<e, and for every ¢ 
[cf., e.g., [5], Chapters 1 and 2]. Thus it follows that the multipliers s,, Be Waa): 
(i.e., the characteristic roots of [X;,,(7, e)}), are the roots of an algebraic equation 
s"+a,(e) s”--+---+-a,(e) =0, whose coefficients are holomorphic functions of € 
for |e|<e) and a,(e) =det[X;,]--0. Hence the multipliers s, can be thought 
of as analytic functions of ¢ with at most branch points of finite order for | e|< ey 
and also we have s,-+0, /=1,...,”. Finally the characteristic exponents 
a, (e) = T+ log s,(e), L=1,...,, can be thought of as analytic functions of ¢ 
for |e|<e, with at most branch points of finite order, once the values «,(0), 
l=1,...,, have been chosen. For instance, for the system Ky = + EX, 
%=%,+%,, T>O0 arbitrary, we may assume o%(¢), %(e)=1+ //e. 

Analogously, for a system as (1.1.2), where ¢ = (&, ..., &,) is a vector parameter 
and the above conditions of holomorphism of Q hold for ¢,, ..., €, in some region 
K of the (é, ..., &,) complex space, the analogous conclusion holds concerning 
the analyticity of the characteristic exponents «,(e) for ¢ in K. 

1.2. The concept of mean value. Let C,, denote the family of all func- 
tions which are finite sums of functions of the form f(t) = e*' p(t), —°<t<+o, 
where « is any complex number and q(t) is any complex-valued function of the 
real variable ¢, periodic of a given period 7, L-integrable in [0, 7]. Obviously 
a function of the form /=e*'(t) admits the infinitely many decompositions of 


the same type which are obtained by replacing « and g by x +71Kw and we ‘*®', 
K=0, +1,.... If m(é) has the Fourier series 
g (t) oes os C ase 
then we shall denote the series " 
(221) Le" Gio) 3 Ce 
n=—OO 


as the series associated with f(t). Following L. Cesar [1], we shall denote by 
mean value, m[f] of /(¢), the number m[f] =0 if inw+a=+0 for all n, m[f]=c, 
if nw +«=O0 for some n. Both the series (1.2.1) and m[f] are independent of 
the chosen decomposition of e*’@(t). Finally let us define m[f/] in the class C,, as 
an additive functional (for more details, see L. Cesart [J] and J.K. Hate [12}). 
Obviously m[/] reduces to the usual mean value for periodic functions f(t) of 


ig 
period 7, and, in this case, we have m[f]=T J f(é) dt. 
0 
The following statements are needed in the sequel. 


(1.2.1) If f(t) EC, and f(t) =e“"'' 9(t), a,b real, g(t) periodic of period T, 
and m|f|=0, then there is a unique primitive of f(t), say F(t), which belongs to C,,, 
and such that m[F|=0. Its associated series is obtained by formal integration of 
the serves associated with f(t) (L. Cesart [1], J. K. Hate [12]). 


(1.2.11) Under the conditions of (1.2.1) we have 


vi 
FO|SNfl[g@|dt, N=N(a+i1b,T), ots, 
0 


where N is a constant depending only on a,b, T and not on g (J. K. Hate [12)). 
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The proof of (1.2.ii) given in [12] is based on the use of convolution integrals. 
The particular case a=b=0 is of interest. Here N(0,0, T)=27, as can be 
Biever oo aa as follows. First suppose f=g, real. Let us observe that 


fil at, fre t)dt=0, and that it is not restrictive to suppose that the 
0 


Ee tinmaois function F is not constant and takes, at £=0 and t=é, 0<é<T, 
respectively, its absolute minimum /, and its absolute maximum F,, R<0<fy. 
Then 


é T 
OSES Bi ae ge rg ot 


0<2(K —F) =({- {)fdts flares. 


Thus 


+ig,, F=G,+i®,, 


ais 
Pi V2, D,, ®, real, then |®,|<27/f |@,| dt, 7=1, 2, and 
0 


3 2 a 2 ss 2)2 B 9 2\1 
[Fl = (®f + OH! = 2>|(Slojai) + (Sloalad | <29J t+ abba 


ih 
Thus |F|<27/f |f|d@¢. This proves that we may assume N(0, 0, 7) =27. This 
0 


fa iso the best €valuation.- Indeed, if j= 1 for OS¢< 244, {= —1 for 2 1ast<a, 

7——m for ast<21(1+a45-f=m for 2'(1+4)St<1, where 0<a<21, 

m =a" (1—a)*, and T=1, then F=?/, a—t, —m(t—a), m(t—1), respectively, 
1 


in the ois intervals above, and |F|,,,,=2724, f |f|d#=a(1—a)+. Thus 
6 


ee ae fifi 21 (fa) > 2-8 as a>O-+. 
0 


§ 2. Expressions for the solutions and the characteristic exponents 


2.1. Notations and hypotheses. Let us consider the linear system of dif- 
ferential equations 


141) y=Ay+teD(t)y 


where y=col(4,,..., Y,), € a parameter, A = [a,,| a constant matrix, ® an nxn 
matrix whose elements q;, (¢) are complex-valued functions, periodic of period 
T =2a/m, L-integrable in [0, 7]. We shall suppose that of the m characteristic 
roots 9;,7=1,...,%, of A, k roots, say 9Q,, VAS eh OS hy rate distinct, 
have real parts R{e;|=0, with 9; 0, (modq@z), and the remaining 1—R, 
say 0;, j=k+1,...,", have R{oe;]<0 and are not necessarily distinct. If we 
denote by 0},..., o,,, the m distinct values of the » characteristic roots @,, ..., Q,, 
and f,...,,, their multiplicities, we have OSkSmSn, y=---=(y=1, 
Matas) Mm 1, Ma bn = 1. 

It is not restrictive to take A =[a;,] in canonical form, that is, A=diag(A,, 
A,,...,A,,), with A;=[o,], 7=1,...,%, A,;=[@,], s=k+1,...,m, where 
R’+1<59, ASR’ +y,, F’=wyt--stos-1, and 4;=e,= des Qji4j—=1, or =0, 
aj,=0 for j7—h=+0,1, k’+157, ASP +h. If 6;11=4+41,;, then 6;=1 or 0, 
for 7=k-+1,...,”, and 6;=0 for 7=2,..., k. 
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Let a (é),.-.,%(€) be the characteristic exponents of (2.1.1). For e=0 
system (2.1.1) has constant coefficients and, by § 1, we may suppose «;(0) =@;, 
j=1,...,n. By §4 we know that the functions «,(¢) are analytic functions of ¢€ 
with at most branch points of finite order (for ¢ complex and finite). 


If we take in the complex plane arbitrary disjoint neighborhoods U; of the m 


distinct points 0}, (aA; me then: for | e| sufficiently small, each point «; (é) 
j=1,...,k, is in the corresponding U,, and the remaining »—k points «;(é) 
(distinct or not) are distributed in groups of mw; in each U, 7=k-+1,...,™. 


Thus for | ¢| sufficiently small, say for | |<) and some ¢, > 0, we certainly have 
a;(€) = a,(€) (modo), ie Na ds aie gee 3 
Rla;(e)] <9, {= BP Ay aes 


This implies that the question of the boundedness in [0, + oe) of the solutions 
of (2.1.1) for |e|<ég is reduced to the determination of the signs of R[a;(e)], 
7=1,...,k. To achieve this we shall determine & solutions of (2.1.1) of the form 
x, =e" z(t, €), where z;(t, €) denote periodic vector functions of period T= 22/w 
in t. To obtain these k solutions we shall use the casting out method of suc- 
cessive approximations mentioned in the introduction. 


2.2. A method of successive approximations. Let us consider system 
(2.1.1) where A has the canonical form described above. Let us consider the 
auxiliary system 


(2.2.1) y= By+ eOli)y, 

where the matrix B is identical with A except that the letters 0,,...,0, are 
replaced by 7,,..., t,. In harmony with a remark made in [1] we shall actually 
assume 

(2:22) T= 0; tor all a eae bee 


where / is fixed, 1S hA<k. Since 0;==0, (mod w1),7=4, ..., k, gaeh, there aces 
circle y, of center 9, in the complex plane such that t,==e; (mod 1), 7-+-h, 
j=1,...,k, for every t€y,. Since R(o;)<0 for 7=k-+41,...,, we can take 
y, sufficiently small that R(g;)<R(t,) for all t€y,, 7 =k+1,...,. Hence 
we have 1;==t, (mod 7) for all 7=1,...,, 7h. We shall consider t, as an 
indeterminate, t,€ y,. 


By the method described below we shall obtain a solution Ve (Vy hae ee 
of a system 


(2.2.3) ie CeO). 


where C is identical with B except that 1, is replaced by t,—e d, and d, is a 
convenient holomorphic function of ¢. In order that the solution of (2.2.3) bé 
a solution of (2.1.1) we must require that T,— €4,= On, T€Y,, having already | 
assumed that t;= 9, for all 7h, 7=1,...,n. 


If we replace y,, in (2.2.1) by 


(2.2.4) 5, ny 6x eet en Pa Vy t, 
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where h is a fixed integer, 1h<k, and equate like powers of ¢ we obtain 


» ‘ 
oh alt jth 

dx 

qh , 

dt == Og gt OY; j=kh+2,...,0 

dam 

a = Th eS Pir % sg j=1,..,.k+4, mz, 
Coe, (m) : 
Se Os Dee eS hh, eet, 


v=1 


whose solutions satisfying the condition “yO )= oi, x? (0) =0, m=1, where 


as usual 6,,=1, 6;,—0 for 7+:h,7=1,...,, are given by the expressions 
ay = Oy Co. 1 = tee 
fae fem Sg, en Naa | Neil Cots circa cate 


x) — oti! fo Oj-4 of +S Gr ane fe BE Ne eel. Os ones 
0 


As is apparent, the functions successively obtained by this process are not of class 
C,, (in general), since it soon becomes necessary to integrate functions of class C,, 
whose mean values are +0 [c/. Theorem (1.2.1)]. Successively, polynomial terms 
in ¢ occur of higher and higher degree. According to the general idea of the 
method as given in [J, 4,9, 13, 14| the process above has to be modified by sub- 
tracting the mean value (at least) from each integrand and by choosing at each 
integration the unique primitive of class C, and mean value zero in terms of 
Theorem (1.2.1). By a convenient formalization the new functions x”) so obtained 
will satisfy modified differential equations and modified initial conditions, and 
the corresponding functions (2.2.4) will satisfy (2.2.3), as will be shown. 

According to one of the equivalent formalizations of the method, namely 
the one given in [9, 14], let us put 


n 

(22:5) Si,=m lem oe FEE) aoe 
j=1 

and define the successive approximations by 


0 t pee 
hp Op, Cs 7=1,...,%, 
m 


n 4 
= 1—1) 
ae a ee th é pu : yy ihe rae ae lh ah he HT at, m™ = 1 5) 


ae Ere {ac 1S Gy ES zi at, 1h, P= ANG2 Nasser, m1, 


Xin 


ff) =e fe-"| Si, V4 daa dt, j=k+2,...,m, m2, 
r=1 


jh 


where the primitives are the unique primitives of (1.2.1). We will show that the 
integrands of (2.2.6) have mean values zero. 
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We first show by induction that for m>0 we have 
a) = et! yt (), ee te ee 
where wi” ({) are periodic of period 7 and m[yph)]=0. For m=1 we have 
xh oie if Sabie: [Pan e™! — m[ Pru] oN] dt = eM iy, (t), 
where yt} is the unique primitive of mean value zero. We have also 
xp el fie! pee dt = ew! wi) (2), Co De Ff = Ai, ees ape 
x) = et? [ e— 7 (py, e™ + 07-4 etal, di =e WO), JH=R12,..-.0 


where yl, j=1,..., m, are periodic of period 7. The integrands have mean value 
zero (see § 1. 2) § since T; hil (mod wt), 7==h, 7 =1,...,”. We complete the in- 
duction by assuming x) = e™ ll) (), (Ayo ty bHa1s42, MA, wherestie 


yh} are periodic of period T and m [wh 1] =0. If we replace the xh, b=4 ).t2, MAS 
n (2.2.6) by the expressions above, we obtain 


m— 


ee Fe 
a re ha, L Pin rh ye pa “Sis Wh 1 | dt = e™! yl? (2) 
r= 


where the w\”) ({) are the unique primitives of mean value zero. The integrand 
has mean value zero since 


n 
Sma =| D gr, | and mp N= 0, 1=1,...,.m—1. 
r=1 


We have also 


ape tt f ON EST oe apa dt em ath ah (Oe ie 7A aes 


r=1 


n 
ae Ser (genre) > Dir Pn a 61 ya dt =e y(t), g=k+2,.... 
r=1 


where the yw) (¢) are periodic of periodic T. The integrands have mean value zero 


since 1;==t, (mod 7), 7==h, 7 =1,...,”. The induction is complete, and we can 
write 


(22:7) Vine xin = MR Sas fee Tee ear 
m=0 m=0 
Remark. The formalization of the method described here follows closely the 
one given in [9] and [14] for nonlinear systems. The proof of its convergence 
given below is a somewhat shorter form of the one given in [9] and [14]. Details 
are pointed out which are needed in §3. For a formalization valid for linear 
and nonlinear problems under the Lipschitz condition alone, see [2]. 


2.3. Proof of convergence. We now show that series (2.2. 7) is absolutely 


and uniformly convergent for | «| sufficiently small and for all ¢ by determining 
majorants &,, such that 


(2.3.4) lw? OLS Ee gi Ageia = Oe 
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for all ¢, and such that > &,,¢” is absolutely convergent for | «| sufficiently 
m=0 

small. We can choose &=1 since yy) = 1 and wy) = 0, 74-4. Let us assume 

that we have the first m—1 majorants, that is, pol <&, fe Naas eb 

2,...,m—1. Then we have 


|S; | =| me pn xy 4 


=|". f Doli pdt 


eS Taf 3 lonllye |dtSn&,MT4, 1=1,...,m, 
0 j=1 
where 
T 
Vie SMa TO OGl we. at, Vice Anant 
0 
and finally, by (1.2), 
|p? ()| = rip x Par Pon P&S phi »| at| 
= Spee suey (n &1M Lies) Sey he 0 = t = fh 
1=1 
We have also, by (1.2), 


wn (| = 


j= t = m—1) 
eT Ta) eae *) Soe oh at 


< elk (tj—m) T| Nn ME, 1=K,,nME 


m—1? 


for all OS/ST, 7=1,...,k +1, 7==h, where N;, are the constants N(t,—1;, T) 
of (1.2.ii), and K;,=N,,e!"%—™l, Successively we have 


n 

ple, n(d)] = [eter ! fem tete | Ys ln + Oa pe a 
r=1 

Cea) [Ni +2,h (n M sina at Ona Tih In M C= a) 


=nMé 


IIA 


m—1 


Feats Ose 
where N,,5, is the constant N(t,—t4+2, 7) of (1.2.11), and 


Kee Neen (tr Oar eee 
Analogously we have 


| wir? (t) Sau eee ,K. jh» 7=Rk+3,....0 


for convenient constants K,, and all 0S/ST. If we choose P=max(nMK;,, 
nM, 1) for all7=1,...,”, then we have 


| pyr (t) )| = Pie + 2h a 


for all OST and 7=1, 2,..., n.-Thus we can choose 


m 


(2;3;5) So = PlE,-1 + 28-1 é,.-1) , 
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and we have |p" (é)|<é,, for all OStST. Relation (2.3.3) defines €,, for all 
Meats 2s 52, oS thus the inequalities (2.3.1) are proved for all OS#<T. Since 
the Pen fone yn (t) are all periodic of period T the same inequalities are proved 
for all ¢. 

It remains to show that Se, |e|” converges for |e| sufficiently small. Let 
us consider the function ™”~° 


(2.3.4) BG, eran ePe (eee. 


e,€ real, and observe that F(&,0)=0, F(&,0)=—1. Thus by the implicit 
function theorem, we conclude that there exists a real function = &(¢) satisfying 
F(é, s) =0, &)=&(0), for all |e| ey and some ¢,)>0 sufficiently small. We have 
also, for ¢) sufficiently small and all | ¢|<ép, 


(2.3.5) ff (era, o 


where this series is absolutely convergent for all |«|<e). If we replace &(e) in 
(2.3.4) by its expansion (2.3.5) and equate like powers of ¢, we obtain 


m 
49 = bo, a, = Pet = is ea ap—1 Ay —k» 
oa. 


and this is the recurrence relation (2.3.3) for the &,,. Thus by taking &;=1 we 
have &,—a,,, and the series 


m m? 


ie Ne: 88 Pe 


m=0 


converges absolutely for |¢|<¢g. As a consequence, each series 


Pin (lt Soe yw? (b) ) cee A IO 
m=0 
converges absolutely and uniformly for all t, — ~o<t<+ ~, 


We now show that the functions y,,(/) given by (2.2.7) satisfy system (2.2.3). 
Let us put 


(2.3.6) by see Sens 


m=1 


This series is absolutely convergent for | ¢| sufficiently small since, by (2.3.2), 
it has the following majorant 


, 


US. n M yy Ep, Chae 
m=1 


and the latter has been shown to be convergent. 


By (2.2.6) we have for 7h, feet i eel 


da) 


(2.3.7) Me pee Gh 1 as > e™ may my Pir gatas 
m=1 r=1 


m=0 m=1 
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where the series are all absolutely convergent for |e|<e 9; thus the first series 
(2.3.7) is the derivative of y,, and hence we have for all ¢ 


dy, a 
Te =, ve oS Pir Vr 
r=1 
Analogous computations hold for 7=k+2,...,. For 7=h, we have 
ax ™m 


” co co n [o-) 

Ch) we m (Mm) m (m—1) m m—l 

dt = >'é Xhh Ne > Par Xrh = ee > Sm e 
m=0 m=1 r=1 m=1 


1=1 


(2.3.8) Oe 


m=0 


The first two series in the right hand member are absolutely convergent. The 
third series is also absolutely convergent as the product of the two absolutely 


convergent series 
co 


oe She h» 2 é x7 9 
m=1 m=1 

Thus the series in the left hand member of (2.3.8) is the derivative of y,, and 
we have for all ¢ 


Naps ley é dy, Sie Se Pir Vr- 

T—1 
Thereby we have shown that the functions (2.2.7) satisfy (2.2.3). The same 
functions will satisfy (2.2.1) also if, having chosen t;=@,, 7h, 7=1,...,”, we 
can determine tT, =T,(€) in such a way that 


(2.3.9) F(t) = Th Fas ©) — 0, = 0: 


Now F is an analytic function of t, and « and Fi(o,,0)=0, F,(o,, 0) =1. 
Thus, by the implicit function theorem, there exists for all |e|<¢) and e9>0 
sufficiently small, a function t,(¢), |e|<é, with t,(0)=@, satisfying (2.3.9) 
and given by a series 


lo) 

a 4 

Th ae On sin oD: Any é 
ae 


which is absolutely convergent for | ¢|< ep. 


Remark. Each t;(e) is an analytic function of ¢ for all | e| sufficiently small, 
and thus the same holds for the numbers S,,, and for the Fourier coefficients 
of the functions yi, 4=1,...,%, m=0,1,..., aS can be proved by induction. 
Wiso we have 4, =-70,,.6(S,,)=—0 for e=0,4—1,.. 2, =A, a0. 


2.4. Expressions for the solutions and the characteristic exponents. We 
may now conclude the considerations of the previous sections (2.1), (2.2), (23.). 
Under the hypotheses (2.1) and for all |¢|<¢) and some ¢)>0 sufficiently small, 
the system (2.1.1) has »— characteristic exponents 


(2.4.1.a) a, (6) =% —0,+0(1) as e>0, Wi ho 3 0; 
and & characteristic exponents 


CAD) (2) = tp = Oo + 6, (6), In == Xo connlts 
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a 


where 


(2.4.2) dy, ar 2 ai Car 


n Tr: n 
2.43) Sam le-™ DS par (t) ¢-®| = TA f e-™D gag lt) 4-2 (0) de, 
( 3) mh a? ; ee 


and 
xy ere, ay 0, (ce eee Pee 


xy) =a Ca eae [> x Pray % a Se a Sin » dt, 
(2.4.4) 
apes CO O "S giet xi) dt, a ee ee See 


Dy at | a w!|> 9, oe V4 6 1%)— |e, aM ie ee 
The corresponding solutions y, =(yy;,---, Yn) Of (2.1.1) are given by 
Vin = a Xp, jai Wee Rae | 


m=0 


where h is any number h=1,..., k. All integrals above are the unique primitives 
of mean value zero considered in (1.2.1), since the integrands also have mean 
value zero. 


If we assume 


0d8 


Pjnlt =e (a;n,c0slwt + b;,,;sinl@t), 1s Boe aoe ee 


1=0 
(4;41, jn. Teal or complex, 0;,9 = 0), 


then 
(2-45) Si,=m lem y Pri a) = (Pen) = Gano» fk be Biel | 
it 
os e™h! pity sah [Pan (t) — m0 (Pyn)] at 


oe) 
ca tae ah (Lm) + (= Orn  Coslat + Ay ny Sin lw t), h = ‘ls niches k, 
i 


ANG SIT Oe == (), 7 2, wa, Me ALSO 
Hh em yl = ott f oe, () dd = oD [(ty— 0)? + Vo}*] X 
1=0 
x fe 0;) 4n1—L@ bj41] coslw@t + [(t, —@ j) Oni L@ a; n7] sin lq th, 


(=1, 50) Fate WA ks 


oy = m| 2% ys) 
(2.4.6) oe Bets = 
Sr 20 “hI 0 Be 
2 Re + ter Cal) Clo) Catto Dayal, 


t=1 
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where h=1,...,h, Cijgi= An 51%; nit On 51 OF nL, LEGS EES ps For e=0 
the same ooh holds with t= ,. If we take e,— T= On; + 2B; On, Ba; 
real, then, for h=1,..., k; we-have 


n 
4 24475 Up ] 
Spite - 7 70 Aino ! 
ale 2} pa OF == B2. 
j+h 
pee oj Bayt? @*) Cy i, +l (a4; By; +? ©) Dy; “}h ; 
(2.4.7) t=1 (ah 5 — Bi, +P)? +403 BR, 
z i 2By ji Anjo 4 jho 
eet Oni + Pay 
7h 
Se ee Aj Es 
= (oi; Bi, 12 @?)? 4a Bi; 
We have also, for the characteristic exponents «,(e), h=1, 2,..., R, 


oy (e) ea On ar Em [Pal ae & (Syl) ais O (e*) : 


§ 3. Boundedness theorems 
We shall prove first two theorems concerning complex systems of first order 
differential equations. We shall then apply these results to real systems of first 
and second order differential equations, since usual transformations reduce these 
to the systems above. 
3.1. Canonical systems of first order linear differential equations 
(3.1.1) Consider the system of first order linear differential equations 


dy 


(3.44) ae =Ayte@y, 
Where oy = Col (Waine- 15. 4),,), A S= diag(t.04)). 9 1 0g)) O15 uy 01, real a, ==, (mod a), 
J+=h,j,h=1,...,n, where Bis annxn matrix whose elements ~;;(t) are complex- 


valued functions, periodic of period T =2n/w, L-integrable in [0, T|, with 
foe) 
Pin (t) ~t Di aj,,Coskot, Anz, veal, 
k=0 


then for ¢ real and sufficiently small in absolute value all solutions of (3.1.1) are 
bounded in (— oo, +00). 


Proof. We shall prove (3.1.1) by showing that the characteristic exponents 
are purely imaginary. Using the algorithm of § 2 with k=n we have, for each 
[P= Moodie 


(3.1.2) Tye top e Dey | Spas 
k=1 
where, by (2.2.5), S;,=™ [Qn] =*4%:n0- By the remark at the end of § 2.3 all 


MAES eet at ie ey Il IN AID ately ATES TOE | e| sufficiently small, analytic func- 
tions of ¢ with at most branch points of finite order. Also, R(t,), R(S,,,) are 
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zero for e=0. Thus each of the functions R(t), R(S,x) is either identically 
zero for all | e| sufficiently small, or O(e”) for some maximal integer y (not neces- 
sarily the same for all these functions). We shall prove that the first alternative 
holds for all t,, S,,,- Indeed, this will be an immediate consequence of the 
statement above and of the following one, which we shall prove. For every integer 
SEG Se ots 


(ay) R(t) PRS) =O? forsall pe A fla, EAN ye VS 


(as) a ee Ne (al, cos kot +4 bin, sin kw t) + un bl, 
k=I 

where all coefficrents ast}. , oi”, are veal, and Up t) ts a periodic function of period 

T, continuous in (— co, + 00) with uy (t)=O(e’) uniformly in t, 7,h=1,...,M, 

[eS nian i= Olean 

We shall prove both (a,), (a2) by induction. 

Obviously R(S,,)=0, t=to,+y,=10,+B,, on, B, real, y,=O(e); hence 
R(t,) =f, =O(e), h=1,..., n, and thus (a,) is proved for y=1. Also, by (2.2.6) 
we have 

xp) a= ene [eee Pah em Se (Pun) ] at 
(3.1.3) oo o0 
=ie™> (kw) 1a,,,sinkot=ie* > bY) ,sinkot, 
k=1 k=1 


where all b") are real. Analogously we have for 7h, 7=1,..., 1, 
(1) it ; S 

ay = et! f Kee D, (t) dt = e&' f cet Sia, cosk@idt 
k=0 


(aca) PAN Ds [(t—@;)?+ (@)?] + [(tTr—0,) 4jnx COSR@t*+ Roa;,,sinkot] 


= etht 


iMe 


; ; 
(ayte cos kat + 1by7, sin kot) + in (f)|. 


where all real coefficients a\}),, 0"), have been obtained by replacing t, =io, + Br 
by zo;, and where, therefore, 


[oe 
os} (t) a [MjnpCoSkat+ VingsSM Rot], 
=0 


B-1.5) rav=i{ (tr —2,) [lt —@,)?-+(Ro)*]4 —4 (0,0) [— (04-0) ?°-+(2O)"] app, 
Ving =t{ko[(t —0))?+ (ko)?]2 kol— (oy Oj (ko)? |} ayy p- 


We shall now evaluate uy (t). To do this let us observe that, as in § 2, we have 
supposed | e|Se for SOME é > 0 sufficiently small in such a way that t,— 0;+0 
(mod w7) for all jh, j,h=1,...,n. We may now determine N=1 so chee 
[(t% — 9)? + (ko?PJASNA+k)%, |[—o,— o;)? + (ko)?] >| S N(4 + k)? 
|(ko)* — i(t, — 9) (6, —o)| SNU+A)2, |o(27,— 20 en ae 


la;nx| = NY 
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for all7==h,7,h=1, ...,n, and all | é|Sé,. By trivial computation we now have 


| ne | =F | an! % |[(Ae)? 4 Tage 0;) (0, — 0;) | (te > 103) | x 
x [(t, — @,)? + (Ro)?] +» [— (6, — o;)? + (Rw)?} 1S N5 (1 + k)? |B, |, 


and 


|>srel == | ane |B, kw (2t — 20; —Br) Mie Oe + (Ra@)?]+ [— (o, —9;)? + (ko)? | | 
<N*k|B,| (4+ k)-4 


Thus the series (3.1.5) is absolutely convergent, is the Fourier series of yu") ({), and 


[| <2NF|Bsl D+ 4-2 = OC) 


uniformly in ¢. Thus both statements (a,), (a,) have been proved for »=1. 


Let us suppose that both (a,) and (a,) hold for all integers 1,...,»—1, and 
let us prove them for ». 


By (3.1.2) at once we have R(t,) =O(e’). This allows a better evalutation 
for x{), ay. Indeed by repeating word by word the considerations above where 
B;,=O(c) is replaced by 8, =O(e") we obtain uw (!) =O(e’) for all 7,h=1,...,” 
Thus (a;), (ag) are proved for #7=1 and ». We have now, for m=2, 


n 
Sh = = it le tt8 Py Pn hj ie »| 
6) 
_ m| >(s 2 1 Ay, ;,COS Rw i F ¥ (a (ann? cosk@mt+10%," sin kot) + wp? (}f 
q=1 


for h=1,...,, where uw; ~” () =O (e’) for m= 2, and the mean value is the usual 


one for ae functions of period 7. Thus this mean value is 77+ times the 
definite integral in [0, 7] of the expression in braces. The sum of products of 
the trigonometrical series is the sum of a purely imaginary even function and of 
a real odd function in ¢. Hence the real part of its mean value, S,,,, is zero. 
The part of the expression in braces depending on the functions mw is the sum 
of products of functions w which are all O(e’) and of L-integrable functions. 
Hence its mean value S,,,—S,,, is O(e’). Thus we have R(S,,,)=R[Sj.2+ 
(Sini— Sma) ] =O +0 (e’) =O(e") for all h=1,...,n, and m=2. Also, we have 
for m=2, 


t Saat (m—1) (m—1) 
em fle 4 Dosa th —>s ih *hh at 


Do =f 


-- > Sul 3 oy (al,? cos kot+iby,” sin kot) + ue? |} at, 


n co 


ap > Fay COS hoot] . da (a? coskot+i bg, Ysinkot) +h Y(t) |— 


j=1\k=0 = 


where, for m=2, all functions are O(«”) uniformly in ¢. We may decompose 
the last integrand into a sum A(é) + B(t), where A(é) is the sum of the products 
of the trigonometrical series in the first line plus the sum of the products of 
the purely imaginary numbers S,,,, above times the trigonometrical series of the 
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second line. Thus A(?) is what the last integrand becomes if 77, is replaced by 
its imaginary part io,. Since in § 2 we have proved that the integrand A+B 
above has mean value zero for any value of the indeterminate t, in a neighborhood 
of a we conclude that m[A+B]=0, m[A]=0, and hence m[B]=0. Now 
A(t) =A,(t) + A, (t) is the sum of a purely imaginary even function A, (¢) and 
of a real odd function, hence m[A,]=0, m[A,]=m[A]=0, and A(é) admits 
of a unique periodic primitive of mean value zero which is the sum of a purely 
imaginary odd function and of a real even function. Since m[B]=0, also B (t) 
admits of a unique periodic primitive of mean value zero which can be evaluated 
by means of (1.2.ii). Actually B(#) is the sum of products of periodic functions 
or of constants which are all O(e”) times L-integrable or continuous functions. 
Thus by (1.2.1i) that primitive of B(t) is O(é’) uniformly in ¢. In conclusion we 
have 
ayn) — ett 2 (al), coskwt +70), sin kot] + ul? (0, 


with al”),, d””), real and wl) (t) =O (e’) uniformly in ¢ for m=2. Finally we have 


n 
ray _ ee? f Birt 2 oni) at 
n 


(3.4.8) = @8i! f elt—eilt D [ Dajrncos hoot) x 


J=1 


x »S (avert leoskat +7102 sin kot) ult 1 (t)| dt, 
k=0 


where, for m=2, all functions w are O(e") uniformly in ¢. We decompose the 
integrand into a sum A(éf)+ B(t), where A(t) is the sum of the product of the 
trigonometrical series times the exponential function. Then m[A(t)]=0 since 
T— Q;=0 (mod) for all 7+. Also, A(t)=exp(t,—@,)t-A,(é), where the 
periodic function A, (f) can be thought of as the sum of a purely imaginary even 
function plus a real odd function, and thus the unique primitive of A(¢) of mean 
value zero has the form 


f eltr—en)? ps (1%jn,COSkawt + BinxSin kot) dt, 
k 


where the trigonometric series is the Fourier series of an L-integrable function 
in [0, 7]. By (4.2) and formal integration we have a series of the form 


Pal (t% — @))? + (ko)?]*[ (t, — @;) 10j4, 0S kot + Roia;,,sinkot— 
—kaBj,,coskat + (t, — @,) Binysinkot]. 
We can now repeat the same considerations as for (3.1. 4) by performing on the 
terms in the brackets the same manipulations as on the terms in the brackets 
of the second line of (3.1.4). The conclusion is that, for m—2, we have 


(3.1.9) she my [afi cos kot + 1d! sin kot] + us (), 


with alt, om, real, and un (t) =O(e’) for all 7h, j,h=1,...,n, and m=2. | 
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The whole process can now be repeated successively for m=3, 4,...,”. Thus 
(a), (ag) are proved for every v. As noticed this implies R (t;,), R(S,,;) identically 
zero for all | | sufficiently small, and (3.1.1) is thereby proved. 

(3.1.11) Consider the system of first order linear differential equations 
(3.1.10) oY =Aytedy, 
where y=Col (yy, .--, Vn), A =diag (a, --., Qn), R(o;) =0, 7=1,...,%, R(o,) <0, 
g=k-+1,...,m, for some 0<kSn, 0; 0, (mod) for all 7+h, ee eae 
e>0 a parameter, and ® an nxn matrix whose elements are periodic functions 
of period T=2n/w, L-integrable in [0,T]. The following conclusions hold: 

(1) 1f0SkSn, and R[m(Mp,)|<0 for allh=1, ..., k, then for e>0 sufficiently 
small all solutions of (3.1.10) approach zero as t>-+ ©. 

(2) If O0<kSn, and R[m(q,;,)]>0 for at least one h=1,..., k, then for e>0 
sufficiently small infinitely many solutions of (3.1.10) are unbounded in [0, +00). 


Proof. Systems of the type (3.1.10) are discussed in § 2 and the characteristic 
exponents are given by 


on = 0, + O(e), h=k+41,...,n, 
Lp, = T = On + €M(Pup) +O (E?), h=1,...,R. 


Since R(o,)<0, we have R(a%,)<0O for all h=k+1,...,m, and «e sufficiently 
small. For h=1,...,k, we have 


R (a) = eR[m(ar)] + O(e?). 
Hence in case (1) we have R(«,)<0 for all h=4,...,; in case (2) we have 
R(«,) >0 for at least one h=1,...,k, and all ¢>0 sufficiently small. Thereby 
(3.4.11) is proved. 
(3.1.111) Consider the system of first order linear differential equations 


(3.4.14) —— et Ay oD y 


where Y= COL, «1-9 Vay A = diag (O,,..~, @,), K(o) =0, 7=1,..., Rk, R(o) <0, 
j=k+1,...,", for some 0Sk<n, 0; 0, (mod wz) for all 7h, 7,h=1,..., k, 
é is a veal parameter and ® 1s an nXn matrix whose elements y;,(t) are periodic 
functions of period T =22/w, L-integrable in [0,T |. The following conclusions hold: 


(1) If either (a) @jr= Mj Vin» Vin= —Paj> ML Pan] =O for all 7 eh, 7, h=1,...,n, 
or (b) Qjp=t*Lj Vins» Vjr=Vrjp, ™ [7%] =0 for-all 4==h, 7,h=1,...,, where all 
functions w are real, all numbers 1; ave veal, ==0, and of the same sign, if OSR<n, 
and for every h=1,...,k, at least one of the functions p,;,7=R+1,...,, os 
not identically zero, then for |e|=+:0 sufficiently small all solutions of (3.1.11) 
approach zero as t—>-+ ©. 

(2) Tf either (a) Pin 4; Wih> M[Prn] = OF Vi n= WYnj> m [Pra =0 for all 1==h, 
j,h=1,...,%, or (Bb) Pjn=tH; Vin» Pir= — Yn; for all 7==h, 1, WA, 1, Whee 
all functions py are real, all numbers qu, are veal, +-0 and of the same sign, if 
O<k<n, and for at least one h=1,...,k, at least one of the functions q,;, 
j=k+1,...,n, is not identically zero, then for | é| + 0 sufficiently small infinitely 
many solutions of (3.4.11) are unbounded in [0, + 0). 
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Proof. Systems of the type (3.1.11) are discussed in § 2 and the characteristic 
exponents are given by 


Oy, = 0,104) as e—> 0, 7 ee es 
Ot, = On + EM (Pan) + Sono + Ole), W=1,--.,%- 


Since R(o,)<0, we have R(%,)<0 for all h=k+1,..., and | e| sufficiently 
small. For h=1,...,%, we have, in all cases (1a), (1b), (2a), (2b): 


R (a) = &® R (Seale—o) + O(8); 
where R(S3;|--9) can be obtained by (2.4.6). By hypothesis we have 
Wih ~»> (4;,,cos lot “+b b;,1,Sin lot), 
1=0 


where 4,7), 0;,; are real, (b;,9=0). In case (1a), from (2.4.6) we deduce 


jh 


riacté feria jo . wy Ij (hg TP +22 @? ) Hj Ha (Gir + On gd) 
(3.1.12) R (Soule 0) ia » e ! Y +2 4 Be 
2 oon Bi at (%j5— BR wm)? ar Oj hj 


pen 
where 0;—0,= enue ,;, By; real, h—=—1,...,k. Thus «,;=0 for7=1,..., 2, 
0%, ;<0 for g=k+4,...,m; and R(Sal ah 20: far each! P=1~.. ks Since mion 
CVETY =I he, i, least one of the functions g,;(#), 7=k+41,...,#, is not 
identically zero, we conclude that R(S3,|--9) <0 for every h=1,...,k. Thereby 


(1a) is proved. The same formula above and the same argument hold in case 
(1b). 

In cases (2a), (2b) the formula above holds where } is replaced by —4. 
Hence (So s|sceye0 forat least onesiif uncer: 


3.2. Systems of first and second order linear differential equations 
(3.2.1) Consider the system 


(3.2.4) y"+Ay+eP(t)y=0, 

where y=col(y,,...,¥,), A =diag (ot, ...,0%), D=(f;,(0), 7, A=1, ..., m, where 
all f;,(t), € ave veal and o;>0, where He functions f;,(t) are periodic of period 
T=2a2/w, L-integrable in fo, T], and all 20,;, 0;+-6,, 6; —o,#0 (mod), 7h, 
p,h=A1,...,. If all functions f;,(t) are even, than for ay |e| sufficiently small 


all solutions of (3.2.1) are bounded in (— c, + ov). 


Proof. As usual the boundedness of a solution y in (3.2.i) refers to the 


boundedness of |y,| +|yi|-+---+|y,|. Statement (3.2.1) is a corollary of (3.1.1). 
Indeed by the transformations y,= w,,_1, yj; =W.;,7=1,...,” We ;—-1= 29;-1+ 29;; 


We; = 10; %;-1—10; 2g;, J =1,...,%, System (3.2.1) becomes 
/ . 
295-1 = 10; 2g;_1 — €(270,) >> hall) ) Zonta Zea), 
/ . 
2a; = — 10; Peet (270;) 2D fal) (294-1 + 4en), Moca siat Wash ars 28 


and (3.2.1) follows from (3.1.1). 
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Remark. Statement (3.2.i) is a particular case of the following, proved by the 
same method and a different argument by L. CEsart, J. K. Hare, R. A. GAm- 
Biri: 7,10 |: 

(x) Under the same general hypotheses as in (3.2.i), if @=@,+ @,, where 
P, = diag (4, ..., Y%,), the non-zero elements of ®, are all below, or all above 
the matrices 4%,...,% in ®,, and for each Y= (y')) either wy) (—2) =y) (0) 
for all u,v, or p(t) =p) for all u,v, then for all |e| sufficiently small all 
solutions of (3.2.1) are bounded. 

For other statements including (3.2.1) [but not (*)] see R. A. GAMBILL [7] 
and J. K. HAE [13, 15]. 

(3.2.11) Consider the system 


Vy ae ed hill Ve 0, 
n=1 


(3.2.2) : 
y+Gy+ay+eLtrOw=0, F=2,...,0, 
h=1 
where c;>0, d;>0,7=2,...,n, € veal, n=, f;,(t) ts a real periodic function of 


period T =22/w, L-integrable in [0, T| with Fourier series 
[oe 
hin (2) ~»> (411 COS lat -|- b;,; $10 lot), 1, h == Al. eeey nN. 
i=0 


If either ayo 0, € %y9>0, OF Ayp=0, M <0, ¢ +0, then for | e| sufficiently small, 
all solutions of (3.2.2) approach zero as t>-+- ~. 

If either a= 0, € Ayy9<0, OY Ayo=0, M>0, c= 0, then for | e| sufficiently 
small, infinitely many solutions of (3.2.2) ave unbounded in [0, + ov). 

In both cases M is given by 


M= Ba D Pp Cyt GiPyp 


j=2 1=0 

Pio= d;*, Gjo= 9, Mi = [d? as 2 wy? (CF SR) ON, 
(3.2.3) pj = 24 m1 (d; — Pw), qji=21m,)1¢;, 

Cyy1 = 4451 411 or Dy 51 Dirt, dD, j1 = 41451 7a Dit M51 


b:n0 = 9, Ge = AH BO AS aoe 


Proof. As usual, boundedness, unboundedness, and approach to zero refer to 
Iya] +]y2]+]y2}-+---+]¥,]-+|y,|- Suppose first cf+-4d, 7=2,...,”. Denote 
Diyos ea) a8, the roots..o1 the equation=7.--¢-74-d,—0, 7==2,...,%. By 
the transformations y,=W,, Y;=W ;~2, Neo eee Dasa 1, A. Wry, 
Wy j;—2=%2j;-2t 29j-1, W2j-1=72j—-2 22j)-2 1 "2j-1 %2;-1> 7=2,...,m, system (3.2.2) 
becomes the following system of 2m —1 first order equations: 


n 
4 = — ef %4— ary: (255571 ena) 
A—2 


Ww 
fi1% + fi (oye Ae , 
b= 


/ 


(3.2.4) %a;-2 = Y2j-2%2j-2 1 Oe Ue Yilg) os 


n 7 
Le = oye Sey eal oe od Cy er es Yo;-9) * Me oleae Di fin (Zen-9-+ ‘a1 , 


h=2 
i = Poop RLCe 
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If we denote this system by 
z=Az+eP(i)z, 


where z= col(z,..-, 2e,-1), 4 =diag (0, 72, %» ---» Yen—2> Yan—1)> P(t) = (Ga), 
fie Nan 2 1, With 


oo . , , 
Pin ~ d. (4x1 CoS lot + b;,,Sin lot), 1 Mal, ent 4, 
i=0 


then we have, for instance, for all s=2,...,m, 4=0,1,..., 


Ay =— A, b= — byy1, 

Ay 95—2,1= — & 515 eer see 1 

(3.2.5) yp a49 = — & 51, bi deat ALAS? 
Ce ee es — 4, 1141= Veea4 — Tasca) Ogah; 
ee = paar: =e 435-2) Veal 


If a (e),...,%,—1(€) are the 27—1 characteristic exponents of system (3.2.4), 
then we have 


ay (€) = — em(fy,) + O(E?) 
aj(e) = 7, O(1) <asy 2-0, fa 2,90 2 4s. 
Thus R[a,;(e)]<0, 7=2,...,2%—14 and o(¢) real for all |e| sufficiently small, 


and «()<0, or >0 according as m(f,1) 0, €m(f,41) >0 or <0. If we suppose 
m(f,1) =0, then 


0% (6) = & (Seile~o) + O(e"). 


We have now, for h=41, 0=o,=7,9 =2, .., 2— 1, and, by 2.4.0) 


2n-1 oo 
M = Sa,|.20 = 27 py 27 i jo4jro+d (? + Po) (—7,€,;:+loD,;)|- 
By (3.2.5) we have 
Creear= = ner “atau aa Cas 
Dione = Dy oun Td cent Pacesl ates S = 2,...,%, 


and, then, by trivial computations, and by remarking that 7,,_.7:,-1=4,, 
Yas-2+1es-1= —¢C,, we obtain (3.2.3). If, for some 7 =2,..., we have c =44d,, 
formula (3.2.3) still holds, since we can regard S,,|,9 as a continuous function 


of c;,d;, and thus we can obtain the same formula by a passage to the limit. 
t A 
Examples. The equation y’+ ef (t)y=0 has solutions y = Cexp(— ef f(s) ds} : 
0 


and thus they all approach zero as t->-+ co if m[f]+-0, em[f]>0, or all but 
the trivial solution are unbounded in [0, + 00) if m[f]-:0, em[f]<0. 


Consider now the system 


yi + (ecost-y, + esint-y,)=0, 
0, 


Yo +3 y2 + 2y¥o+ (— esint-y,+ ecost-y,) = 
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Here’ 7,;= cost, f/,=sint, f.4=—sint, /oo= cost, d,=2, .¢.=3,. and hence 
Ci21= —1, Dyy3=0, and M=—207<0. All solutions of this system approach 
zero as t>-+ oc. 


Consider now the system 
y, + (ecost-y,+ esint-y,) =0, 


Vo +343 + 2¥_+ (ecost-,-+ ecost- ys) =0 


Wwithiey 1 —C0S7,./,,—sint, -f>,— cost, j,,—cost, and hence C,.,=0, D,,;—4, and 
M=3.-201>0. Thus infinitely many solutions of the last system are unbounded 
in [0, + ov). 

(3.2.11) Consider the system 


wtau+eX hry =, 
(3.2.6) ve 


n 
Vp C5 Vj + ayy; bed finl Vireo; fame es 3 Mg 
n=1 


where a>0, ¢;>0, d;>0, 7=2,...,, € real, n=1, f;,() real periodic functions 
of period T=22/w, L-integrable in [0, T] with Fourier series 


fin (t) ~D (ax1,008l@t + 6,,;sinlot), Ra peasy ne 
l=0 


If either ayy 94-0, €4,49>0 OY a49=0, M<O, € 0, then for | e| sufficiently small, 
all solutions of (3.2.6) approach zero as t—>-+ ov. 

If either a, 9-0, €449<0 O07 449=0, M>O0, e+, then for | e| sufficiently 
small, infinitely many solutions of (3.2.0) are unbounded in [0, + cc). In both 
cases M 1s given by 


co 


M=)) Dd (bp Git GPa), 


Sis 
Pig — 2, Pur= — (2a)7 (a + Po’), 
ho=Hni= 09, pe =a d;*, Gjo—9, 


(3.2.7)  my= [a7 + 2 (c? — 24,) + bot}, Gee as Mh 


Ss 
i) 


b= — (2a) m,;,(2o% — d;), d= (2a) m, 1c; lo, 


(hee Ui a Oo 


oi n0= 0, Qii= Ay 1 Gait by 512 ;11, Dy j,;= Gyr Dy 51 = Dist 41, 


jh=1,....m, 1=4,.... 


If n=1, and a4)=0, then we have necessarily M <0, and the first alternative 
holds. 

Proof. As usual boundedness, unboundedness, approach to zero refer to 
lyn} toi] +--+] yn] +o]. Suppose of ++ 4d;,7=2,...,. Denote by 79;-1, 72; 
the roots of the equation 7?+c,7+d;=0, 7=2,...,”. By the transformations 
Vj = Wej-1) Vy, = Wey, qe tate, and #,=2,123, We=@2o, We;—1= 223-1 + %2;, 
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Wy; =Voj-1 22j)-1+ 72; 27, System (3.2.6) becomes 


n 
/ —_ > 
Bit = SLAW iis (2o4-1+ Zen) » 
15k 
n 
/ - y 
ty = — 424+ 67 D) fy, (Zon-1t 220); 
ea} 
(3.2.8) ‘ " 
| 
Zap —1 = Vaj—1 22j—-1 € (taj Pega) D fin Zen—1 + 224) » 
h=1 
n 
; 
295 = 195205 + E (2; — Yoj—1) )* Lafial 24-1 1 Zen) 7 =2,...,%, 
h=1 
or 2=Az+e@z, where aie Ze,), A=diag (0, —4, re, ---,%2,), P= 
CO 
(G0), 1 4= ,2n, with g;,(E a2 (4;,, cos lwt+6;,,;sin lot), and, for 
instance, bias 
Qy9s—11= Gast — U's r= Oe ee ee 
As 95-1,1= 42,951= 4 Ys, Sa 12,225, 0, 
Cig ef (Coa Tocca) Orie S=2,...,%, 
Ay 51,1 (Yes — as—1) Asi S=2,...,N. 
Analogous relations hold for the coefficients 6,5. If %(é),...,%2,(é) are 


the 2” characteristic exponents of system (3.2.8) we have 
a (6) = —eatm(f,) + Ole), aa(e)=—a+ oll) as e>0, 
aj (es) =7,+0(4) as e-0, a ee 


Thus, for all | ¢| sufficiently small, we have R[a;(e)] <0, 7=2,...,, and a (e) 
is real and corresponds to a simple multiplier. If m(f,,)=+-0, then R[a(e)|<0 
[> 0] for all |e| sufficiently small with em[f,,]>0 [<0]. If m[f,,]=0 then 


a (€) = €7(Soile.o) + O (é%). 


We have now, for h=1, 0.—0,= —a—0=—4, % = —4, B,-=0, 0,—01=7; 
(== 3, ms 20, and. DY (2.4.0). 


M= R (Soil.—o) = ER {204 Ay29 419 aa (a? + Pw?) (a cae 7 lo Dy 1)} a 


+3RE fr" ass odno +d 04 ++ Pw) (—4, Cy tlob,,)}. 


Since 
Cun = —a*(aiy,+ 0f4)), Dia= 0, 
Cine == Cresi= @ Ca. eto) ae 
Die oe Dy os1= aga ia = Ve eae 0 oe Ss SS 04, Sao. 


by trivial computations we have (3.2.7). 


If =1, then M is given in the first line of (3.2.7) and thus M <0, and a (e) <0 
for all | e| 4:0 sufficiently small. | 
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Remark. For n=1, system (3.2.8) becomes 


m= ea (—fi%y—hi%), %=—a%,+ea (Ai%1 + fir), 
and the contention of (3.2.iii) is a consequence of (3.1.iii), case (1a), with y= W3—1. 


Examples. The equation y’’+ay'+ef(t)y=0, a>0, has all solutions y(t) 
approaching zero as t—>-++ oo [with y’(t)]. Consider the system 


yitay, = éf,1(é) Vi t+ éfie(t) ¥.=0, 
Vo Co Vo + dy Vz é fay (t) Mit €feo(t)¥e=0, 


with m [7,14] =Cy19+-0. Then for | e| sufficiently small, ¢C,,)>0, all its solutions 
approach zero as t—>-+ oo. 


Consider now the system 
yi +, + ecost-y,+ ecostsy,=0, 
Yo + Vo + ¥,+ ecost-y,+ ecost-y,=0. 


We have here m[cos t] =0, M = —+4<0, and all its solutions approach zero as 
t—> + oo provided | ¢| is sufficiently small. 


Consider the system 
y, +4, + ecost-y,+ esint-y,=0, 
Yo + Vo + Yo + ecost-y,+ ecost-y,=0. 
We have here m[cost]=0, M=+4>0 and infinitely many of its solutions are 
unbounded in [0, +00). The same holds for the system 
vy, +4, + ecosi-y,+ ecost-y,=0, 
Vo + Yo + Yo — Sin? -y, + ecost-y,= 0, 
with m[cost]=0, M=4>0. 
Both (3.2.ii1) and (3.2.iii) are particular cases of the statements (3.2.v), (3.2.v1). 


We give below a few statements which show how the previous ones can be ex- 
tended. They are all corollaries of (2.4). 


(3.2.1v) Consider the system 


y FO W+ 2 fall K=O, 0; > 0, p=1,...,k, 
4=1 


yi + oy; +49, 4 eLifin) n= 0, ¢;, d>0, J=R+14,...." 
—- 


where 1k<n, 20;, 6,+0,, 0;—9,=0 (moda) 7 -=1),-1, W— 40 ay Wy & Ped, 

fin (t) veal periodic function of period T=22/w, L-integrable in [0, I], with Fourter 

series kes 

tin) ~ Da (ajn1cosl@t + bj,,;sinlwt),  4,h=1,...,n. 
I=0 


If M, <0, for allh=1,...,k, then for all | |= 0 sufficiently small all solutions of 
(3.2.9) approach zero as t>-+ cc. If M,>0 for at least one h=1,...,k, then for 
all | e|+=0 sufficiently small infinitely many solutions of (3.2.9) are unbounded in 
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[0, + 00). In both cases M, 1s given by 


lee) 


M, =); D (Page Cag Gj Drji) > i. ery 
j=1 1=0 
j+h 


where 


Pujt = 0, Vi = 24210 [40% oF (/? w a; o4)?]4, 


(21, 5k, (+h 


Pijo= Gn j0— 0, 


Prjo= — 246,[(d; — 04)? + 7, Gl Kujo 90, Mids ee 

Prji= — 42 6j(08j1— 4705) Mii» Mj = 2210 (Brjr — 4Onj1) Maj, 
My j1 = (Bi j1 ai (c} oan 4d,) O; He c (44 c ie d; aie on) 2 — (cj eM 4d;) or (2B, j;0— 65 )%, 
p= t2tcd—d+Po%, Py =—o,t4;+ Po?, 

vnj1=Bot+ 41 Pw+ of (444 —d), 

bij = UO +4 eRe +o (4td—d), gk 1,-0.0, b= 12.28 
Ci ji = yj) Ant + Dy 51 Dnt» Dyj1 = An Dn jt — Dining» 

Wed = Ag ancy hy t=O 1 

Onjo = 0, hf Ne 


(3.2.v) Consider the system 


nit 6 dhl) w=0, 
= 
(3.2.10) yi +07 y+ ed (5, (1) Paes Saas One eee 
st 
yeh YW +a + ED tial Ma =0, 6; dO, Gee ee 
i 


where 25k<n, 20;, o;+0,, 0,—0,+0 (moda), 7h, 7,h=2,...,k, & real, 
fn (t) 1s a real perrodic function of period T=22/w, L-integrable in [0, I] with 
Fourier series 


fin (t) ~ Dd (4jn1 Cos wt + b;,,;sin lot), A me Ole oe 
t=0 


Tf either a; 90, €4449<0, M, <0, h=2,...,R, OF a49=0, M,<0, h=1,2,...,, 
e+0, then for | e| sufficiently small all solutions of (3.2.10) approach zero ast—+>+ o., 
If either ay;9==0 and we have at least one of the inequalities ea,,,><0, M,> 0, 
h=2,...,R, OY a49=0 and we have at least one of the inequalities M,>0, h=1, 
2,...,k, then for |e|=- 0 sufficiently small infinitely many solutions of (3.2.10) 
are unbounded in [0, + ov). 
In both cases M,, is given by 


ae 2 22 (bust Cait + GnjxDa pi), eR Ra leven o 
j= - 
7 +h 
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where 

Pepe) =e, ak: Puen 4] SRE 1 hen; Gigs On (2a 
Priji= 2 aay = Ores Cig 0, ADRS Priji= (4; — Po) My1, 
hji= 6,40 M;1, M;,= 24 [l4wt + d? + Pw? (c# — 2d,)}>, 1h A, 
Pro — 276%, Tro= 0, h=2,...,83 Pr =44(— of + Po), 


Gi17— 0; h=2,...,R; Prjo= 9, Gnjo— 9, [he Ds sip he f 22 dere s 
Prjo= — 246, [(d; — of)? + of EG Opn IOs a Re = ke ae 


The formulas for j,;1, Cif Nay Dk IR 2, Fae hy beat; 2). averting 
same as those for p,;; and qy,;, of Theorem (3.2.iv) above if we put c,=0, dj= 0; ; 
(hae Pee 


(3.2.v1) Consider the system 


Ww +amteDhiOm=0,. a0, 
al 


(32244) y; + 07 9; + €>) fn() %=0 , o> 0, JD AS, Re 


H+ 659; + yy; + eX tall )¥n=0, 6, d>0, J=R+1,...,0 


where 2=kSm, 20;, o;+0,, 6;—0,220 (mod @), 7==h, 7, h=2,...,k, & real, 
fin(t) ts a real periodic function of period T=22/w, L-integrable in [0, T], with 
Fourier series 


Me 


fin (2) OG (a;,,coslwt + b;,;sinlat), 1, Wes TO 


Dy CUMCTEG rors OnE. Gyio > 0 0, AEDS. wh Orta eg 0, 0, WA, 25. ie 
e=£ 0, then for | e| sufficiently small all solutions of (3.2.11) approach zero as t+ . 
If ether ay49+:0 and we have at least one of the inequalities €a,,),>0, M,>0, 
h=2,...,k, OF A449=0 and we have at least one of the inequalities M,>0, h= 
1,2,...,k, then for |e|0 sufficiently small infinitely many solutions of (3.2.11) 
are unbounded in [0, + co). In both cases M, is given by 


=>, Ds PagtCajut M1 Daj | =a ea gonice 


ith —0 
Piuo=—2*, Ni10= 9, Peace: hjo= 9; ae eee 
Prjo= 4G, Gjo=0, F=R+14,....0 
Pui= — (2471 (V+ Po)4, 11=0, 


Oi (2a) 1G, Poy j= 0, = 2, ees, b= 1,2, 
prj = (2a)4 (4; —P ow?) m1,  ijr= (24) * log mz, 

m= [di + Pw? (ci — 24) + bot] ae he a gle lee le 2s as 
Dy19 = (24) a2 (a? + of) 1077,  %10= 9, 

Pa Oe pial! igs, tle 1, D5 58:5 
Pijo= hio=%, jae NESS a Rigel (icc aren 
Pajo= — 22 6,[(d; — 09)? + 04.65), Gio 10) ji een oy es 


a 
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the formulas for Prjts Wnjir 7 = 25-09% Wed ok Rs i=1,2, ..., are the same 2 
the formulas for pyjis Injr Of Theorem (3.2.1v) above if we put c;=0, d;=a;, 
OND Pe 

The proofs of (3.2.iv), (3.2.v), (3.2.vi) are omitted since they are obvious 
modifications of the previous ones. 


Examples. Consider the system 


yi + y, + 2ecost-y, + Ecost: 2+ Ecost-y3= 0 
(3,2:42) Vo + 2%, + ecost-y,+ Ecost- V+ Ecost:yz3= 0 
y3 +43 + y3+ ecost-y, + ecost- 2+ Ecost-y3=0 


where ¢=+-0 is real and is thought of as being sufficiently small. By (3.2.vi) 
with k=2, n=3, a1 9=0 we have M,=—0.01<0, M,= —0.052<0, and all 
solutions of (3.2.12) approach zero as t+ co. By replacing in (3.2.12) fis, 
or f;,, by —cost, we have M,=+0.01>0. By replacing in (3.2.12) fg, oF fe, 
by —5 cost we have M,=0.26>0. By replacing in (3.2.12) fig by —2sinz#, 
or fy, by +2sinz, we have M,=0.24>0. In all these cases infinitely many 
solutions are unbounded in [0, + co), no matter how small ¢>0 is chosen. 


Consider the system 


yy, -E 2y,-esin?-y, + ecost-y,= 
(3.2.13) . 
Yo + 242 + Yo + ECOST: Vy + ECOSE: Yo = 


where ¢=+0 is real and is thought of as being sufficiently small. By (3.2.1v) 
with k=1, n=2, we have M,= — 7;<0, and all solutions of (3.2.13) approach 
zero as t->+ oo. By replacing in (3.2.13) /,,=cos¢ by f,;=sin#, we have 
M,=+75>0. By replacing in (3.2.13) fig by —sint, we have M,=+ 44. In 
both cases infinitely many solutions are unbounded in [0, +-oc), no matter how 
small ¢ +0 is chosen. 
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Bounds for Solutions of Nonlinear Differential Systems 


MICHAEL GOLOMB 


Communicated by L. CESARI 


1. Introduction 
Bounds have been developed for solutions of differential systems of the form 


& = A) x +F(t,x(d) (1.1) 
by many authors ([3], [9], [10], [11], [12], [13]; for more complete bibliography 
see [4], [5a], [6]), usually based on the assumption that the nonlinear terms, 
represented by F(t, «(é)) in (1.1), decrease faster than the linear terms as || «|| >0 
and have the effect of small perturbations of the solution to the linear part. 
In this article bounds are established for solutions of systems (1.1) with various 
kinds of nonlinearities, some that increase more slowly than the linear terms 
when ||x||-> oo and others that decrease faster than the linear terms when 
|| ~||+0. In each case the bounds are derived from bounds to the general solution 
of the linear part assumed as known. In the case of (relatively) slowly increasing 
nonlinear terms the bounds hold for solutions with large initial values, in the 
case of (relatively) rapidly decreasing nonlinear terms the bounds are valid for 
solutions with small initial values. The bounds reflect at least the growth of the 
various terms involved and are rather sharp for some equations. They are 
developed from the solution to a certain nonlinear integral inequality [equa- 
tion (2.6)]. This solution [see (2.17) and (3.18)] is an extension of a result for 
the corresponding linear integral inequality, which has been used frequently for 
the purpose of estimating the growth of solutions of differential systems and was 
first formulated as an explicit lemma by R. Bettman ([J]; see also [4], p. 35, 
[7], [14]). 

2. Large Initial Values 

The systems considered are of the form 

a* = A(t) + F(t, «(f) (2.1) 
where x and F are vectors with » real components and A is a matrix with n? 
real components. The norms used are denoted as || x|| and || A|| and are defined 
as the sums of the absolute values of the components. The components of A(t) 
and F(t, x) are defined for all t=¢, and all real x, locally L-integrable functions 


of ¢ for fixed x and continuous in x for almost all #=t). Moreover, it is assumed 
that 


Fe,)\|Sa@ (lel), t2% (2.2) 
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- where a(?) is an integrable function, y(u) defined for w>0 as a positive non- 
decreasing continuous function which is submultiplicative, that is 

gp (wr) So) et). (2.3) 
These assumptions will be made throughout this article and will not be restated. 


For any real vector x) equation (2.1) has a solution x= x(¢) which exists 
on some interval to the right of t=‘) and assumes the value x) at t=f). This 
solution satisfies the integral equation 


AD) eal) Xn +1 XU DIAG ee) ar (2.4) 


where X(t, t) is the matrix that solves the equation 


dX 
7 =AOX 


and for which X(t, t) is the identity matrix; X(¢) is an abbreviation for X(E, to). 
Another assumption to be added to those stated above is 


IX@D|SeO(), StS! (2.5) 
where g(é), b(t) are locally L-integrable functions and 
WXOll<e@, OOS||X20|, wt. (2.5) 


Since X(é, t) = X(#) X1(z), inequality (2.5) is always satisfied for g(¢) =||X(4]||, 
b(t) =||X+ (||. However, sharper bounds are often found. For example, if the 
matrix A (i) in (2.1) is constant and has the m characteristic numbers A,, Ag, ..., A, 
then || X(é, t)||=||X(é—1)||Sc(«) expa(t—t) for any «>max Re 4; so that we 
may choose g(t) =c(«) exp(a?), b(f) =exp(—a#). However, in Sections 2 and 3 
we use g(t) and || X(é)|| interchangeably, with the understanding that where g (?) 
is identified with || X(4)|| there b() is identified with || X+(@]||. 


From (2.2), (2.4), (2.5) one derives the inequality 


0) <e(0 [f+ fH) 9 (Ho) dz], 2.6 


which is satisfied in a certain interval to the right of t=¢). The following ab- 
breviations are used here: 


AQ =|leO|],  fo= [loll (2.7a) 

h(t) = a(t) dé). (2.7b) 
If we put 

r(t) =f h(t) pf) at, (2.8) 


then by (2.6) and (2.3) 


 — h(t) (FO) Sh ele (fo +7) ShO o(eO) olfo+7)- (2-9) 
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Thus, if r’=dr/dt, y/ 
Se GO (2.10) 
P(fot+7) (8) 
Let p(x) be an integral of p(w) defined by 


va) =| AL hata HO (2.11) 


Then (x)>0 for x>1, y is continuous and strictly increasing. Let x=y*(y) 
denote the inverse function defined for y(0)<y<y(o). Because of the sub- 
multiplicative property (2.3), w(x) satisfies the inequality 


aty 1+y/% 


ve+s)—vlx)=f Say au= vz yl(i+2), wy>o. (2.12) 


p(x) “\ 


x 


Let inequality (2.10) be integrated between the limits fy and ¢. Since 7 (f) =0 
one obtains 


p(fot+7) —vlfo) Ss, (2.13) 


where the abbreviation 


t 


=f) p(g(t)) dt =f a(x) b(t) (|| X(a)||) ae (2.14) 


0 


is used. Using (2.12) in (2.13) one obtains 


v(t OY < Pe (), 


fo} fo 
and if 
PO) s() <y(oo) (2.15) 
then 
1+ 10 << ye(200) og), (2.16) 
and because of (2.6) 
f() Stov"( 2" so) e0. (2.17) 


For the solution «(¢) of (2.4) this means 
le OI Sloll v*(2 22? s@) XI (2.18) 


If condition (2.15) is satisfied then the right-hand term of (2.18) gives a finite 
upper bound for the solution x in the interval from f, to ¢. It then follows that 


the solution can be extended beyond ¢. Thus, the solution exists in any interval 
ty) StST for which 


lal (Chee es 
‘5 oo) el | . 


It should be noticed that if w/(w) increases indefinitely with w then T increases 


, that is the solution exists in an arbitrarily large interval 
ats ie if only xo is sufficiently large. 
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If the integral J dulp (wu) is divergent, the solution exists for all t=¢, and any 


aan Ane the Pel f)StST, the right-hand term of (2.18) represents an 


upper bound for bac )||. In particular if the integral Wie (uw) is divergent and 
the integral es - 
s(cs) = Fae) (>) (XC) 21 (2.29 


is convergent, then ||x (é)|| is bounded on fj) <t< oo by K\|X(t )||, where depends 
on ||xo||, but not on ¢t. Thus, if ||X(é)|| is bounded, that is, if the linear system 
x’—=Ax is stable, then every solution of (2.1) (with x)-+-0) is bounded, and if 
x'—= Ax is asymptotically stable, then every solution of (2.1) tends to 0 as t+ ow. 
These results are summarized in 


Theorem 2.1. The equation 


ax 
=A) x + Fb) 


has a solution x= x(t) with x (ty) = %9+0 in any interval ty StS<T for which 


ee ae 7 I|¥o\l du 
s(T) eee b(t) @ (\|X(z) ||) dt < aaa o (u) 

and an upper bound 1s given by 
le Sllzoll ye (Le se) ix@ll, sm SEST. 


If the integral Faull) is divergent, then the solution exists for all t=t, and any 


ep BOs See the integral s(co) 1s convergent, then ||x (t)||/||X(¢ )I| is bounded 
(on t=t,) and every solution x(t) 1s bounded if the system x'=Ax ts stable, and 
tends to 0 as too if x'=Ax ts asymptotically stable. 


For illustration, assume 


Aba = a0) elle oe OSs Baa aS (2.21) 


Then the integral fdu/p(u) =f u*—'du is divergent, and equation (2.1) has a 
i af 
solution x(t) existing for all ¢2¢) for any %)=% (fo) +0. Moreover 


x OllS (aol! + ks)" XO 


eae (2.22) 
if R+0 and 

Ol Sollee XO|,  ¢2% (2.23) 
ieee.) tere 


s(t) =f a(x) b(t) ||X(z)|P* at. 


If s(co)<co and the system x’= Ax is stable (asymptotically stable), then the 
solution x =0 of (2.1) is stable (asymptotically stable) if k=0, but need not be 
so if k+0. This latter occurrence is illustrated by the single equation 

dx 


roe —x+e'|x|!sgn x, t>0 (2.24) 
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whose solution for x(0) = x )=-0 is 
a(t) = e*(4 + |xo|? — e-#')? sgn x9. (2.25) 
It is bounded, but the solution x—0 is not stable. The general bound (2.22) 
gives in this special case 
[x(@)| Se*(1 + |x|? —e-*)’, 

which is the exact least upper bound. 

It was seen that if fdu/p(u) =o then the solution «= x(t; %») of (2.1) 

1 


exists for all £=¢), but does not necessarily tend to 0 as ||x9||~0. However 
x(t; %») remains bounded (for fixed #) as ||x||—0 provided p(w) satisfies addi- 
tional conditions. To prove this, it is sufficient, by (2.18), to show that 


y (x) = ryt(P) s}, x>0 (2.26) 


x 
is bounded near x =0, for any fixed s. This is obviously the case if 
p(x) S Ax, 0s4751 (2.27) 


where A is a constant. (2.26) is also bounded if 


p(x) SAX, OS "S14 
2.28 
ay Ve see See 
p(w) 


BS 0, weer Pads (2.29) 
For then 


or ny (2.30) 


which is bounded because of (2.29). 
Moreover, it is seen that lim yy (x) =0.if (2.27) is satisfied and also if (2.28), 
(2.29) with «+ B>1 are satisfied. Thus we have proved 


x 


: = 
Theorem 2.2. If either (a) x1q(x) is bounded near x =0 and | Saujp (| 
af 


. x —1 
ts bounded near x= co or (b) x *(x) is bounded near x =0, x6 fdulg (u)) 
ne 


1s bounded near x= co and B>0, «a+ B=1, then every solution x = x(t; x9) of (2.1) 
with X= % (to) 4-0 can be extended to the interval t = ty and \|x(t; x9)|| is a bounded 
junction of xy near x =0. Moreover, ||x(t; x») ||+0 as ||xo||->0 if either (a) is 
satisfied or (b) ts witha +B>1. 
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3. Small Initial Values 
The system considered in this section is again 
d 
= = A(t)x + F(t, x(d). (3.4) 


The same general assumptions on A, F are made as in Section 2; in particular, 
_ Inequalities (2.2), (2.3) and (2.5) are assumed to hold. In addition we assume 


a 


[a= (3.2) 


p (u) 


for any positive a. We consider solutions x= x(t) of (3.1) for which x)= x (t)) 
is sufficiently small. The symbols X(2), 0(4), /(, 4(, g(), A, 7, s(, oO, vO 
are to have the same meaning as in Section 2. 


We introduce numbers R, T as follows. If w(co) = du|p(u)<co, then we 


put R= oo, T=oo. Otherwise we choose, for given nae numbers R, T 
such that 


<7, oll yr(* Ue” o(7)) <R. (3.3) 


Let the function y(x), O<*<R, be defined as 


R 
70) =0, art] foe. O<«<R. (3.4) 


This function is continuous and strictly increasing, positive for 0<*<R, and 
the inverse function y*(y) is defined for OS y< co. By (2.16) and (3.3) 


loll +7 Sllaall v*(2 20s) <R,  Se<T. (3.5) 


Xo 


Thus the differential inequality (2.10) may be integrated between ¢, and ¢ with 
the use of the integral 1/y of 1/p. One obtains 

1x (\lxoll) — Ax (loll +70) S80, %St<T 
and . 

Pe AE ee = : 

“(lxoll +7@) S 14 — x(I|*oll) s (2) z tSt<T (3 6) 

provided 
yell) SQ <1. tSt<T. (3.7) 


If [dulp(u)< ow, then R= o and (3.7) is a condition on the smallness of |||]. 
i 


If fdulp(u) = 0 then x(x) +0 as R->oo and (3.7) is a condition on the size 
iI 
of R, no restriction on ||xo|| being necessary. 


From (3.6) and (2.6) we have 


le Ol S28 py) IO 


1— x((|*ol|) s 
Arch, Rational Mech. Anal. 19 


\3 Sue (3.8) 
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Since TI was arbitrary and since conditions (3.3) and (3.7) can be satisfied — 


with positive ||x9|| it follows that equation (3.1) has solutions ~(/)=-0 which 
exist in arbitrarily large intervals tj <¢<T and for which the estimate (3.8) is 
valid. These results are a ree in 


Theorem 3.1. If the integral { duly (uw) (a>0) is divergent, then the equation 
0 - 


ax 


a AW e+ F(t, x (t)) 


has a solution x= x(t) im the Se large interval ty<t<T, with the upper 
bound 


z 1('\ 0) Xt pia D 
leDSx(G— 2 ais@)IXO, — @St< 
provided x (ty) = %» ts so small that 
* d 
(|[X(z) ||) @ See 
n=fe () (2) (Xt) 4 < | ey 
Here x(x )=1/ fanjp (u), and R is any positive number 


fo.) 


), of R=co if [dulg@) < oo: 


1 


Fé +) ||Sa@|[x|**,  k>O0, (3.9) 


> ||xoll v*(H (I xoll) s( 


For example, if 


then equation (3.1) has a solution x= x(?) in 4) <¢<T, with the upper bound 


||%o|| 
IRAN seep erat] stl CaF (3.10) 
if ‘ 
Riltolik s(t) = lla J a(r) b(t) ||X(z)|P edt <1. (3.11) 


If the integral s(co)<oo then we can take ||x»||>>0 small enough so that 
the solution x= x(t; x9) exists and estimate (3.8) holds for all t=). Since 
%(*)—>0 as x-+0 it follows in this case that ||x(¢; %9)||/||X(@|| 0 as ||x|| +0 
uniformly for ¢=¢). Thus, in this case, x =0 is a stable (asymptotically stable) 
solution of (3.1) if the system «’= Ax is stable (asymptotically stable). Again 
assuming s(co)<co,a more precise statement on the order of oP | (£3 xo) ||/|| x (AI 


as ||x9||—>0 can be made for the case where Y (Up) <<Uy for sae Uy (necessarily 
O<u%<1). By (2.3) for OS*%<R, 0<yX1, 
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that is 

x(ny) 5 PY (AC OES RM es yeas (3.12) 
It follows that 

+(e )s4 O0<«x<R = 

ae Oa oo Wee acy = (3.13) 


The function (y)/y has by (2.3) a value =1 for y=1, and by assumption 
a value v9<1 for y=u,. Thus v=g(y)/y takes on every value between v, and 1 
as y varies from %, to 1. Some inverse function is defined for vy<Sv<1, and we 
have us y<1. Hence, the equation 


PO) — 1 — y (lll) 52°) G14) 


has a (not necessarily continuous) solution y= 1/a(||xo||) with © (|x|) <4/ 


Us 
if ||x|| is so small that 


x (oll) (20) S14 — y= 1 — Pio (3.15) 
Using (3.13) we then have 
%(\|%ol)) 
1 ei ray) = (llell) lll (3.16) 
and by (3.8) 
sup Ty 2 (lal) [ol 6.17) 


if ||x9|| satisfies condition (3.15). We thus have proved 


Theorem 3.2. If the integral f dulp(u) (a>0) ts divergent and the integral 
Ja(c) ) (|X (x)]|) de as pene then there exists a solution x= x(t; X») of 
GB. 1) on the whole interval t= ty for all sufficiently small ||x9\|, and 

sup HAE 0) O(1) as ||x|| 0. 
t>to 


If moreover p(uy)<Up for some Uy, then 


lle (3 ol =O(||xo|l) as ||xo|| > 0. 


SEP Lee) 


t>to 


In these cases x =0 is a stable (asymptotically stable) solution of (3.1) if x'=Ax 
is a stable (asymptotically stable) system. 


The following example may serve to illustrate this theorem. The system 


(3.18) 
= = = (sin logt + coslogt—2a)x%,+%4**, k>0, #>0 


19* 
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is used (with k =1) by O. PERRON ({12]; see also [4], p. 87) as a counter example 
for a proposed stability theorem. The solution of (3.18) is 


—at 


X14 = C1 @ 
.; Ci aeons (3.19) 
garg BE 8e Ca fe t sin log ax|, 
0 
where c,=%,(0), co=%2(0). Here 
\|X (4) || —_ Pesta V soaes Si) 
||X + (¢) || = ge eg tsinlogt+2at 
and the integral 
$ (00) a i (e-*7 + ag 3 cad a (7+ Pa sh Saka dt (3.20) 
0 
is convergent if either 
J St an ee (3.21) 
2 a—--zs 
or 
a=1 and k>1+—. (3.24) 


In either case, application of Theorem 3.2 reveals, and the explicit expression 
(3.19) confirms, that the solution x =0 of system (3.18) is asymptotically stable 
and 

|x, (t)| + | x2 ()| = C (Je,| ae | cz|) (e~9F deat halite 2 a 0 


for some constant C independent of ¢ and c,, c,. Although x=0 is not a stable 
solution of (3.18) if k=1 and a is sufficiently small (PERRON’S example), it is 
an asymptotically stable solution for any a> if only & is sufficiently large. 


4. Almost Constant A 


More specific results are obtained for the case where the matrix A(é) of the 
system 


&% = Alt)x + F(t, x) (4.1) 


is almost constant. By this we shall mean the following: There exist complex 
numbers 4,,A,,...,4, [the “asymptotic characteristic values” of A(t)] and 
matrices S(t, t), B(t, t) such that ||S(, 7)||, || S+(é 2)||, || B@ 7)|| are bounded 
for fj tSt<oo and the fundamental solution matrix X(f, t) of the system 
«%' = Ax is given by 


SX S+4= A(I + B) (4.2) 


where / is the diagonal matrix with the elements exp A, (¢—7) (R=1; 2.226%) 
and J is the unit matrix. Some conditions which assure that the matrix A(t) 
is almost constant according to this definition can be deduced from known 
results on asymptotic expansions of solutions to linear systems (see, ¢.g., [8], 
[71], [12]) and from a direct discussion (see, e.g., [2], [8], [4]). General conditions 
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have been given in connection with Crsari’s method of reduction of a linear 
system to “L-diagonal form’ (see, ¢.g., [5], [15], and for more references [5a]). 
From (4.2) one derives the bound 


IX@ Dc, grt, (43) 
where « is any real number > A= max Re A; (i =1, 2, ..., 2) and C is independent 
oft, t. Thus we may identify 6 (d), g(¢) of Sections 2, 3 with Ce~*', e*' respectively. 


If now 
IF (é, ») || S a (d) |] 


ie; betel he (4.4) 


and we choose (uw) =" then by (2.13) 
t 
s(t) =C f a(t) eM" dt (4.5) 
ta 


and s(oo)<oo if a(é) has an exponential type number <(1—A)A, that is 
a(t) exp [(k—1)A+ e]tis bounded on ¢ =f, for some ¢> 0. In this case Theorem 2.2 
applies if k<1 and Theorem 3.1 if k>1. The results are formulated in 


Theorem 4.1. Assume ||F(t, x)||Sa/(8) ||x||* (k=0), AQ) ts an almost constant 
matrix in the sense defined above with the asymptotic characteristic numbers d,, 
Ag, +++, A,, and a(t) has type number <(1—k)A where A=max Re d;. Then the 


system 
a 


7 =AM*e+Fltx()), t2% 


has a solution x= x(t; x9) existing for all t=ty for any x)= 0 tf RX1 and for all 
sufficiently small x) 1f k>1. In either case the solution has type number SA. 

In the case k=1 the theorem requires that a(#) has a negative type number. 
However, weaker conditions would suffice. Since this case has been treated 
extensively (see, ¢.g., [3], [4], [5], [6], [9]), it is not investigated any further here. 


This work was sponsored by the United States Army under Contract DA-11-022- 
ORD-2059. 
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Coefficient Problems 
in Systems of Partial Differential Equations 


ERWIN KREYSZIG 


Communicated by S. BERGMAN 


Dedicated to Professor Dr. H. BEHNKE on his 60" birthday 


1. Introduction 


Various coefficient problems have been considered for partial differential 
equations of the form 


1.1) 4 [49 + C(x, vy) v] = uz t+ c(z, 2*)u =0 


where z=x-+7y, z*=x—1y. That is, relations have been obtained connecting 


basic properties of a solution 
fo, 2) 


(1.2) (eee) —— ae 
m,n=0 

of (1.1) with properties of the coefficients w,,,,. These relations yield information 
about the domain of regularity, the location and nature of singularities, the 
growth of w(z, z*) and other properties. Cf. BERGMAN [1, 4, 5], INGERSOLL [6], 
KREyYSzIG [7]*. For this purpose certain Bergman operators (see below) serve as 
a powerful tool. These operators, which transform analytic functions f(z) of a 
complex variable into solutions of (1.1), preserve many properties of /(z), so that 
theorems on analytic functions yield corresponding theorems on solutions of (1.1). 
These Bergman operators have the property that all information about x (z, 2*) 
is obtained from the sequence (u,,9), m=0,1,..., of the coefficients in (1.2). 
A recent investigation [7] studied the extent to which information about u (z, 2*) 
can be obtained from other sequences of coefficients in (1.2). 

The unification of methods is one of the aims of the theory of linear partial 
differential equations. The consideration of coefficient problems is an approach 
in this direction. Indeed, various relations between properties of u(z, z*) and the 
coefficients in (1.2) are independent of the coefficients of the partial differential 
equation. It is therefore important to extend these considerations to systems of 
partial differential equations. 

In the present paper we shall investigate coefficient problems for solutions 


co 
(1.3) u(%, 21) 25 23) = > Umnpq 21 ees zp! 
m,n, p,g=0 


* Cf. the references at the end of this paper. 
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of systems of partial differential equations 
(1.4) 1 [Av t+ Cy (%q, Ve) 0] = Uapet + Cn (242%) ¥ =O, k=1,2, 


. = . 
Zp—=Xp+tVe, 2% =%,—1*Ve» 


é 1 ra) ae? a) 1 re) ~ oC 
————(—— ae —it.— AE — = }- 4 = , 
OZ, 2\OxX, OVR Oz; 2 Ow, oY 


U (21,24, 29,22) = 0(%, Ver Xa» Vo)» Gy (Ce ae) = 


Ble 
YY 
= 
— E>) 
R 
= 
<= 
= 
“~~ 


whose coefficients c,(z,, 24) are entire functions. 

BERGMAN & SCHIFFER [2] proved the existence of operators transforming pairs 
of analytic functions of two complex variables into solutions of the system (1.4). 
Coefficient problems for these solutions were considered by BERGMAN [3]. Those 
operators have the property that all information about u(z,, S529) Catia 
obtained from the sequences 


(Umopo)> WM; P= 0, Ny eign ADA Mesgng)s Mer =O, 4, cee, 


of the coefficients in (1.3). Hence arises the basic problem as to what extent 
information about the behavior of (2%, 27, 22, 22) can be obtained from other 
sequences of these coefficients, say, from the sequence 


(Herarba eter MP Dine pd ze One DORM, UOC ta 30t, PaO Nees 


This problem will be considered in the present paper. 


2. Bergman Operators 
We now consider some basic properties of the Bergman operators needed 
for what follows. For better understanding we mention also the case of a single 
partial differential equation (1.1). In this case, every solution regular at the origin 
can be represented in the form 


ad 
(2.1) u(z,2*) = f E(z,2*, 2d f(4z[4 —#)) 1-2)“ at. 
il 
This representation can be cast into the form 
(2.2) u(z,2*) = g(z) + d Qn(z,2*) J (2 —C)"-*e(@) de, 
n=1 0 
where 


Q, (2, 0) = 0, WA Ze 
and therefore 


(2.3) u(z,0) = g(2). 


The associated function g (z) of u(z, 2*) isan analytic function of a complex variable, 
regular at the origin. The functions Q, (z, *) depend on c(z, z*) in (1.1) but not 
on g(z). They are entire if c(z, z*) is entire. The relation (2.3) yields the starting 


point for considering coefficient problems. Cf. BERGMAN [1, 4, 5], INGERSOLL [6], 
KREYSZIG [7]. We note that | 


_ \2z_ dbg(2z) 
ig ces  et 


cf. (2.4) and (2.2). 
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When passing from a single partial differential equation to a system we 
encounter considerable complication. While in the case of one equation we apply 
the Bergman operator to analytic functions of one complex variable, in the case 
of a system we have to apply it to analytic functions of two complex variables, 
whose theory is not sufficiently developed for our purpose. In particular, little 
is known about relations between coefficients of the development of functions 
/(%, 2) and the location and type of singularities of these functions. Thus for 
studying systems of equations it is natural to investigate first solutions whose 
associated functions are so chosen that relations of the afore-mentioned type can 
be obtained. This will be done systematically in the subsequent sections. 

Each solution of the system (1.4) which is regular at the origin can be re- 


presented in the form 
aL 


Pete eas 25529) =| fz (4, 21,4) Ea (zo, 22, ts) {f: (4, [1— 4], 2.[4— 8) 


=il = 
+ felalt — 4], 22 [1 — #))} (1 — A)-4 (1 — &) Fd, dy. 
The associated functions f, and /, of u are analytic functions of two complex vari- 
ables, regular at the origin. The generating functions E, and E, are entire functions 
of the respective variables z,, z{ and z,, 23, and regular functions of #, and %, 
respectively, for |t,|<1, k=1,2. Furthermore, 


(2.5) Ee 0-u =) ENO, 4.) free 12 


E, and E, depend on c, and c, in (1.4) but not on the choice of f, and f,. Cf. BERG- 
’ MAN & SCHIFFER [2]. 

We note that if x,, y,, x,, and y, are real then z, and z; are conjugate; other- 
wise z, and z; are independent. 

From the complex solution u(z,, 2°, Z., 22) we may obtain another solution 

f (1.4) by applying to (2.4) the operator R defined as follows: 
co co Cc 
R( pb ann th 2h”) as / » Om etn > ann “tm 2h) . 
m,n=0 m,n=0 m,n=0 

The solution U=R(wu) of (1.4) is real for real values of the variables x,, y;, 
k =1, 2; we shall call such a solution a veal solution. Obviously, for real values 
of x,, y, the operator R is identical with the operator Re (“take the real part’’). 

From (2.4) BERGMAN [3] derived a representation suitable for considering 
coefficient problems of (z,, 2}, 22, 22); it has the property that the relation 
between the solutions of (1.4) and their associated functions becomes very simple. 
BERGMAN proved that for each system (1.4) there exist four functions S, (z,, 2%, C,), 
T,, (Zz, Ze, Cx), R=1, 2, such that each complex solution of (1.4), which is regular 
at the origin can be represented in the form 


(2.6) Wty a 5) 2a) = P(g,, 22) = &1(&, 22) + 80 (2, 22) 
lg if Sy 24 (C1, 22) 44 ae So £1 (2152) Ce 
0 


ach Ty &5 (Cx, 22) dty+ J Ty 89 (24,09) dee 


== fi ms So 81 (01,02) 40, 40 + (Ge Ty 89 (01,09) 40, dla, 
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where the associated functions gy (%, 2) and go (41, z#) of the solution w are analytic 
functions of two complex variables, regular at the origin. 


For a real solution of (1.4) we then have 


(2.7) U (21,21 22) 23) = R (P(g, 8s)) 
= 4{P (g, (21,22), 8e(%1» 22) +P (&, (et, AWE Co? Zy))}. 


In (2.6) and (2.7) the associated function g,(%, z}) can always be chosen so that 
(2.8) £2 (0, 22) = 82(4, 0) =0 
holds; cf. [2]. From (2.5) it follows that 

Sy (%4,0,6)) = Soe, 00) =0, Ta, 0,0) = LO, 22,02) =0. 


Using (2.8) we thus obtain from (2.6) the important relations 


(2.9) u(%,0, 29,0) = 1 (41,20) 
and 
(2.10) u(z, 0, 0, 22) =f, (2,22) + £1 (4, 0). 


In the case of a real solution we obtain from (2.7) and (2.8) 


(2.41) U(%, 0, 22,0) =3 (21 (21, 22) + 8 (0, 0)) 
and 
(2.12) U(%, 0, 0, 22) a 3 (g (4, 2) a0 81 (4, 0) = & (0, 22)) : 


Hence from the sequence (#959) OT (Umoog) Of the coefficients in the development 
(1.3) we may obtain information about basic properties of the complex solution u 
of (1.4), cf. BERGMAN [3]. The investigation is much the same as for a single 
partial differential equation. 

Clearly, we should be able to obtain information about the behavior of 
also from other sequences (tm 54), %>0, g>0, both fixed, m,p=0,1,.... We 
thus arrive at the problem stated at the end of the preceding section. To solve 
this problem it suffices to derive relations between an arbitrary fixed sequence 
(Umnpq), M>0, g>0, m, P=0,1,..., and the sequence (Umopo)- In doing this 
we may concentrate on complex solutions, since the generalization of the sub- 
sequent results to the case of real solutions will be immediate. 


3. The Functions uy, 9(2, 22) 


We wish to derive relations between an arbitrary given sequence (nine 
n>0, q>0, both fixed, and the sequence (#99) of the coefficients in the re- 
presentation (1.3). Using these relations and (2.9) we can obtain information 
about the associated function g, (2, 2.) of the corresponding solution w (z,, 2%, 2», 23) 
of (1.4) and thus about the afore mentioned and other properties of this solution. 
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In order derive such relations we introduce the functions 
co 
(3.1) WET, Aaliam DU ne tr ee, Tis ea = 1S We he 


The coefficients of each of these developments constitute just one of the sequences 
above. Relations between these functions are thus equivalent to the desired 
relations between the sequences. 


From (1.3) and (3.1) we have 


[o.e) 
(3.2) U(21; 21 Laas) a » Uy, q (21> 22) 4y fags 
n,q=0 


and in particular 
(3.3) Ugo (21, 22) = 81 (21, 22); 


cf. (2.9). We represent the coefficients of the system (1.4) in the form 


Cc 


ie) 
(3.4) Cy (4, zi) oo 3 ony, ay re Co (22, 23) =) oF) a 29 2%. 
m,n=0 ‘eine 


Introducing the functions 
(Bes) an ee (a) ssa (22) = Dep De SDM tocpabeGs==lOK tin ae, 


we may write the representation (3.4) in the form 
(3.6) Cr aees,) oe ) (2 Col(Znta tO. (zs). 22 2. 


We now insert (3.2) and (3. 2) into (1.4). In order that the two resulting equations 
be satisfied identically in 2, 23, the coefficient of each product zf” zx? in these 
equations must vanish. We thus obtain two systems of linear differential equations 
of the first order involving the functions u,,,(2,, 2.) and their first partial deri- 
vatives. Indeed, from the first of the equations (1.4) we find 


n—1 


ov 
ee y—1(21) % = 0; ND q=0,1,..., 


(3.7a) 


and from the second equation (1.4) 
ee 
7 ged SS = ee 
(3.7b) a +> 0, 1 Oy ly eaag =1Oe 


Hence our problem leads to the investigation of these two systems of differential 
equations. We denote the equations (3.7a) and (3.7b) corresponding to a pair 
(n, q) by [m, q], and [m, q],, respectively. 
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4. Singularities of Type S 

In the subsequent sections we shall derive relations between singularities of 
an arbitrary function ,,4(%, 22), 7 >0, g>0 [ef. (3.1) ], and those of the function 
Uo (21> 22) =81(%, 22) [ef (2.9)]. For simplicity we shall refer to an (arbitrary) 
point P: (z,, 2) in all the statements. It is clear how this is to be understood 
from the standpoint of the theory of analytic functions of two complex variables. 
We confine ourselves to some remarks in this direction which are necessary for 
what follows. 

Let t.(%, 22), ™>0, g>O [ef. (3.1)] be an arbitrary analytic function, 
defined in a simply-connected bounded domain D of 22,-space, containing the 
origin. Let P€®D be an arbitrary point and U(P)<® a sufficiently small neigh- 
borhood of P. 


Definition. A singularity of U,4(%, 22) at P is said to be of type S if u,.(%, 2) 
can be represented in U(P) in the form 


(4.1) Ung (21; Zo) = Nya (%1, 22) + Se (22) 1 (%) + $1 (%) 72 (22) 


where the functions 49,7, and 7 are regular at P while the functions s, and Sz 
are singular at P. 

Obviously, if u,,(%, 22) has a singularity of type S at P then its singularities 
in U(P) lie on portions of finitely many planes z= const. or zg= const. We 
should stress the fact that this is an essential restriction on the type and location 
of the singularities. Indeed, in the case of an arbitrary analytic function defined 
in ® and singular at P the singularities of this function in 1(P) may he on por- 
tions of arbitrary analytic surfaces. For example, if f(z,, 2.) is a meromorphic 
function in D it can be represented in the (sufficiently small) neighborhood U(P) 
in the form /[(%, 2.) =, (%, 22)/he(%, 2), Where h, and h, are regular functions 
in U(P) without common factor. Let hs =0 at P. Then if #,+0 at P, the func- 
tion / has a pole; if h,=0 at P, the function / has a non-essential singularity of 
the second kind at P; these singularities are isolated while poles lie on analytic sur- 
faces. For further details see BEHNKE & THULLEN [8] and BOCHNER & MARTIN [9]. 


5. Behavior of woo at Singular Points of Ung. Behavior of u, 4 at Regular 
Points of Uoo 

Relations between m9 (21, 22) = 8, (%, 22) and an arbitrary function u,, q (21> 29), 
n>0, g>0, may be investigated systematically as follows. We shall determine 
the behavior of #9 at singular points (Section 5) and at regular points (Section 7) 
of u,,, and, conversely, the behavior of u,,, at singular points (Section 6) and at 
regular points (Section 5) of uo. 

We start with the first of these four problems. The result will be that singul- 
arities of type S are the only ones of u,,, (21, 22), 17> 0, g>0, which may correspond 
to regular points of %9(z,, 2,). We shall prove the following simple basic 


Theorem 1. If a function u,4(21, 22), 0>0, g>0, has a singularity at a point P 
which is not of type S, then ugo(%, 22) is singular at P. 


Proof. Since the coefficients c,(z,, 2¢), k =4, 2, of the system (1.4) are assumed 
to be entire functions, all the coefficients c®(z,) and c(?) (z,) of the differential 
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equations (3.7) are entire functions. Let u,,(%, 2), >0, g>0, be an arbitrary 
function which has a singularity, not of type S, at a point P. We consider the 
sub-systems of (3.7) consisting of the 2ng+-n-+-q equations 


[v, Xa, Peat N OEM ne xe Aa ge 
[v, 2 ]p, Pan OM, ae 1), ne os 


At least one of the coefficients c™(z,) and c®)(z,) of each of the two equations 
[”, 7], and [m, g], cannot vanish identically, since otherwise u,, q(41» 22) would be 
constant. Since the singularity under consideration is not of type S, at least 
one of the two first partial derivatives of w,,,(%, 2) is singular at P. Hence, at 
least one of the functions u,,, 0S»<m or u,,,, 0<x<q, occurring in the respective 
equations [”, g], and [n, g], must be singular at P. The corresponding coefficient 
appears also in other equations of the sub-systems of (3.7) under consideration. 
By repeating the preceding conclusion we arrive at the intermediate result that 
at least one of the two following statements holds. (Case I) The function u,(%, 29) 
is singular at P. (Case II) The function wu, 9(z,, 2) is singular at P. In Case I 
the function 0%9,/0z, is also singular at P. In Case II the function 0u,,5/02z, is 
also singular at P. We consider Case I. From equation [0, g], we find that at 
least one of the functions u,,(%, 2), OS%<q, is singular at P. Using suitable 
equations from among [0, x],, 0<*<q, we obtain the final result that m9 (z, 2s) 
is singular at P. In Case II we have to consider the equations [v,0],, 0<vSn, 
and arrive at the same result. 


Conversely, let us assume that up 9(z,, 2) is regular at a point Q. Then, by 
integrating the equations of the above sub-systems of (3.7) we see that #,,,(%, 22) 
may be regular at Q or may have a singularity of type S at this point. Using 
(2.9) we thus have 


Theorem 2. I} the associated function g,(z%, 2%) of a solution u(z, 21, 29, 22) 
of (1.4) ts regular in a simply-connected bounded domain D which contains the 
origin, then each of the corresponding functions t,,4(21, 22), %,9=1, 2, ..., 18 regular 
or has at most singularities of type S in B. 


Hence the regularity of g,(z,,2,) in D imposes strong restrictions on the 
nature and distribution of possible singularities of each corresponding function 
Unq(%» %) in Q; the only singularity surfaces which u,,(%, %,) may have in D 
are finitely many portions of planes z,= const. or z,= const. If u,,(%, 2s) is 
singular at one point of such a plane, it is singular at all points of this plane 
which lie in 2. 


It is worthwhile to mention that in the case of a solution 
foe) co 

toy) u(2, 2*) = >) 4, (2) 2°") 4, (2) 2 Hine” 
n= m= 


of a single partial differential equation (1.1) the following similar but stronger 
result holds. If %(z)=g(z) [ef. (2.2)] is regular in a simply-connected bounded 
domain % which contains the origin, then the functions u,(z), w=1, 2,..., are 
regular in 8. This follows from (2.2) and (2.3), as was proved in [7]. 
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A similar improvement of Theorem 2 cannot be obtained in general but only 
for certain trivial cases where the functions c, (z,, z¢) in (4.4) are of such a special 
form that some of the functions #,,(%,%2), %>0, g>0, become completely 
independent of (1, 22)- 


6. Behavior of uy, at Singular Points of Ugo 
If u,,(%, %) has a singularity of type Sata point P, then (2, 22) is singular 
at P while some of the functions ,,(%, 22), O<v<m", O<x<q may be regular 
at P. For example, if the function (4, 2%) corresponding to a solution 
u(2,, 21, 2, 22) of the system 


(6.1) Uz, + 3 uU=0, pat, 2 


has a singularity not of type S at a point P, then the associated function g, (2, 22) 
of u is singular at P while u,,(%, 2) is regular at P. 


Hence we have the situation that singular points of m9 (z,, 22) may correspond 
to regular points of #,,(%, 22), #>0, g>0. 


For example, the following statement holds: 


Lemma 3. Suppose the series developments (3.4) of the coefficients cy(%, 24) 
and Cy (zy, 23) of the system (1.4) do not involve constant terms or terms which depend 
only on the respective variable z, and z,. Let z£” and 234 be the smallest powers 
of zi and zx, respectively, occurring in (3.4). Then, for every solution u(z,, 21, 29, 22) 
of (1.4), regular at the origin, the functions u,,, (2, 22), O<vSin, O<*SXq are entire 
functions. 

The proof follows from (3.7). 


The fact that singular points of u%99(%,, 2) may correspond to regular points 
of U,4(%, 22) is of principal importance. It has the consequence that the infor- 
mation about the domain of regularity and other properties of a solution u of 
(1.4) which can be obtained from a sequence (t#m,54), %>0, g>0, both fixed, 
is less than that which can be obtained from the sequence (#59). There are 
two different cases, as follows: 


(a) If for every solution of (4.4), regular at the origin, certain functions 
Uyy (2, 2), O<¥Sn, O<xXq, are entire, these functions must be completely 
independent of uo (%1, 22) = 8 (%, 22). Obviously this holds if all the coefficients 
of the differential equations (3.7) relating these functions and w%9(z,, 2) vanish 
identically. In order that all the functions w,, (z,, 25), Pe 0p Ase i2. besindes 
pendent of each other, the system (1.4) must reduce to the form 
(6.2) Uz, 2% =O, oe 

(b) Certain singularities of some special associated functions £1 (21, 2) May 
coincide with zeros (of sufficiently large order) of the coefficients of (3.7). Because 
of this coincidence, a function w,,,(%, Z), #>0, g>0 may be regular at such a 
point. However, since each coefficient of (3.7) depends only on one of the two 


variables 2, or z, this may happen only for singularities which lie on singular 
planes z,= const. or z,= const. 
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We may sum up our result as follows: 


Theorem 4. Suppose that the function Uy q (21,2), M>0, g>0, both fixed, and 
the function Ugo(%1, 22) ave related by some of the equations 


[v, 2],, Pale ke MP OF, ack a: 
and 


cf. (3.7). Let Ugo (%, 25) be singular at a point P so that at least one of the singularity 
surfaces of Ugo(%, %) passing through P is not a plane z,= const. or z= const. 
Then U,,9(%, 2) ts singular at P. 


7. Behavior of Ugo at Regular Points of Ung 


We consider an arbitrary function w,,, (21, 2%), ~>0, g>0, both fixed. Suppose 
this function is regular in a simply-connected bounded domain © containing the 
origin. We ask for information about the behavior of wu 9(z,, 2) in D which 
results from this assumption. If w,,, (2, 29) and u99(2,, 2.) are independent (cf. the 
previous section), we cannot draw any such conclusion. In order to exclude this 
trivial case we assume that u,, (21, 2) and u99(%, 2) are related by some of the 
Equallonsn |), 5° P41, Deen My sO, 1y ee! , gs andy) zip P= 0,4) Cane, 
%~=1,2,...,¢q, ef. (3.7). Ifa function c®(z,) or c®(z,) does not vanish identically 
it always occurs in several of these equations. This fact has the following con- 
sequence, The equations which relate u%9(%, 2) and u,,,(%, 2) always contain 
at least one of two special classes (A) or (6), and the equations of that class 
relate u%9(%, 2) and u,,(%, 2). Class (A) consists in the systems 


(7.1) (A,) a oe 1(%) CNG, = 0, Mea AG 2) saa Wy 
1 
and 
OUo x 
(7.2) (Axle 2x =e =) as (22) Upp = 9, WRECEAN (OME Oe 


Class (B) consists in the systems 


ee 
73) (By) =e 2.5) dal if ay) Ug = 0; i Ne see 
and 

a 
(7.4) (B,) Sins 45) oh cl Says eos 


Whether class (A) or class (B) relates the functions u9(%, 22) and u,4(%, 2s) 
under consideration depends on the form of the coefficients c,(z,, z@) of (1.4). 
We may consider class (A) only, since for class (B) the manner of reasoning 
is exactly the same. 

We eliminate from the system (A,) the functions u,,, v=1,2,...,%—1, and 
their derivatives. For this purpose we consider (A,) as a system of m homo- 
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geneous equations in the 27 —1 unknowns 


OUy g OUgg OUn—1,9 d1 
— . 2q? as po ees Un—1,q> a > 7 an , 
Oz, O24 4 


where the coefficient of the ‘‘unknown”’ 1 involves all the terms containing the 
functions “%,, %,q and their partial derivatives with respect to %. If we dif- 
ferentiate each equation of (A,) (7—1) times with respect to z,, we obtain a 
system (A,) which consists in »? equations in the same number of unknowns; 
the determinant of this system will be denoted by D(A,). If D(A,) =0, the further 
procedure of elimination has to be applied to a suitable sub-system of (A,) for 
which the determinant of the corresponding enlarged system (obtained by 
suitable differentiation) does not vanish identically. It is not necessary to pay 
further attention to the exceptional case D(A,)=0, and we may assume that 
D(A,) =0. The equation 


(7.5) D(A;) =0 


is a linear differential equation involving only the functions m,(%, 22) and 
Unq(%4,%) and their partial derivatives with respect to z,. We shall write this 
equation explicitly. Let [1,1], [2,1], ..., [w, 1] denote the equations of (A,) and 
[v, @] the equation obtained by differentiating equation [y,1] o—1 times. We 
may write the enlarged system (Aj) in the form 


O tty g v—1 sl Rata cy an PP thee 
B 


6 » ONS a 
(7 ) [y 0] ye dzg— p- i az8 0, 


We arrange each of these m? equations so that the terms involving the functions 
Uy, and u,, and their derivatives occur first. If we develop the determinant 
D(A,) by its first column then (7.5) takes the form 


7.7 ¥ (ante) tee + O(a) Ree) = 0; 


s=0 OR 


in this representation, 


v,0=1 — 


s 


aia OS s=<n), 0 (S220) 
W Di (Oss 1), 


D,, is the cofactor in D(A,) of the element in the first column and in the row 
corresponding to equation [o, t], and 


0 (s= oe), 
Pros= Oe 1\ a oS 
| s da OS ae 


is the coefficient of the function Puy,/0z1 in the equation [y, o]. 


a 
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For example, in the case of the functions Uzq(%, 2), (7.7) has the form 


(c6” — of? Reema A) u =o) eit) Mog 4 5 aa Oleg 5 a) Cag 
dz, dz, ot Oz, BV 22 


= 0. 


We may derive some general properties of the coefficients of (7.7) as follows. 
If all the coefficients b, vanish identically, the same holds for the coefficients a,, 
since otherwise “), would be the same function for solutions u of (1.4) corre- 
sponding to any function w,,, whatsoever. That is, D(A,) =0, but this case was 
excluded. If the functions a, vanish identically but at least one function }, 
does not vanish identically, then (7.7) is a differential equation in u,, which 
does not involve u,,. That is, #,,, and uw), are independent ; this leads to a contra- 
diction to the above assumptions. Hence, a,=:0 for at least one value of 
s=0,1,...,%—1, and also 0,30 for at least one value of s=1,2,.<,n. If 
s=m (<n) is the largest value for which a, 0, then (7.7), considered as a dif- 
ferential equation in %,, has the order m. 

In a similar manner we can transform the system (A,) into a single linear 
differential equation 


(7.8) > [est 


The preceding general properties hold also for the coefficients of this equation. 


w, (29) Ze Ate (ae) ugg 
Oe 023 


Now, the solution w(z) of an ordinary linear non-homogeneous differential 
equation L(w) =/(z) with analytic coefficients and f(z) analytic is regular in the 
intersection of the domains of regularity of these coefficients and /(z). Since the 
coefficients of (7.7) and (7.8) are entire functions of the respective variable we 
obtain 


Theorem 5. Suppose the coefficients c,(z,, z¢) of (1.4) ave such that, for a 
solution u, the corresponding function U,4(2, 22), M>0, g>0, both arbitrary and 
fixed, is related to the corresponding function Ugg(%, Z) by some of the equations 
Bey 1 Der) ee Ol Ae... g and [Nix], PHO, 1, ancy PO dy eg 
[cf. (3.7) ]. Then the following statement holds. If u,,(%, 22) 1s regular in a simply- 
connected bounded domain D containing the origin, then Ugo(%, 2) 18 regular in a 
sub-domain D' of D which differs from D by at most fimitely many portions of 
the planes z= const. and z,= const. 

This means that if ,, and m9 are not independent and wu, is regular in a 
domain ® of the above type, the function “)) may have some singularities in 
}, but only singularities of a very simple type. These singularities, considered 
in the intersection of D anda plane z,= const. or z,= const., are either isolated 
or each point of D in this plane is a singular point of U9. 
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